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ABSTRACT 


This  report  is  a  manual  for  using  the  two  computer  programs: 

1.  "Unbalance  Response  of  a  Rotor  in  Fluid  Film  Bearings" 

2.  "The  Stability  of  a  Rotor  in  Fluid  Film  Bearings" 

The  report  gives  the  analysis  on  which  the  programs  are  based,  and  the  in¬ 
structions  for  preparing  the  computer  input  and  for  interpreting  the  computer 
output . 
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SYMBOLS 


A  Cross  sectional  area  of  shaft,  in2 

Ajk  Major  and  minor  axis  of  ellipse,  in  (or:  lbs,  lbs. in) 

Influence  coefficients  for  shaft  section,  see  Eqs.(4)  to  (7) 

C  Radial  bearing  clearance,  inch 

C Bearing  damping  coefficient  for  translator^  whirl,  lbs. sec/in 

Bearing  damping  coefficients  for  conical  whirl,  lbs . in. sec/radian 
E  Youngs  modulus,  lbs. in2 


CL  Rotor  mass  eccentricity,  inch 

F*,R,  x-and  y-components  of  bearing  reaction,  lbi . 

I  Croaa-aectlonal  moment  of  inertia  of  shaft, 

(tti  between  atations  n  and  (nfl)Jt  in^ 

Tp  Polar  maas  moment-of  inertia  of  a  rotor  mass,  lbs. in. sec2 

(in  input:  lbs. in  ) 

Iy  Transverse  mass  moment  of  Inertia  of  a  rotor  mass,  lbs. in. sec2 

(in  input:  lbs.ln^) 

Bearing  spring  coefficients  for  translatory  whirl,  lbs/ln. 

Kr  Rotor  stiffness,  lbs/ln 


Rotor  stiffness,  lbs/ln 
Ln  Length  of  shaft  section  between  station  n  and  (ntl) ,  inch 

L  Bearing  length,  inch 

i  Rotor  length,  inch 

M  Bending  moment  (M  to  the  left,  M'  to  the  right  of  station  n) 

n  n  J 

x  and  y-components  of  pedestal  mesa,  lbs.a«c2/ln  (ln  input:  1 

Mr  Total  rotor  masa,  lbs.sec2/in 

Bearing  spring  coefficients  for  conical  whirl,  lbs. in/ radian 

2 

My  Maas  at  rotor  station  n,  lba.sac  /in  (in  input:  lbs) 


Rotor  length,  inch 

Bending  moment  (M  to  tha  left,  M'  to  the  right  of  station  nllbsln. 

n  n  * 

2 

x  and  y-components  of  pedestal  mass,  lbs. sec  /in  (in  input:  lbs). 

2 

Total  rotor  mass,  lbs. sac  /in 


x  and  y-components  for  force  transadtted  to  base,  lba. 
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Time,  seconds 
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Cosine  and  sine-components  of  unbalance,  lbs. sec  (in  input:  oz .  ir.) 
Shear  force  (V  to  the  right  of  atatlon  n),lbs. 

Bearing  reaction,  lbs. 

Components  of  the  rotor  amplitude,  inch  (in  output:  mils) 

Coordinate  along  length  of  rotor,  inch. 

Phase  angle  between  amplitude  radius  vector  and  unbalance 
see  Fig.  5,  radians 

Angle  between  major  axis  and  x*axis,  see  Fig.  5,  radians 
Pedestal  damping  coefficients  for  conical  whirl,  lbs. in. sec/rad. 
3*  j  3%.  j  components  of  the  slope  of  the  deflected  rotor, rad. 
Pedestal  spring  coefficients,  lbs/in. 

Pedestal  damping  coefficients,  lbs. sec/in. 

Angular  speed  of  rotor,  radlans/sec. 

Critical  rotor  speed,  radlans/sec. 


Subscripts  and  Superscripts 


X 

**i**M*iM 

c 

s 

n 

F 

O' 


in  x-dlrectlon 
in  y-dlrection 

first  index  gives  force  direction, second  index  gives' amplitude 

direction. 

cosine  component 

sine  component  ~ — 

applies  to  station  n 

pedestal 

relative  between  Journal  and  pedestal 
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INTRODUCTION 


A  rotor  supported  in  fluid  film  journal  bearings  is  a  complex  dynamical  syst 
which  exhibits  a  variety  of  physical  characteristics:  critical  speeds,  insta 
ility,  unbalance  vibrations,  etc.  In  designing  a  rotor-bearing  system  for  a 
given  application  it  is  necessary  to  have  methods  available  from  which  these 
performance  characteristics  of  the  system  can  be  predicted  and  thereby  ensur 
that  the  design  is  adequate  for  the  specified  operational  conditions.  It  it 
the  purpose  of  the  present  report  to  describe  two  computer  programs  by  which 
a  particular  rotor-bearing  system  can  be  investigated.  The  first  program: 
"Unbalance  Response  of  a  Rotor  in  Fluid  Film  Bearings"  calculates  the  whirl 
amplitudes  induced  by  a  specified  unbalance.  The  second  program:  "The  Stabl 
of  a  Rotor  in  Fluid  Film  Bearings",  calculates  the  threshold  of  instability 
for  the  rotor-bearing  system. 

In  the  dynamics  of  a  rotor-bearing  system  the  fluid  film  journal  bearings  pi 
a  very  importarit  role.  They  are  normally  the  predominant  source  of  damping 
such  that  without  this  source  it  would  be  Impossible  to  run  the  rotor  throug 
any  of  its  critical  speeds.  Secondly,  the  bearing  film  is  flexible  and  thet 
it  may  lower  the  critical  speeds  drastically.  The  film  flexibility  also  cau 
the  bearing  to  act  as  a  vibration  Isolator,  attenuating  the  dynamical  forci 
transmitted  to  the  pedestals.  Finally,  if  the  speed  gets  sufficiently  high 
the  bearing  film  looses  its  ability  to  dampen  out  any  transient  motion  of  tl 
rotor  and  transfers  instead  energy  from  the  rotation  of  the  rotor  into  a  wfc 
lng  motion  of  the  rotor  mass.  This  is  called  fractional  frequency  whirl 
("oil  whip")  and  is  a  self-excited  instability  of  the  rotor  bearing  system, 
speed  at  the  onset  of  the  instability  is  called  the  threshold  speed  and  as  f 
rotor  speed  is  Increased  beyond  the  threshold  speed  the  whirl  amplitude  lnci 
rapidly,  preventing  further  operation  of  the  machine.  It  is,  therefore,  net 
at  the  design  stage  to  ensure  that  the  selected  rotor-bearing  sy*t?ia  does  nc 
experience  Instability  in  the  operating  speed  range.  Likewise,  it  is  also  < 
sirable  to  evaluate  the  magnitudf  of  the  rotor  amplitude  due  to  a  residual  t 
balance  such  that  too  large  amplitudes  will  not  be  encountered  in  the  actual 
machine.  The  two  computer  programs  described  in  the  present  report  provide 
method  for  performing  these  calculations. 

1 


The  unbalance  response  program  la  very  general.  Ic  calculates  the  rotor  whirl 
amplitude  and  the  force  transmitted  to  the  base  due  to  a  given  rotor  unbalance. 
The  rotor  is  flexible  and  may  have  any  arbitrary  geometry.  Also,  there  can  be 
splined  couplings  in  the  rotor  and  several  bearings.  The  bearing  pedestals  can 
bi  assigned  both  flexibility  and  damping.  Since  the  bearing  film  forces  are 
not  the  same  in  all  directions  the  whirl  motion  of  the  rotor  is  treated  as  two- 
dimensional  such  that  it  becomes  an  orbit  around  the  equilibrium  position.  The 
orbit  is  elliptical  and  its  dimensions  and  orientation  vary  along  the  length  of 
the  rotor.  The  computer  program  calculates  the  whirl  orbits  for  a  number  of 
points  along  the  rotor  and  gives  also  the  components  of  the  force  transmitted 
to  the  foundations. 

The  program  for  investigating  the  stability  of  the  rotor-bearing  system  applies 
to  an  arbitrary  rotor  geometry.  There  may  be  several  bearings  and  the  stiffness 
and  damping  of  the  bearing  pedestals  can  be  included.  The  program  calculates 
the  speed  at  onset  of  Instability  (the  threshold  speed)  and  the  corresponding 
whirl  frequency. 

In  both  programs,  the  dynamic  properties  of  a  fluid  film  bearing  are  expressed 
in  terms  of  8  coefficients:  4  spring  coefficients  and  4  damping  coefficients. 

The  coefficients  depend  on  the  bearing  type,  the  bearing  dlmenalons ,  the  vis¬ 
cosity  of  the  lubricant,  the  bearing  load  and  the  rotor  speed.  The  values  of 
the  coefficients  are  given  in  a  previous  report  (Ref.  6)  for  a  wida  variety 
of  bearing  types,  geometries  and  operating  conditions. 

The  report  sets  forth  the  analyses  on  which  the  computer  programs  are  based. 
Detailed  instructions  are  given  for  preparing  the  computer  input  and  for  inter¬ 
preting  the'  output.  . 
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DISCUSSION 


aj-JS&aasal 

The  computer  programs  determine  the  interaction  between  the  bearings  and  the 
rotor.  The  unbalance  response  program  is  concerned  with  the  synchronous 
amplitude  of  the  system  under  the  action  of  unbalance  forces  and  the  stability 
program  is  concerned  with  the  free,  self-excited  motion  at  the  onset  of  hydro- 
dynamic  instability  ("oil  whip")*  Whereas  the  dynamic  properties  of  the 
bearings  derive  from  lubrication  theory  (Ref.  4) ,  the  analysis  of  the  rotor  it¬ 
self  derives  in  its  principle  from  the  Myklestad-Prohl  method  (Refs. I,  2,  3). 
However,  in  its  original  form,  the  Myklestad-Prchl  method  is  set  up  only  for 
calculating  the  critical  speeds  of  the  rotor  and  is  further  limited  to  plane 
vibrations.  In  t**e  present  analysis  the  motion  is  treated  as  two-dimensional , 
damping  is  included  in  the  bearings  in  addition  to  stiffness  and  the  analysis 
is  valid  for  any  spee^,  not  just  the  critical  speed. 

In  general  a  rotor's  cross-sectional  dimensions  and  its  mass  distribution  variei 
along  the  length  of  the  rotor.  Thus,  for  calculation  purposes  it  is  conven¬ 
ient  to  break  the  rotor  up  into  short  sections,  each  section  having  a  constant 
cross-section.  Furthermore,  when  there  ere  many  sections,  the  mass  of  each 
section  can  be  divided  into  two  parts  and  lumped  at  the  end  pointa  of  the  sectlc 
Concentrated  masses  like  wheels,  impellers,  etc.,  can  be  made  to  coincide  with 
an  end  point  of  a  section.  In  this  way  the  rotor  is  replaced  by  an  Idealized 
model  constating  of  a  number  of  mass  points  connected  by  weightless,  flexible 
bers.  The  model  can  be  brought  as  close  to  the  actual  rotor  as  desired  by  mak¬ 
ing  the  subdivisions  small  but  in  practice  only  a  limited  number  of  divisions 
is  needed  to  obtain  a  very  good  accuracy. 

Since  the  bearing  film  properties  to  a  large  extent  control  the  whirl  motion 
and  the  stability  of  the  rotor,  it  is  necessary  to  represent  the  dynamical  bear¬ 
ing  film  forces  as  accurately  as  possible.  The  method  of  representation  is 
based  on  the  aesumption  that  the  whirl  amplitude  is  small  compered  to  the  bear¬ 
ing  clearance  such  that  tha  dynamical  forces  can  be  replaced  by  their  gradients 
around  the  steady  stgte  journal  center  position.  Xn  this  way  the  dynasdcal 
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forces  become  proportional  to  the  whirl  amplitude  and  to  the  corresponding 
velocity,  and  the  factors  of  proportionality  are  called  spring  and  damping 
coefficients.  They  differ  from  conventional  mechanical  sprlng-dashpot  systems 
by  also  containing  cross-coupling  terms  in  addition  to  direct-coupling  terms,  i.e., 
the  dynamical  force  in  a  given  direction  (say  the  x-dlrectlon)  is  not  only  pro¬ 
portional  to  the  amplitude  ar.d  velocity  components  in  that  direction  but  is  also 
proportional  to  the  amplitude  and  velocity  components  in  the  mutually  perpen¬ 
dicular  direction  (i.e.,  the  y-dlrection) .  Hence,  in  an  arbitrary  reference 
coordinate  system  with  x  and  y-axds  the  two  dynamical  force  components  can  be 
expressed  by: 

Rr  =  “Kor*  “CfwX  *  ^  “* 

Rf  *  ~  “*  Ctpj 

where  x  and  y  are  the  amplitude  components,  x  and  y  ere  the  velocity  components 

K  and  K  are  the  direct  coupling  spring  coefficients,  C  end  C  are  the 
xx  yy  r  o  r  »  xx  yy 

direct-coupling  damping  coefficients,  K  end  K  era  the  cross-coupling  spring 

xy  yx 

coefficients,  and  end  C^x  ere  the  cross -coupling  damping  coefficients.  These 
8  coefficients  ere  functions  of  the  bearing  Sommerfeld  number  defined  through 
the  rotor  speed,  the  steady  state  bearing  reactions,  the  lubricant  viscosity  and 
the  bearing  dimensions  (for  gas  bearings  the  coefficients  are  functions  of  the 
compressibility  number,  the  bearing  eccentricity  ratio  and  the  whirl  frequency). 
Thus,  the  coefficients  vary  with  speed.  A  method  for  calculating  the  coefficients 
is  given  in  Refs.  4  and  5  and  values  of  the  coefficients  for  several  bearing 
types  nay  be  found  in  Raf.  6. 

Frequently  the  pedestals,  on  which  the  bearings  are  mounted,  are  as  flexible  as 
the  bearing  film.  In  such  cases,  the  pedestal  stiffness  must  be  Included  In  the 
calculations.  For  completeness  the  analysis  allows  for  both  stiffness,  damping 
and  inertia  In  the  pedestals.  Furthermore,  as  the  rotor  bends  under  the  in¬ 
fluence  of  the  unbalance  forces,  the  journals  become  cocked  in  their  bearings. 

The  fluid  film  resists  the  tilting  and  this  can  be  expressed  by  a  set  of  8  spring 
and  damping  coefficients  in  analogy  to  the  previously  discussed  coefficients. 

The  unbalance  response  analysis  includes  this  effect,  both  la  the  bearings  and  In 
the  pedestals.  The  resistance  to  tilt  normally  effects  the  rotor  motion  only  at 
speeds  above  the  second  or  third  critical  speed  but  If  the  pedestals  ere  made 
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soft  for  alignment  purposes  resonance  conditions  may  exist  which  can  only  be 
explored  if  the  effect  of  tilt  is  included.  For  the  stability  analysis  this 
effect  is  in  almost  all  cases  very  small  and  it  has  been  ignored. 

Occasionally  the  rotor  is  not  a  single  member  but  consists  of  several  rotors 
connected  by  splined  couplings  (e.g.,  a  turbine-generator  set  connected  by 
a  splinec  coupling).  The  unbalance  response  program  allows  for  including 
splined  couplings  anywhere  in  the  rotor  and  assumes  that  no  bending  ooacnt  is 
transferred  through  the  coupling.  This  feature  is  not  included  in  the  stability 
program. 

In  the  unbalance  response  program  the  whirling  motion  of  the  rotor  is  generated 
by  unbalances  built  into  the  rotor.  In  general  the  unbalance  varies  in  magnitude 
and  circumferential  location  along  the  rotor  such  that  under  speed  the  unbalance 
forces  may  bend  the  rotor  into  complicated  shapes  (e.g.,  resembling  a  "cork-screw") 
The  bend  rotor  whirls  around  its  steady  state  position  (i.e.,  the  position  the 
rotor  would  occupy  if  there  were  no  unbalance  forces)  with  each  point  of  the  rotor 
axis  describing  an  elliptical  path.  The  dimensions  and  orientation  of  the  ellipse 
varies  along  the  length  of  the  rotor. 

If  the  rotor  runs  at  high  speed  and  has  large  disks  (e.g.,  turbine  wheels,  etc.) 
mounted  on  the  shaft  the  gyroscopic  moment  becomes  important,  especially  if  a 
wheel  is  overhung  at  one  end  of  the  rotor.  The  gyroscopic  moment  is  proportional 
to  the  mass-moment  of  inertia  of  the  wheel,  the  square  of  the  speed  and  the  de¬ 
flection  angle  of  the  rotor.  If  the  rotor  motion  is  considered  as  a  transverse 
vibration  of  a  beam  (i.e. ,  the  whirl  orbit  is  a  straight  line)  the  gyroscopic 
moment  tends  to  "soften"  the  rotor  and  lower  the  critical  speed.  On  the  other 
hand,  if  the  bearing  spring  and  damping  coefficients  are  the  same  in  the  vertical 
and  the  horizontal  direction  the  rotor  whirl  orbit  becomes  a  circle  and  the  gyro¬ 
scopic  moment  stiffens  the  rotor.  Actually,  the  whirl  orbit  is  elliptical,  i.e., 
somewhere  between  a  straight  line  and  a  circle,  and  the  effect  of  the  gyroscopic 
moment  csn  only  be  assessed  by  performing  the  complete  rotor  analysis.  Zt  is  a 
non-linear  effect  since  it  depends  on  the  dimensions  of  the  elliptical  whirl  orbit. 
In  the  present  analysis  the  gyroscopic  moment  is  taken  into  account  in  the  un¬ 
balance  response  program  and  is  calculated  by  an  iteration  procedure. 
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The  greatest  difficulty  encountered  In  performing  the  numerical  calculations 
Is  the  magnitude  of  the  numbers  and  the  loss  of  significant  figures.  These 
difficulties  become  pronounced  when:  a)  there  is  an  excessive  number  of  rotor 
mass  stations,  b)  the  rotor  la  very  stiff  and,  c)  the  bearings  are  very  stiff. 
There  is  no  universal  remedy  for  the  problem  but  If  trouble  arises  two  possi¬ 
bilities  may  be  tried:  a)  reduce  the  number  of  rotor  stations  to  the  essential 
minimum  and,  b)  apply  a  scale  factor. 

Let  the  scale  factor  be  0l  .  Then: 
multiply  the  speed  by  o(. 

multiply  (El)  by  e(X  (e  g. ,  multiply  E  by  d*  ) 
multiply  the  bearing  spring  and  damping  coefficients  by  <& 

(i.e.,  multiply  K  ,  wC  ,  M  ,  cuD  etc.  by  d*  ) 
multiply  the  pedestal  stiffness  by  of  and  the  pedestal 
damping  coefficients  by  A. 

(i.e.,  multiply  Xj,  and  Xtj  by  d*  ,  6^  and  by  st  ) 
leave  the  rotor  masses,  the  rotor  length,  the  pedestal  masaea 
and  the  unbalance  unchanged. 

The  numerical  results  will  give  the  amplitude  unchanged  whereas  the  bending 
moment  and  the  transmitted  force  must  be  divided  by  dl  to  obtain  the  actual 
values. 


•  and  Biamiflnlm-IaMtlgM 

Referring  to  the  sign  convention  given  in  Fig.  2  and  considering  first  a  con¬ 
tinuous  rotor  the  three  basic  equations  for  determining  the  rotor  motion  are: 

it,  4*  t 


(1-a)  Force  balance  for  a  shaft  inc 
(2-a)  Moment  balance  for  a  shaft  lncri 

(3-a)  Shaft  deflection: 

where: 


M*  ET 


x  -  amplitude  in  vertical  direction,  inch 
y  -  amplitude  in  borlsontal  direction,  inch 


/ 


X 

a 

A 

I 

E 

P 

<w 

0) 

M 

V 


coordinate  along  the  rotor  length.  Inch 

eccentricity  between  mass  center  and  shaft  center,  inch 

2 

cross-sectional  area  of  shaft,  in 

4 

cross-section  moment  of  inertia,  in 

2 

Youngs  modulus,  lbs/in 

2  4 

mass  density,  lbs. sec  /in 

mass  moment  of  inertia  per  unit  length,  which  is  effective  in 

2 

gyroscopic  moment,  lbs. sec 
angular  speed,  radians/sec 
bending  moment ,  lbs . in 
shear  force,  lbs. 


For  the  stability  analysia  set  a  ■  0  and  redefine  u>  to  mean  the  whirl  frequenc) 
These  three  equations  may  be  combined  to  give  the  familiar  4th-order  differen¬ 
tial  equation  governing  the  unbalance  vibrations  of  a  rotor: 

<4-«)  -fax  (£1  fe)  =  fAu»'(x+*)+«>%£ 

and  the  same  for  the  y-dlrection. 

(see  Ref.  3,  page  330) 

For  a  circular  whirl  orbit: 

For  a  straight  line  orbit: 

For  an  elliptical  whirl  orbit,  aee  Eq.(28)  and  (29)  in  this  report. 


At  the  bearings  there  is  an  abrupt  change  in  the  shear  force  and  the  bending 
moment  due  to  the  baaring  reactions.  Let  e  bearing  be  at  a  ■  •  Then: 

(5-a)  Km**~VX**;  =  + 

(6-a)  -M^r  S « 
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where  K  ,  C  ,  H  ,  D  etc.,  are  the  bearing  spring  and  damping  coefficients. 

XX  XX  XX  XX 

Actually,  the  effect  of  the  pedestal  should  be  included  in  the  above  equations 
as  shwon  InEq.  (12)  and  (13)  in  the  analysis. 

The  numerical  method  uses  Eqs .  (1-a) ,  (2-a)  and  (3-a)  by  rewriting  them  into 


finite  difference  form: 


AV=  CoVjAaz)  (x+a) 

aM=  v«  wft, -£,.)«](£) 

ax--  (&)«  +  jfyfcfc 


Together  with  Eqs.  (3-a)  and  (6-a)  these  equations  form  a  act  of  recurrence 
relationships  which  can  be  solved  step  by  step,  starting  froai  one  end  of  the 
rotor  until  reaching  the  other  end.  The  details  ara  given  later. 

Occasionally  it  is  desired  to  perform  a  dimensionless  analysis.  The  two 
governing  quantities  are:  ' 

<*-•>  *  rf- 


(8-a) 


where: 


tHj  Mr 


reference  value  of  IX,  lbs. in 

rotor  span  between  bearings,  inch 

/*  PAdx,  total  rotor  mass,  lbs. sac2/ in 
o 

rotor  stiffness,  lbs/ln 

equal  to  or  proportional  to  a  critical  rotor  spaed,  rad/sac. 


For  a  uniform  shaft  (El  ■  constant,  A  ■  Constant): 


where  n  designates  the  order  of  the  critical  speed.  Thus,  for  the  first  mode: 
n  ■  1  i.e. , 


(ei),=  ^pVr=?.-)47-(-Er 


(Uniform  shaft,  first  mode) 


However,  It  is  not  necessary  that  be  a  critical  speed  but  Eq.  (7-a)  must 
be  satisfied. 


The  dimensionless  parameters  become: 

x'  =  x/a0 

z/t 

(El)'*  £I/(GI)0 

V'~  vu.kr 

M‘ '  M/^  If  ^ 

K./K--  <%)(&)(&) 
UJ--  {&)  <Hf.)  («g-) 

K*  m„/k ;/•*  (<*)  (&)(*)*($?) 


where: 


(fiA)  •  JtfA/Mr 


reference  value  for  the  rotor  mass  eccentricity,  inch 

reference  value  for  the  radial  bearing  clearance,  inch 

V 

actual  radial  bearing  clearance,  inch 
reference  value  for  the  bearing  reaction,  lbs. 
actual  bearing  reaction,  lbs. 
bearing  length,  inch 


The  dimensionless  bearing  coefficients  are  given  the  form  above  since  the  values 
obtained  from  lubrication  theory  are  CK^/W,  QnC^/W,  etc.  Normally,  a  dimen¬ 
sionless  analysis  is  only  performed  for  a  simple  system  where  all  bearings  are 

identical,  i.e.  C  «  C  and  W  ■  W  .  In  that  case  the  basic  dimensionless 

oo 

parameters  are: 

speed  ratio:  ) 

dimensionless  rotor  stiffness:  K^,  ■  CK^/W 

dimensionless  bearing  coefficients:  CK^/W.CmC^/W ,  ^(CM^_/WL^) «tc. 


Thus,  Co  perform  a  dimensionless  calculation  for  a  given  value  of  K'^use  aa 
Input  to  the  unbalance  reaponae  computer  program: 


Speed  -  (—)/.  1047 1976 


0) 

«  f 

Mas*  at  station  1 


n  3.86069*10  (m*atation  weight,  lba; 
m . 


L  “J 


n»nuaber  of  stations) 
5 


(I  -IJ 

(Ip-I^)  at  station  i  ■  — “2  3.86069*10 

**I 


Crosa-sectlonal  moment  of  inertia  for  section  i-(i+l):  1000*1/1^ 
Young  modulus  ■  1 


Length  of  section  i-(i+l)  •  t.ll 


CK 


Bearing  spring  coefficient  ■  (*r-**) 

r 

.  OuC 

Bearing  damping  coefficient  •  “T  ( — rj^) 


Unbalance  such  that:  £U  (os. in)  ■  6177.1  EU  (or. in)  -  6177.1 

i  *  i  y 


Then  the  computer  output  will  give: 
amplitude  **  and 


bending  moment  ■  M'  ■  M/a^^i  ■ 

a 

transmitted  force  *•  (actual  force) /•0*r  ■  (actual  force)/  WK^ 


COMPUTER  PROGRAM:  UNBALANCE  RESPONSE  OF  A  ROTOR  IN  FLUID  FILM  BEARINGS 


This  section  of  the  report  describes  the  basic  analysis  and  the  detailed  in¬ 
structions  for  using  the  computer  program:  PN0011:  "Unbalance  Response  of  a 
Rotor  in  Fluid  Film  Journal  Bearings"  for  the  IBM  704  digital  computer.  The 
program  calculates  the  rotor  deflection  and  bending  moment,  the  pedestal  de¬ 
flection  and  the  transmitted  force  resulting  from  a  specified  rotor  unbalance. 

It  differs  from  conventional  programs  by  taking  into  account  the  variation  of 
support  flexibility  and  damping  along  the  whirl  path  of  the  rotor. 

The  supports  for  the  rotor  consist  of  a  fluid  film  bearing  on  a  pedestal,  both 
members  possessing  flexibility  and  damping  for  translatory  and  rotational  motion 
The  flexibility  and  damping  are  linear  in  displacement  and  velocity  respectively 
the  proportionality  factors  denoted  as  spring  and  damping  coefficients.  The 
fluid  film  is  represented  by  4  spring  coefficients  and  4  damping  coefficients 
for  translatory  motion  and  similarly  for  rotational  motion,  thus  allowing  for 
coupling  between  the  motion  in  two  mutually  perpendicular  directions.  The 
pedestal  has  no  such  coupling  and  is  represented  by  2  spring  and  2  dmsplng 
coefficients  for  both  translatory  and  rotational  motion  with  corresponding 
pedestal  mass  and  mass  moment  of  Inertia.  Hence,  each  point  of  the  rotor  will 
whirl  in  an  el’iptlc  path  around  its  steady  state  position. 

In  addition,  the  program  Includes  the  effect  of  gyroscopic  moment  and  provides 
for  couplings  in  the  rotor. 

TjffiQRBHSAL^HALYg;? 

The  analysis  is  an  extension  of  the  Kyklestad-Prohl  method,  see  Ref.  1,2  and 
3.  The  rotor,  which  is  actually  a  continuous  system  with  an  infinite  number 
of  degrees  of  freedom,  is’  replaced  by  a  finite  number  of  lumped  masses  conn¬ 
ected  by  weightless  springs.  The  computer  program  calculates  the  vibrational 
response  of  this  equivalent  system  exactly. 
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{ 


Thus  the  accuracy  of  the  results  depends  only  on  how  closely  the 
idealized  system  resembles  the  actual  rotor. 


Starting  from  the  left  end  of  the  rotor,  the  program  calculates  step 
by  step  the  bending  moment,  shear  force, -  slope  and  deflection  along 
the  rotor.  Neglecting  the  shear  force  contribution  to  the  deflection, 
we  get  from  Fig.  1: 

“>  Mr,,.  •  Af»  ♦  i,  V, 

<J>  0M  *  Or,  *  Cl*ML 

m  Xrl,l  *  X,  +  U  Or,  +  c „  Mr!  ♦  dr,  Vn 


where: 


The  program  asaumes  El  constant  between  mass  points.  At  the  mass 
polnta,  the  forces  acting  on  the  rotor  are  introduced.  Four  contrlb- 
utlons  exist:  (1)  inertia  force,  (2)  unbalance  forces,  (3)  bearing 
reaction,  and  (4)  gyroscopic  moment.  Zn  general,  not  all  4  contrib¬ 
utions  apply  to  each  mass  point. 


Inertia  force.  The  rotor  performs  harmonic  vibrations  at  the  same 
frequency  as  the  rotational  speed.  Thus  the  inertia  fores  is: 


V 
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(8) 


-  m 


(9) 


dzx 

<nl 


muJlx 


-m  ~L ,  - 
<nl  1 


Unbalance  forces.  To  allow  for  change  In  circumferential  position  of 
the  unbalance  along  the  rotor,  the  unbalance  is  given  two  components 
Ut  and  Uy  .  This  gives  rise  to  an  X  and  y  component  of  the  unbalance 
force: 


(10)  (v„-  Vi^,)uilb  =  0JzU%ca*oit 

(11)  (Vy.rt-  =  oJ  COSUJt  +VZUx3 

Bearing  reaction.  The  bearing  supports  have  flexibility  and  damping 
for  both  translatory  and  rotational  motion  of  the  rotor.  Since  the 
equations  for  the  two  types  of  motion  are  analogous,  only  the  equat¬ 
ions  for  translatory  motion  will  be  derived. 


The  bearing  support  is  shown  in  Fig.  3.  It  consists  of  a  pedestal 
with  mas 8  (  ) .  supported  by  springs  (tta,4y)  *nd  dashpots 

(<la,<fj).  There  is  no  coupling  between  the  X  and  tj  direction,  i.e. 
no  transfer  impedance,  nor  between  the  translatory  and  rotational 
motion.  The  pedestal  supports  the  bearing  fluid  film  which  is  rep¬ 
resented  by  4  springs  and  4  damping  coefficients.  If  the  relative 
motion  between  the  journal  center  and  the  bearing  housing  Is  denoted 
(X)  W*),  then  the  bearing  reaction  becomes: 


t V"  "  Kyx'-Cy tX  -  tffj'-.Cft  j 1 

Setting:  .  .  .  , 

X  »  Xc  COiUit  +  StW  art 


»  "if 


c  COScJt  ' 


sin  u)t. 


we  get  from  Newton's  second  law  for  the  pedestal  mass: 
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(Kxx  *  i\ *“ )xc  *  Oi((ax  +$x)Xs  +  Kxy  Ljc  +  U)Cxy(jlim(i>x~UJLM4x))(ct uJCiXs 

~u(C> or  +  (Tx)Xc  +()(u*i\i~lJtM*x)x>s-U)Cxijljl  f  ^  -  "UXJir 

KijxXc +  +w(Cj^f<Ti|)  ij4  * 

-  U)  Cijx Xg  +KijxX 5  <jc 

Solving  the  equations  we  obtain: 

CVw  “  (“AVwXc  "A  \4jfX<  ~A  Vex  ljc  ~AV*x  co*  cut 

(14)  +  (^ViJtXe-AV'*jiX*^^\4t^c‘‘^'4r^s)5«»»lJt 

AVb|  cj*)  cosu/t 

♦  (  *  V^Xe-MfyXj  +A  Y^  -  a  Vf  «jt)  smut 

where: 

AVajt  *  K«cf  +  OJ  Car^  4  K^Jty  £  4  <*->Ca^r 

*“  ^(jCnj  | 

AVfcjf  "  Kxx li  +  U>C«ft  +  S  ♦  ItlCwjt 

(15)  ^Vig  — K«  i  +  w  C*li  -Kaj  +  u)Ca|$ 

aV«^  »  KytK^u^Cvpitl^s^  k)Cff  tr 

AV^  »  -  Klj*i+'  coCf.  k  ~  fcfftr  +u>C<n|  • 

A  Vtij  *  K«fir  J  +U)  Cyi  J  ♦  t^f|  fg)  C^|  ^ 

A\4l^  *  ♦(*)  Cifif  -  <fuC«f|^ 
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and: 


The  equations  for  rotational  motion  are  analogous  to  eq. (14) except  for 
a  sign  reversal  (sign  convention,  see  Fig.  2): 

mtq*  (* M<ur0e  +  f$)  CO&Uft 

+  ♦  4  Me*  (fa )  si  nut 

(16) 

(aM^0c  +  ^5  )  CO**lt 
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where  the  coefficients  etc.  are  computed  from  eq.  (15) 

as  AtAurtV**  s  iMs*  •  AVk*  etc.  by  replacing  the  Cranalatory 
•pring  and  damping  coefficients  by  the  corresponding  rotational  co¬ 
efficients. 

Since  the  fluid  fils  coefficients  are  functions  of  speed,  directly 
through  the  SoaBsrfeld  number  and  Indirectly  through  the  dacraase  of 
eccentricity  ratio  with  increasing  speed,  the  computer  program  prov¬ 
ides  for  expressing  the  coefficients  as  a  function  of  spaed,  e.g. 

(18)  Kx*  »  Kjur  jO  ♦  K*x,i  *  U)  ♦  K«,|  '  ft)1 

and  similarly  for  the  other  coefficients.  (J  is  the  rotor  speed  in 
radians/sec. 

Gyroscopic  Moment.  The  gyroscopic  moment  derives  from  the  change  of 
the  angular  momentum  vector  of  the  rotating  rotor  mass  as  it  whirls 
in  an  elliptical  path  around  the  steady  state  position  of  the  rotor. 

For  two  special  cases  the  gyroscopic  moment  is  known: 

circular  whirl  path:  *  (.Ir“Xr)uXS 

(19) 

straight  line  (transverse  vibrations): 

where  $  la  the  slope  of  the  rotor  deflection  an dXf  and Xr  are  the 
polar  and  transverse  mass  moments  of  inertia,  for  an  elliptical  path 
the  gyroscopic  moment  is  no  longer  linear  with  rasps ct  to  the  slope  of 
the  rotor,  indicating  that  an  elliptical  path  is  actually  not  possible. 
However,  in  general  the  effect  of  the  gyroscopic  nenont  is  not  too 
big  and  for  the  present  analysis  an  elliptical  path  will  ho  assumed. 

The  coordinate  system  is  shown  in  Fig.  4,  where  9  is  the  steady  state 
shaft  center  position  and  0*  is  the  whirling  shaft  center.  The  moving 
coordinate  system  (frfif  )  is  defined  by  its  unit  restore: 
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(20) 


_  ( e/ieW 

et  s  [  foe1*?1 


=  (  \j  9*1  fx 


* 


♦ 


i 


f/^F 

•fTTFTp 

9 


{imp 


ins**  * 


\~(-JL-  -jLrr,  jfZIXpl) 
*  4»9l7fi-j  ys5^5  »i fe^P  r*  v) 


*  o) 

TiJw)  *  l6’*'1) 


The  angular  velocity  vector  becomea: 


(2D  CJ*  ( 


The  moment  needed  to  sustain  the  motion  is  given  by  Eulers  equations: 

Mj  *  rrO>|  *  (l,p~  IT)  UifUJ] 

(22)  Mn  *  JT  +  ^0>| 


where  I  denotes  mass  moment  of  Inertia  and  Xj  **  Xy.  ?  *nd  Tj*Xp 


Let  us  first  assume  that  (X,^  )  corresponds  to  the  directions  of  the 
major  and  minor  axis  in  the  elliptical  variation  of  the  rotor  slope. 
Their, 


(23) 


6,  •  E  cos  (u>t+&) 
<f,  •  6  91*  (u*t**) 


Combining  eq.  (20),  (21)  and  (22): 

M,  •  «**  *  i,u‘ti(ies  J 

which  clearly  shows  that  the  gyroscopic  moment  is  not  linear  with 
respect  to  the  rotor  slope.  However,  -only  the  first  harmonic  can  do 
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work  on  Che  rotor.  Hence  a  Fourier  analyaia  will  be  performed.  The 
following  integral*  apply: 


$fn/  cofA  fa. 
E^Cos1/  +6*™** 


0 


sin*  cos*  dX  _  n 
+6*si*\ r)1 


iff 


ur\lx. 


E'cos*X+&,SirtiX 


.iff 


SitilL  4* 


_ 4* _ 

£xcoiV  +  6Y*^ 


di _  . 

(f  *fS**+&*li*\)* 


.IS... 
a  (e+6) 


a? 


M L 
r6 

e‘6‘ 


Then  the  firac  harmonic  becomea: 

-*»-(f4I'''rr)wV. 

(24) 

M«j  *  (-fTS  Zp'Ir)  ufx9i 

In  the  limit,  eqa.  (24)  agree  with  eqe.  (19). 

Eqa.  (2  4)  muat  be  tranefomed  back  to  the  actual  (^<f) -coordinate 
ay a tern.  Setting 

0  •  Bt.  C&s  uii  *  Bi  tim  tut 

(25) 

m  fc  t#«  cut  t  <ft  ft*  uJt 
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describing  an  elliptical  variation  of  slope,  we  get: 


where  £  is  the  angle  from  the  position  X-axis  to  the  major  axis  E,  pos¬ 
ition  in  the  same  direction  as  U)  .  Then: 


(27) 


$,  «  9  cos  p  *<f> 

< f,  m  -  0  +  (j>  COS  p 


Substituting  4q.  (26) - (27)  into  eqs.  (24)  gives: 

(28)  M,  -  (m,!,-  tofaMtglp- 9cTr]casuft 

(29)  -  Mr  *  ( MyJi  ceswt  *  o;*[2AMA*Iip*f.jJ s^uit 

where 


(30)  a  Ad  - 

(31)  A  M+a  * 

(31)  *MdY«  (*♦*)■♦<*-*)« 


Since  eq.  (28)  and  eq.  (29)  are  not  linear,  an  iterative  method  ia 
used.  For  each  rotor  speed,  the  program  performs  a  number  of  iterat¬ 
ions.  The  first  iteration  is  done  without  gyroscopic  amment.  After 
the  first  iteration,  the  gyroscopic  moment  is  calculated  from  eq.  (28)- 
(29)  and  these  values  are  used  in  the  second  iteration  and  so  on.  The 
calculation  has  converged  when  the  relative  change  in  rotor  amplitude 
and  slope  between  two  iterations  is  smaller >han  a  specified  limit. 
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EQCAT ION'S  FOR  ROTOK  CALC'- LAI  ION 

The  bending  moment,  the  ahear  force,  the  slope  and  the  deflection  are 
expreased  by: 

M/  ■  Mjrc  cos  wt  +  Mr,  ainwt 
V*  =  Vue  cos  out  *  Vx$  s/iujfc 

B  *  6c  cos  vat*-  Os  tmuit 

X  »  Xc  cos  w-t  ♦  X  t  siwuJfc 


and  similarly  for  the  y •direction.  Then  eq.  (1),  (2),  (3),  (8),  (9), 

(10),  (11),  (14),  (16),  (17),  and  (29)  say  be  coablned  to  give  the 
equations  used  in  the  rotor  calculation  (see  Fig.  2): 

M«n  *  Mitch  ^sMunBcn^^MuriBpi*  dMxa4ti» 

Mini)  *  Mart  "  A 

a* 

Mtqndca+A  +A May 4Li * t~ (M ^cm - ^ 

»  Mq«a  -A  Mdq»  &.*  ♦A  Mci^^Ja'AMfcyi,fwtAMai^<f>,l|+(MqfH'trt<j5«)Air. 
VgtM  ■  Vxc(eM+[#l*u>  *d\iis*] )fc«"A V^aXea'AVcaa^u.-AVaBi^M  ♦  u»2U*« 

(32)  Vwe  *  *AVmm  Xc*  ♦  VaeaJXje  VAa^rt'il/trt^rt  -0JlUy» 

Vq«*  "Vqqan  *  oVuf*  Xa*  -  A  V^y,*  ♦  -A  Vfc^  y*, 

Vqea  •!***•  ♦**&!•*••  -dtfnt*  X*  ^Uan 

Ms^ihT  Ma«a  ♦  L*V** 

Mjca«**  Ma*e  ♦  Ln  Vrta 


Myt  M  >j<M  4LnV«j«* 

*  LuVtjin 


8c*.* 

den  +  An  Mkh  bn  Vwa 

0.n 

♦  An  Mru  ♦  bn  V*»n 

4  CU\  M  qcn  ♦  bnVpn 

f*' 

>4  AnMtpa  4  bnVt^fM 

X*,*rl  * 

+  ♦’bn Mren  td^Vscn 

* 

** 

+  Lnd«n  4bnM<m  +cUVm* 

y«,nn* 

^cn 

4  U  4  LnMqt*  +flln 

4*^1* 

+  Ln^H+  UkMqsn  4«lnUpn 

In  th«  above  aquations  da,  fet,  d*  are  (Ivan  by  sq.  (4),  (3)  and  (7), 

- dA'Ufn  and  a^aim,  - dVdT  by  aqa.  (15)  and 

(  Mm*m  M*n  )&  — - )(^f#  by  sq.  (28)-(29). 

Boundary  Conditions.  The  rotor  la  assuend  to  have  free  soda.  Mo  loss  In 
gsnarallty  occurs  by  this  condition  sines  it  nay  be  changed  by  letting 
the  end  points  have  bearing  support.  A  proper  choice  of  support  coef¬ 
ficients  will  then  allow  for  any  type  of  end  conditions. 

For  a  rotor  with  free  ends  the  bending  scent  end  the  sheer  force  are 
sero  et  the  end: 

(33)  M»ti  •  Mui  •  Mti*  Mhm*  KifU  *M|m  *  0 

(34)  VWr*V(|4r"V(fer  *0 
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Starting  from  the  left  end  of  the  rotor  (aee  Fig.  2),  eq.  (33)  Is  uaed.  - 
However,  the  slope  and  the  deflection  ere  unknown.  Using  the  super- 
position  principle,  each  unknown  le  applied  separately.  A  summation 
gives  the  combined  effect.  Ten  calculations  are  performed,  using  eqs. 
(32). 


1. 

Ba  •  1 

Bu  *i«*i*< * U'SUtTO 

2. 

Bs>  •« 

Bn  •  fn  *  A,  *Xw  =  Xti  —  4ei  =^fjt =  U*n=  Ui ifQ 

3. 

<Pc,  •  1 

Ou  mf*i  * Xai ■  )(f i  =  ^4  =  =  UiM=  1(^=0 

4. 

<fb.  - ' 

9ci  *Xn*Xll 

5. 

tfc*  1 

Be  i  =^ii  =(ft  i  =^u=y*i 

6. 

tfv*  I 

Bc\  =&v  =$i  =Xti  =t(‘i  =4ii 

7. 

1 

5ti s  Bw  — “<01  =X«  -Xm=4*i 

a. 

V* 1 

9. 

Ujr*  S  Ur*  (jftp,  Vl/q* 

6n  -Bit  =<Pii  — X»i  =ifu  =lfn  —  0 

10. 

Gyroscopic  svment  applied 

6k,  =  A,  =Um=V° 

For  each  calculation  eqe.  (32)  are  used  to  calculate  the  bending  mom¬ 
ent,  the  shear  force,  the  elope  and  the  deflection  along  the  rotor.  At 
the  right  rotor  end,  station  r,  eq.  (34)  must  be  satisfied,  l.e. 


Maopt - 

Mmr,i 

a, 

Misti- - 

Mkmt* 

A, 

— Mi»f  f 

. 

»o 

(35)  J 

MyiS,  x - 

►  4 

*  =*  J 

-Wqss'.l 

- 

U*,  * 

& 

—  y*/*,  1  -  Vy«v,/» 

!Mt> 

H *r,\ - 

¥tK» 

Ml 

-WrpF,  f -Vyi/,1# 

Yu*,' 

|^«5t - 

y* 

"1(4*',  1  -^qiy|(e 

1 

Ki«r.» 

V. 

- 

■w 

.J 

k.  j 

-K|sr,1 

Eqs.  (33)  are  then  solved  for  Ai ,  A* — — and  the  actual  values  of 
bending  moamnt,  sheer  force  etc,  along  the  rotor  can  be  determined. 

A*-  a  given  rotor  speed,  eqs.  (33)  are  first  solved  without  gyroscopic 
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moment,  i.e.  Mxcr^o*  a\^4rlifOi  Then  the  gyroscopic  moment  is 

applied  according  to  eq.  (28)-<29)  and  new  values  are  found  for 

V‘<  from  eqs.  (35).  This  process  is  repeated  until  at 
the  k' th  iteration: 


(36)  - +1^%::  i  , 


i0^vift!oMrK - hri 


<6< 


where  is  the  convergence  limit  specified  by  the  computer  input. 

If  the  calculation  does  not  converge  within  a  specified  number  of  iter¬ 
ations,  the  program  goes  on  to  a  new  rotor  speed. 

In  the  computer  output,  the  rotor  deflection  is  given  by  the  dimensions 
of  the  elliptical  whirl  path.  We  have: 

X*  Xe  COSUlt  +  Xa  Si’iMut 

(37) 

l 

U«  ifc  tosujt  +  l^sinwt 

As  shown  in  Fig.  5,  the  (X,y )  -coordinate  system  is  rotated  an  angle  (3 
in  the  same  direction  as  u)  to  become  Than 


X,  -  Clccj  (ut+a) 
-  lo  si*(iJt  td.) 


where  a  and  b  are  the  major  and  minor  axis  respectively  of  the  ellipse. 
From  Fig.  5: 

X,  *  xco*p  simp 
*  — '  *  -XSiHp+Lf  Cosp 

Then: 

(39)  (Xc.  *X»  Hjc+tft)  * 
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Here  1C  Is  necessary  Co  allow  b  Co  become  negative.  The  reason  la  Chet  the  trans 
formation  from  Che  x-y-coordinates  Co  the  ellipse  oust  be  able  to  discern  be¬ 
tween  forward  end  backward  whirl  (l.e.  the  shaft  center  may  travel  in  the  seae 
direction  or  in  the  opposite  direction  of  the  direction  of  rotation  depending 

x 


on  Che  values  of  x  , 

c  s 

the  instantaneous  radius  vector  be 


y  end  y  ).  Let  ths  angle  between  the  x-exis  erd 
C  8 


Then: 


l.e. 


Therefore: 

(39a) 


y : 

Y  -iAtr'd) 

J  *  ‘ 


(Xglft-Xstyt )  ^0  •*  hrwtrd 

(xti wWI 

*  itrefjit  I  ms 


To  find  d  end  expand  Eq.(38) 

(a)  acoid  •  XcUSj 

(b)  -ASf*4»  XfMfi  4  tjf  Si'h/8 

(c  )  b  *«>4  *  -X«  Sinfi 

bcsM  * -X,  Si + 

M'+to\U)\  U)' :  b‘« 

-{(f-Ut  i  al-4*  * 

i.t.  <o%2  jb  -  *&!&=* t-  tMfi* 

cl*(4)*  i  («*-!»*)  J»iit  s  -  -f 

(a'-b*)  «  *?-*$ «■ 

«s2w  • 

Thus  in  total: 


(d> 

Then: 


Next: 
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(40) 


to-n  Z  £ 


.  ited&*XLHt )x 


(41) 


ta*  **  Xi wZtjMft 


Co$-2 A  = 


d'-V 


Thu*  ^  Is  the  angle  from  the  positive  X-axis  to  the  major  axis  of  the 
ellipse  (positive  witho))  and  is  the  phase  angle  for  the  radius 
vector,  measured. positive  from  the  maj  .r  axis  in  the  direction  of  ut  . 
The  computer  output  gives  dtb t  $  and  (X  for  both  the  deflection  ard  the 
bending  moment. 

Coupling  Stations.  The  programs  allow  for  couplings  in  the  rotor.  At 
these  stations,  the  bending  moment  vanishes,  i.e.  Mnm0  (the  coupling 
point  is  taken  just  to  the  right  of  the  mass  station).  When  the  prog¬ 
ram  encounters  a  coupling  station,  say  station  l  ,  the  following 
equations  are  set  up: 


Mari,i  Mks1 - 

’*»■ 

Mtteif - AXaclp 

V 

AfjK*/  sMittifit 

Mwy  Mui(i - fAai.g 

A' 

St  — 

Mai,*  —  Mai,  t 

Xs. 

M* +  Miiijio 

- 

Mfei,S - Afp*,# 

</« 

*  MicL'ie 

AM,1  ***"  tit 

— 

fs 

• 

/ifas 

or  upon  solving 

(43)  &L  *  Ail  ( i, j  - f. 2, 5,4  ) 

where  9ct  etc.  and  Xt"Xu  tX\9Xn  etc.  Then  the  bending 

moment,  shear  force,  slope  and  deflection  before  the  coupling  station 
become  functions  of  y*i  and  Lf%i  only.  As  an  axaaiple,  let 

the  shear  force  at  a  station  be: 

Vita  mW9*t  6*1  +V*fa  +  VfXci  *  14/*  +  *■  /f  +\/l0 

Introducing  eq. (A3)  gives: 


V** = [v, + +  «»  V2  +  «„V,  +  a,,V4  ]yc, + [ V( + 4„V; + «„  ]  14, 

+  [V,*a„y+a„V,  +  «„V,  ^1%+fV,  +aHV;  -  4*^+^ 

+  [V,  -*V„+  k,Vi  +t,V,+  t,V,+ l,V4] 

and  similarly  for  V,c„  - ^.Instead  of  oW  $, 

we  have  aa  new  variables  Che  slopes  Just  Co  Che  right  of  Che  coupling 
station  m,  i.e.  ,  j«,.Thsn  the  celculetion  proceeds  as  be¬ 

fore  until  either  a  new  coupling  station  or  the  right  end  of  Che  rotor 


is  reached. 


Transmitted  Force  and  Pedestal  Motion 

Let  the  force  transmitted  to  the  pedestal  at  station  ft  be  denoted  F. 
From  Eq.  (12)  it  is  seen: 

Hire  *  V* c+-i  ~  Mrce  +  1%  <*?Xqi  U?  U** 

(45)  FVJ  *►  W'Xj„  Wl  Uij» 

Rj«  s  %*,*••! 

^  i  *  %*(•>-•  ~  ni" w 

Denoting  the  amplitude  of  the  pedestal  massa  Vp  and  (sea  Fig.  3) 
we  get: 

X  -  Hrc  -  i 
*•  "  Xy-MxOjS-iU^ 


u  x  le  «5 

“  X^MjCo’+uody 

v  -  Hu  (**  ~  ^W)  -Fw  u>t 


= 


x.i  = 


Hf<  * 


«jf*  * 


(X,-Mru)k)V(<odk)1 

_  F»c  <od»  j;  f>»  (Xn-Mscu*? 

<*c*— 

*'f  *t  (uid^)1 
-  fw  Cpd^  Hr  ^is(>tf 


The  force  transmitted  to  the  base  becoaMa: 


P*  *  X*Vp  i-tOidy  Xp 

P«>  *  **  ‘fr  +  iw  4 


26 


(47) 


Pxc  -  F,c  4 

Pyj  ’  Rr,  4  KU)\S 

P«,c  *  *V<  4  “V 

Pjs  8  R,i  4  wl(^M 

Energy  Balance 

Let  the  relative'  amplitude  between  rotor  and  pedeatal  oat a  be  X>= X~Xp 
and  tj,=  ij-^p  at  a  bearing  station.  Then  the  energy  dlsai* 
pated  in  the  bearing  and  the  pedeatal  per  revolution  becomea: 

Energy  Diaaipated  - 

(48>  7r{  0)Cw  (x'c+X'f)  ig)  +  ((o(ni+U}CF)(x‘ci1't+X‘s  if't) 

(x'c  -  x’  i/c)  4  (xfc+xjs )  +  ) } 

+7i{<oD„(e,cI+e(1)  +(4b,,Mv‘M<ts ) 

A  aunaation  over  all  bearinga  givea  the  total  diaaipated  energy. 

At  each  unbalance  atation  there  la  an  energy  input: 

(49)  Energy  Input:  7T(  Ui  (Xj-lfc)  +  Cu*  l/y  ( *£+  lj%)  } 

Suomlng  over  all  unbalance  atatlona  givea  the  total  energy  input 
which  muat  equal  the  diaaipated  energy. 

gfflgaEiLjmpyi 

The  input  data  la  prepared  according  to  the  following  inatructlona. 

Mote  that,  unleaa  specifically  stated,  no  input  card  may  be  omitted. 

Card  1  and  2:  (72  cola.  Hollarlth)  Identification:*  Any  descriptive 

text  may  be  punched  in  cola.  2-72.  These  two  cards  must  always  be 
Included. 

Card  3:  (1015)  Control  parameters  - 

Word  1.  Number  of  rotor  mass  stations  *  The  number  of  naas  stations  is 
determined  by  the  above  considerations.  Also,  there  must  be  a  mass 
station  at  each  rotor  end,  at  each  bearing,  at  each  unbalance  and  at 
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each  coupling  point.  The  mass  at  a  station  may  ba  zero.  The 
maximum  number  of  maaa  stations  Is  80. 

Word  2.  Number  of  bearings  -  This  lntagar  denotes  the  total  number  of 
bearings  along  the  rotor.  A  maximum  of  28  bearings  is  possible. 

Word  3.  Number  of  unbalance  stations  -  This  integer  gives  the  total 
number  of  mass  stations  at  which  unbalance  Is  applied.  A  maximum  of 
80  unbalance  stations  is  possible. 

Word  4.  Number  of  coupling  stations  *  This  Integer  gives  the  total 
number  of  coupling  points.  It  cannot  exceed  20. 

Word  5.  Pedestal  flexlble/rigld  -  If  this  integer  is  zero,  the  program 
assumes  the  pedestal  to  be  rigid  for  both  translator?  and  rotational 
motion  and  no  pedestal  data  is  Included.  If  the  Integer  is  1.  the 

pedestal  has  flexibility  and  damping  and  pedestal  data  must  ba  provided. 

\ 

Word  6.  Support  tilting  -  If  this  Integer  Is  zero,  neither  the  bear* 
Ings  nor  the  pedestals  resist  rotation.  In  that  case,  neither  the 
Input  for  the  bearing  dynamic  coefficients  for  rotational  motion  nor 
the  pedestal  data  for  rotational  nation  can  ba  included.  If  the 
integer  is  1,  the  bearings  and  the  pedestals  have  flexibility  and 
damping  for  rotational  motion. 

Word  7.  Gyroscopic  sonant  -  If  this  integer  la  zero,  no  gyroscopic 
moment  is  Included  In  the  calculation.  If  gyroscopic  axxaent  la 
desired,  the  integer  should  be  1. 

Word  8.  Number  of  computations  -  It  was  indicated  above  that  the  eight 
bearing  parameters  were  dynemic  coefficients  and  so  could  account  for 
the  variation  of  parameters  with  running  speed  in  an  approximate  manner. 
However,  If  e  more  precise  representation  of  these  parameters  is  des¬ 
ired,  these  values  can  be  entered  each  time  a  new  running  speed  is 
designated.  In  order  to  facilitate  this,  there  is  provision  in  the 
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program  for  entering  only  the  bearing  or  bearing  and  pedestal  data 
and  the  corresponding  running  speed  without  re-entering  the  rotor, 
coupling  or  unbalance  data.  Then  this  word  8  of  the  control  para¬ 
meters  designates  the  number  of  sets  of  parameters  which  are  to  be  run. 

If  this  value  is  1,  the  program  assumes  that  the  bearing  data  is  being 
entered  as  coefficients  of  quadratic  equations  in  Oi  .  Note  below  that 
the  input  format  of  the  bearing  data  differs  depending  on  whether  this 
value  is  equal  to  or  greater  than  one. 

Word  9.  Diagnostic  -  If  this  integer  is  zero,  no  diagnostic  will  be 

performed.  A  value  of  1  will  provide  the  diagnostic  output:  the 

diagnostic  increases  the  amount  of  output  a  considerable  amount  and 
is  provided  primarily  for  use  in  trouble-shooting  the  program  and  so 
this  value  should  alway  be  zero. 

Word  10.  Input  -  If  this  integer  is  zero,  the  program  will  return  to 

read  in  anew  set  of  input  upon  completion  of  the  computation.  For 

the  last  set  of,  input  this  value  should  be  1. 

Card  4.  (1P4E15!?) 

2 

Word  1  is  Young'  modulus  E  in  lbs/in  .  It  is  constant  throughout  the  rotor. 
Since  the  program  never  uses  E  by  itself  but  always  in  the  product  El 
(I*cross-sectlonal  moment  of  inertia)  any  actual  variation  in  E  can  be 
absorbed  by  changing  I  accordingly. 

Word  2  is  the  scale  factor  for  the  determinant  in  the  simultaneous  equation 
subroutine.  In  general  this  item  is  1.0.  It  is  a  factor  by  which  the 
determinant  is  multiplied  to  control  computer  over/underflow.  The  simul¬ 
taneous  aquation  subroutine  is  used  4  places  in  the  program:  oner  when 
solving  for  the  unknown  end  deflections  (i.e.Eq. (35))  and  3  timos  when  solv¬ 
ing  for  the  unkown  slopes  in  the  coupling  calculation  (i.e.  Eq.  (42)).  If 
an  over/underflow  occurs  during  the  calculation  the  program  output  will  con¬ 
tain:  "OVER/ UNDERFLOW  IN  XSIMEQF  AT _ (integer)"  where  the  value  of  the 

Integer  is  1  to  4.  If  it  is  1,  2  or  3  the  error  is  in  the  coupling  cal¬ 
culation.  If  it  is  4  the  error  is  in  solving  Eq.  (35).  Changing  the  scale - 
■  i 

factor  may  eliminate  the  trouble. 
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If  che  determinant  ia  aingular  the  output  givaa:  "MATRIX  IS  SINGULAR 

IN  XSIMEQF  AT _ (integer)".  If  either  of  the  two  errors  occur  the 

program  proceeda  with  a  new  rotor  apeed. 

B2I2&-M1A 

The  rotor  data  will  differ  depending  on  whether  the  effect  of  the 
gyroacopic  moment  ia  included  in  the  computation.  For  the  caae  where 
no  gyroacopic  moment  la  Included;  l.e.  wh»re  word  7  of  card  3  ia  zero, 
the  rotor  data  la  entered  aa  foil  owe: 

Card:  (1P3E14.6)  -  An  input  card  muat  be  given  for  each  aaaa  atatlon. 
Each  card  haa  3  itema. 

Word  1  -  che  aaaa  at  the  atatlon  in  lba. 

Word  2  *  the  length  of  the  ahaft  section  to  the  right  of  the  atatlon 
in  inchea. 

Word  3  ~  the  cross-sectional  moaMnt  of  Inertia  of  the  ahaft  section  to 

4 

the  right  of  the  station  in  in  . 

For  the  laat  maee  atatlon  the  ahaft  length  and  the  croes-sectlonal 
moment  of  inertia  has  no  meaning  and  may  be  set  equal  to  aero. 


If  gyroscopic  motion  is  Included  and  word  7,  card  3,  la  not  equal  to 
zero,  then  each  card  contains  two  more  Items  In  addition  to  the  3 
items  indicated  just  above.  Also  for  this  case,  the  rotor  data  cards 
are  immediately  preceded  by  a  card  which  contains  two  values  defined 
as  follows: 

Card:  (1S,1PS23.6) .  Gyroscopic  moment  parameters  - 

Word  I  *  Number  of  iterations  -  For  each  rotor  speed  the  program  first 

calculates  the  unbalance  response  without  gyroscopic  moment,  lased  on 
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the  thuj  obtained  rotor  slopes,  the  gyroscopic  moment  Is  computed  and 
applied  to  the  rotor,  resulting  in  new  values  of  the  slops  and  the  process 
is  repeated.  The  program  counts  the  number  of  iterations,  excluding  the 
ca'culatlon  without  gyroscopic  moment.  If  the  count  exceeds  this  input 
item,  the  results  obtained  are  printed  out,  the  iteration  count  is  reset 
to  1  and  a  new  rotor  speed  calculation  starts. 

Word  2  -  Convergence  limit  -  After  each  gyroscopic  moment  Iteration,  the 
following  relative  error  is  calculated: 

iaP£*m*r!e£*V . 

where  6ci,  $31,  and  (Pf,  are  the  slopes  and  X*i ,  Xt> ,  4  a  and  ^31  are 

the  deflections  at  the  left  rotor  end.  The  superscript  is  the  Iteration 
number.  For  each  iteration  the  computer  output  gives  the  iteration  number 
and  the  error.  When  the  error  is  less  than  or  equal  to  the  input  conver¬ 
gence  limit,  the  program  prints  the  results,  resets  the  iteration  count  to 
1  and  proceeds  with  a  new  rotor  speed. 

Following  this  card  are  the  rotor  data  cards. 

Card:  (1P5E14.6).  An  input  card  is  required  for  each  mass  station.  Each 

card  contains  5  items;  the  first  3  words  are  the  same  aa  those  for  the  non- 

gyroscoplc  moment  case  above  and  the  remaining  two  are: 

7 

Word  t  -  the  polar  mass  moment  of  inertia  in  lbs. in 

2 

Word  5  -  the  transverse  mass  moment  of  inertia  in  lbs. in 
LOCATION  OF  BEARING  SUPPORTS 

Card:  (1415).  This  list  provides  the  numbers  of  the  mass  stations  at 
which  there  is  a  bearing. 

UNBALANCE  DATA 

Card:  (15,  1P2E15.7).  A  card  is  provided  for  each  of  the  unbalance 
stations.  Each  card  contains  3  values: 

Word  1  -  an  integer  which  denotes  the  number  of  the  mass  station  at  which 
the  unbalance  applies. 

Word  2  -  the  cosine  component  of  the  unbalance  in  os.  in. 

Word_3  -  the  sine  component  of  the  unbalance  in  os.  in. 
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Bv  providing  two  unbalance  component!,  it  la  poaaible  to  take  lato 
consideration  the  circumferential  variation  of  unbalance  along  Che 


rotor. 


COUPLING  DATA 

If  the  rotor  does  not  contain  couplings,  (Word  4,  card  3  la  aero),  than 
no  coupling  data  is  necessary.  If  word  4,  card  3  is  not  sero,  the 
following  card  must  be  included. 

Card:  (1415)  •  This  is  a  list  of  integers  denoting  Che  suss  stations 
at  which  there  is  a  coupling. 

PEDESTAL  DATA 

If  the  pedestal  is  considered  to  be  infinitely  rigid,  then  no  pedestal 
data  is  required.  In  this  case  word  5,  card  3,  pedestal  flexible/ rigid 
is  sero.  Otherwise,  pedestal  data  is  required.  The  pedeetal  data, 
like  the  bearing  data,  Is  separated  into  translator?  and  rotational 
parameters.  Also,  as  before,  the  control  parameter  is  the  word  6,  —  - 
card  3,  support  tilting. 

Card:  (1P6E12.4)  -  A  card  must  ba  provided  for  each  bearing,  (ha  It 
are  6  values  as  follows: 

Word  1  -  the  weight  of  the  pedestal  in  the  X  coordinate  la  lbs. 

Word  2  -  the  pedestal  stiffness  along  the  X  coordinate  in  lbS/ln. 

Word  3  -  the  pedestal  damping  along  the  X  coordinate  in  lbs>sec/in. 

Word  4  -  same  as  word  1  but  for  the  If  coordinate. 

Word  5  -  same  as  word  2  but  for  the  14  coordinate. 

Word  6  -  same  as  word  3  but  for  the  14  coordinate. 


If  word  4,  card  3,  support  tilting,  is  not  sero,  then  all  of  the  cards 
concerned  with  the  tranalatory  parameters  are  followed  bjr  the  cards  for 
the  rotational  parameters.  Again  there  are  4  values  on  each  card  as 
follows: 

Word  1  -  the  mass  moment  of  ineTtla  of  the  pedeetal  sums,  aaaociated 

2 

with  the  X  coordinate  in  lbe.ln 
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Word  2  -  Che  pedeatal  aprlng  coefficient  for  rotational  notion,  aaaoc- 
iated  with  the  X  coordinate  in  Iba-in/rad. 

Word  3  -  the  pedeatal  damping  coefficient  for  rotational  notion,  aaaoc- 
latcd  with  the  X  coordinate  in  Iba-ln. -aec/rad. 

Word  4  -  aarae  aa  word  1  but  aaaociated  with  the  coordinate. 

Word  5  •  aame  aa  word  2  but  aaaociated  with  the  ^  coordinate. 

Word  6  -  aame  aa  word  3  but  aaaociated  with  the  coordinate. 

SPZED  AND  BEARING  DATA 

Each  bearing  la  represented  by  16  dynamic  coefficlenta,  8  for  trana- 
latory  motion  and  8  for  rotational  motion.  Of  the  8  coefficlenta,  4 
are  aprlng  coefficlenta  and  4  are  damping  coefficlenta.  Since  the  co¬ 
efficlenta  In  general  change  with  apeed,  each  coefficient  la  expreaaed 
by  three  componenta;  e.g. 

^  Kx*,i  W  +  Km,a  U)* 

where  4)  la  the  rotor  apeed  In  rad/aec,  Kw,«  in  lba/ln. ,  K«,i  In 

2  2 

Iba-aec/in. rad.  and  in  lba-aec  /in. rad  .  Similar  equatlona 

hold  for  the  other  IS  coefficlenta.  Aa  indicated  earlier,  if  it  la 
dealrad  to  enter  the  bearing  data  at  aach  value  of  frequency,  there  la 
provlalon  for  thla  in  the  program. 

If  word  8,  card  3  la  1,  the  program  aaaumea  the  bearing  data  ia  prov¬ 
ided  aa  frequency  dependent  coefficlenta.  In  thia  caae,  a  card  la 
provided  with  the  apeed  range  and  increment  and  thla  la  followed  by 
the  beering  data.  An  Input  card  ia  given  for  each  coefficient  at  aach 
bearing.  Each  card  contalna  three  itama,  namely  the  above  mentioned 
three  apeed  componenta.  The  aaquence  of  the  Input  carda  ia:  flrat 
all  the  carda  for  the  tranalatory  motion  and  then  all  the  carda  for 
the  rotational  motion.  The  carda  for  tha  rotational  motion  are  not  , 
required  if  word  6,  card  3,  aupport  tilting,  ia  aero.  The  carda  ahould 
be  given  in  the  following  order:  £*  ,(wfV|  ,0^ 

for  bearing  1,  Kn  j4*Crv - for  bearing  2,  etc.  to  the  laat  bear¬ 

ing,  then  (If  word  6,  card  3  #  0),  Mm, 
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tor  bearing  1,  etc.  to  and  including  the  laat  bearing. 


Card:  (1P3E14.6)  -  Speed  date. 
Word  1  -  initial  speed. 

Word  2  -  final  speed. 

Word  3  -  speed  increment. 


Cnrd:  (1P3FU.6)  -  Bearing  data  -  in  the  order  defined  abcaa  with 
three  values  on  each  card  as  follows: 

Word  1  -  the  icient  A0  of  the  expression  +A,  u)  +  ui* 

Word  2  -  the  coefficient  A(  of  this  expression. 

Word  3  -  the  coefficient  At  of  this  expression. 


If  word  8,  card  3,  is  greater  than  1,  the  program  ass  jams  the  bearing 
data  will  be  provided  for  each  value  of  spaed.  In  this  case,  a  card 
is  provided  with  a  single  spaed  value  and  this  is  followed  by  the 
bearing  data  as  follows:  all  of  the  transl story  stiffness  end  damping 
coefficients  are  provided  in  the  order;  two  cards  for  each  bearing. 

The  first  card  contains  ths  X  coordinate  translatory  coefficient* 

Kxx  jUiCrr,  Ketf  and  (uC^* nd  ths  second  card  the  if  coordinate  translatory 
coefficients  /oCyj;  ky*  and  toCyg,  both  cards  for  bearing  one  followed 
by  two  cards  for  bearing  two,  etc.  ,  to  the  total  mabar  of  boarlnge  . 
Again,  if  word  6,  card  3  le  zaro,  no  rotational  parsneters  are  required, 
otherwiae,  they  are  provided  in  a  siailar  meaner:  one  card  of  velnee 
and  a  aacond  card  of  for  bearing 

1  followed  by  two  card*  each  for  the  resMlnltg  bearings. 


The  card  format  in  this  caae  la  than: 


Card:  (1PE14.6)  -  Speed. 

Card:  (1P4E14.6)  -  Bearing  data  in  the  order  defined 
values  on  eech  card  aa  follows: 


Card  1  -  Word  1  - 

Word  2  -t-Cw,* 
Word  3  -  ^ 

**>«»  4  -<oC,,te 


The  computer  output  is  largely  self-explanatory  and  each  output  item  la  identified 
by  a  descriptive  text.  Two  sample caae a  are  shown  in  Appendix  A.  ,The  output 

first  lists  all  the  inputdata,  i.a.  the  two  heading  cards,  the  control  words, 

I 

Youngs  Modulus,  the  rotor  data,  the  bearing  stations,  the  unbalance  data,  the 
coupling  stations,  the  pedestal  data,  tha  speed  data  and  finally-  the  bearing 
data.  Thereafter  follow  the  results  of  the  calculations  with  on*  aat  of  output 
for  each  specified  rotor  speed.  First,  the  spaed  is  given  in  RFM  which  say  be 
followed  by  the  input  bearing  data  if  it  is  new  fpr  every  spaed.  Next,  a 
9  column  list  gives  tha  rotor  amplitude  and  bending  moment  at  each  rotor  station. 
Both  the  amplitude  and  the  bending  moment  require  four  quantities  for  their 
description.  Since  each  rotor  station  whirls  in  an  elliptical  orbit  it  is  con¬ 
venient  to  express  the  four  quantities  in  terns  of  tha  dimensions  of  tha 
ellipse.  Then  the  four  quantities  become: 

'I 

1.  the  major  axis  of  tha  ellipse:  a(i.e.  the  maximum  aapllmuda  or  the 

maximum  bending  moment  during  one  revolution.  ’ 

2.  the  minor  axis  of  the  ellipse:  b 

3.  the  angle  between  the  x-axls  of  the  overall  reference  system  and  the 
major  axis  of  the  ellipse:  0,  degrees  (in  output  identified  by: 

ANGLE  X -MAJOR) 

4.  the  phase  angle  with  respect  to  tha  cosine -component  of  the  unbalance: 
cA  ,  degrees. 

The  amplitude  is  given  in  thousands  of  an  inch  (mils)  and  tha  bending  moment 
is  given  in  lbs. in. 

The  selected  method  of  presentation  is  illustrated  by  Fig.  5  end  is  given  in  de¬ 
tail  in  the  analysis  by  Eqs.  (37)  to  (41).  However,  e  general  description  will 
also  be  given  here. 

The  presentation  is  based  on  two  referanca  coordinate  systems.  The  first  refarenca 
system  is  the  stationery  x-y  system  fixed  with  respect  to  grand,  end  which  has  at 

I 

each  rotor  at'  on  its  origin  in  the  center  of  the  statically  deflected  rotor 
(l.e.  the  deflection  due  tu  gravity).  The  x-y-syetem  is  "communicated"  to  the 
rotor  via  tha  specified  values  of  the  bearing  spring  end  damping  coefficients 

36 

V 


(  K«, t  OjCmfitc, )  and  Che  correapondlng  padeatal  data.  In  other  words,  the 
directions  of  the  x~axls  and  the  y-axls  are  chosen  when  preparing  the  computer  In¬ 
put  and  the  choice  reflects  in  the  input  values  used  for  rtc.  Then  the 

elliptical  rotor  orbit  ia  centered  in  the  origin  of  the  x-y-syrtem  (i.e.  the 
steady  state  shaft  center),  it  has  a  major  axis  a,  a  minor  axis  b,  and  the 
orientation  of  the  ellipse  is  defined  by  the  angle1  3  between  the  x-axia  and 
the  major  axis,  measured  in  direction  of  rotor  rotation.  Note,  that  both  a,  b 
and  0  vary  along  the  rotor.  A  negative  value  for  b  signifies  backward  whirl. 

Thus  a 4b  and  3  specify  the  dimensions  and  the  orientation  of  the  elliptical 
orbit  but  one  more  quantity  ia  needed  to  specify  the  position  of  the  moving  shaft 
center  on  the  ellipse  at  any  given  time.  The  phase  angle  s(  serves  this  purpose. 
Let  the  major  axis  be  the  x^-axis  and  the  minor  axis  the  y^-axis  (see  Fig. 5), 
i.e.  the  x^-y^-systea  is  obtained  by  rotating  the  x-y-systera  an  angle  3  in 
the  direction  of  rotor  rotation.  Then  the  instantaneous  position  of  the  shaft 
center  Is  given  by: 

tj,*  fcSrtfait  +st) 

Note  that  the  orientation  of  the  x^-y^-sysjtem  changes  along  the  ro*ror  since  3 
does  with  respect  to  the  x-y-syatea  the  instantaneous,  shaft  center  position  is 
given  by: 

X  *  \/ UctsfiP+lbsiHltf  Csj  (tut*  ♦  WYl 

If  *  V  (•'»  (tut +4  +  W'ff  ia*  /»)) 

In  addition  to  determining  the  instantaneous  position  of  the  shaft  center  with 
respect  to  a  stationary  coordinate  system  Jt  also  may  be  desired  to  know  the 
position  with  respect  to  the  rotor  unbalance.  The  location  of  the  unbalance  in 
the  rotor  is  defined  by  a  coordinate  system  which  ia  fixed  in  the  rotating  shaft 
and  whose  axes  are  called  "the  cosine  axis"  and  "the  sine  axis".  Hence,  esch 
unbalance  consists  of  two  components:  a  coalne  component  ends  sine  component 
(in  the  analysis  the  symbols  Ux  and  are  used,  respectively ,  see  Eq.s(10) 
and  (11)).  The  instantaneous  orientation  of  the  cosine-system  is  defined  by  the 
angle  (  lot  )  between  the  fixed  axis  and  the  coalna-axis.  Thus,  the  instsntsneous 
phase  angle  between  the  amplitude  vector  (i.e.  the  radius  vector  from  the  center 
of  the  elliptical  orbit  going  through  the  instantaneous  shaft  center  position) 
and  the  total  rotor  unbalance  vector  is: 
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jrtglt  by  which  amplitude  vector  lags  unbalance  vector  ■ 


Here,  £  UtH  and  ilUj*  indicates  the  summations  of  the  cosine-components 
and  the  s ine -components ,  respectively,  of  all  unbalances.  It  is  seen  that  the 
lag-angle  is  not  constant  as  the  shaft  center  moves  around  its  orbit.  It 
attains  its  maximum  and  minimum  values  when: 

Although  the  discussion  above  is  prlatarlly  aimed  at  describing  the  motion  of 
the  shaft  center  (i.e.  the  computer  output  for  the  amplitude)  the  same  des¬ 
cription  applies  to  the  output  for  the  bending  moment.  However,  for  each  rotor 
station  the  output  lists  one  line  for  the  amplitude  but  two  lines  for  the  bend¬ 
ing  moment.  Whereas  the  output  for  the  amplitude  applies  at  the  rotor  station 
itself  the  bending  moewnt  has  one  value  isBedlately  to  the  left  of  the  station 
and  another  value  immediately  to  the  right  of  the  station.  The  output  gives 
the  left  hand  value  first  (i.e.  the  output  gives  H*  and  M*  respectively,  se» 
Fig.  2).  The  two  values  are  in  general  the  same  unless  the  particular  station  is 
a  bearing  station  which  resists  tilting.  The  last  listed  value  of  the  bending 
moment  should  always  be  zero  (i.e.  corresponding  values  of  the  major  and  minor 
axis  should  be  zero).  In  general  they  are  not  exactly  taro  but  vary  small.  The 
amount  by  which  the  values  differ  from  zero  gives  an  indication  of  the  accuracy 
of  the  calculation.  Note,  that  for  this  reason  the  last  values  of  the  angles 
and  el  are  meaningleo* . 
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Following  Che  output  for  the  amplitude  and  for  the  bending  moment  come  tl 
suits  for  the  force  transmitted  to  the  bearing  housing  (equal  to  the  dynai 
bearing  reaction).  If  the  pedestals  are  flexible,  the  force  transmitted  I 
the  foundation  and  the  amplitude  of  the  pedestal  mass  are  also  given.  Eac 
the  three  quantities  are  presented  in  two  ways:  first  in  terms  of  the  coi 
ponding  ellipse  (i.e.  in  analogy  to  the  rotor  amplitude)  and  secondly  by  l 
x  and  y-components .  Thus,  if  the  transmitted  force  is  F  the  output  gives 
quantities;  Che  major  axis,  the  minor  axis,  the  orientation  angle  the 
angle  of  j  If^l  j  If^l  and  <(y  where  the  last  four  items  are  defi 

force  in  x-direction:  Fx  “  IF»I  Cdsfwt+d*) 
force  ir<  y-direction:  ^  I  Sin6ut**sfj) 

The  transmitted  force  is  given  in  lbs  and  the  pedestal  amplitude  in  thousar 
of  an  inch  (ails). 

The  next  line  of  output  serves  as  a  check  on  the  calculation.  It  gives  the 
energy  per  revolution  put  into  .the  system  by  the  unbalance  forces  and  the  e 

v 

dissipated  per  revolution  in  the  bearings  and  pedestals.  Theoretically,  th 
two  values  should  be  equal  but  numerical  inaccuracies  in  the  calculations  a 
cause  discrepancy.  Normally  they  'differ  in  the  fifth  or  sixth  decimal  placi 
The  energy  is  given  in  lbs . Inch/revolution. 

To  convert  it  into  HP  multiply  the  energy  value  by  the  speed  in  RPM  and  divi 
by  3.96  •  105. 

If  the  input  does  not  include  any  gyroscopic  moment  effects  the  calculations 
repeated  for  a  new  rptor  speed  and  the  output  will  follow  the  description  gi 
above.  If  the  gyroscopic  moment  is  included  there  are  two  sets  of  output  fo: 
each  rotor  speed,  each  set  having  the  format  as  explained  above.  The  first 
applies  to  a  rotor  without  any  gyroscopic  moments,  and  the  second  set  givei 
the  final  result  for  the  calculation  with  the  gyroscopic  moment  included.  T. 


two  sets  are  separated  by  a  two  column  llat  giving  the  sequential  results  of  the 
iterations  needed  to  perform  the  gyroscopic  moment  calculation.  The  first 
column  idencifiea  the  iteration  number  and  the  second  column  gives  the  relative 
cor'. ergence  of  the  iteration  procedure  as  explained  in  describing  the  computer 
input . 
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COMPUTER  PROGRAM:  THE  STABILITY  OF  A  ROTOR  IN  FLUID  FILM  BEARINGS 


This  section  sets  forth  the  basic  analysis  and  the  detailed  instructions  for 
using  the  computer  program:  PN0017:  "The  Stability  of  a  Rctor  in  Fluid  Film 
Bearings"  for  the  IBM  704  digital  computer.  The  program  calculates  the  speed 
at  onset  of  instability  (the  threshold  speed)  and  the  corresponding  whirl 
frequency. 

Each  rotor  support  consists  of  a  fluid  film  bearing  mounted  in  a  pedestal,  both 
members  possessing  flexibility  and  damping.  The  bearing  fluid  film  is  represent! 
by  8  dynamic  coefficients  and  it  is  the  value  of  these  coefficients  which  prim¬ 
arily  govern  the  instability  mechanism.  For  a  given  application  they  vary  with 
the  speed  of  the  rotor  and  they  be  specified  in  the  computer  input  for  each 

speed  to  be  tested.  The  bearing  pedestals  are  represented  by  2  spring  coefficiei 
and  2  damping  coefficients  (in  the  vertical  and  the  horizontal  direction,  re¬ 
spectively)  and  the  corresponding  masses  may  also  be  given. 

The  program  is  to  a  large  extent  compatible  with  the  unbalance  response  program 
such  that  much  of  the  input  data  used  in  the  latter  program  also  applies  to  the 
stability  program. 

THEORETICAL  ANALYSIS 

The  analysis  is  an  extension  of  the  methods  used  in  the  previous  section  to  deter- 

mine  the  unbalance  response  of  the  rotor.  Thus,  the  following  discussion  assumes 

familiarity  with  the  earlier  given  analysis.  The  rotor  is  again  represented 

by  a  finite  number  of  mass  stations  connected  by  weightless  but  stiff  shaft 

sections  which  can  be  brought  to  approximate  the  actual  rotor  to  any  degree  of 

accuracy  depending  on  the  number  of  mass  stations.  Each  bearing  is  represented 

by  8  dynamic  coefficients:  K  ,  C  ,K  ,  C  ,  K  ,  C  ,  K  and  C  ,  which 
1  1  xx  xx  xy  xy  yx  yx  yy  yy 

depend  on  the  operating  speed  of  the  rotor:  o>,  radians/sec. 

< 

The  purpose  of  the  analysis  is  to  establish  the  onset  of  instability  of  the  rotor- 
bearing  system.  No  external  forces  act  on  the  rotor  (i.e.  there  are  no  unbalance 
forces),  Instead  the  dynamical  equilibrium  of  the  steady-state  operation  of  the 
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system  is  tested  for  a  series  of  discrete  speed  values  over  a  speed  range.  This 
is  done  by  disturbing  the  rotor  with  an  assigned  frequency  V  radians/sec.  To 
this  end  the  previously  established  equations  are  used  as  summarized  in  Eqs.(32) 
by  replacing  U)  with  V  or  since  CO  is  given,  the  disturbing  frequency  is  speci¬ 
fied  as  a  ratio  of  the  speed: 

(note:  V  =(uj)cd  ) 

The  rotor  unbalance  components  and  U  are  eliminated.  Applying  Eqs.(32) 
and  the  boundary  conditions  Eqs .  (33)  and  (34)  yields  Eq.(35)  where  the  right 
hand  side  is  now  equal  to  zero. 

Since  the  assigned  frequency  is  a  pure  frequency  and  does  not  contain  a  transient 
term  the  outlined  procedure  applies  to  the  threshold  of  instability.  In  other 
words,  the  calculation  determines  the  state  of  neutral  stability  at  which  the 
effect  of  any  disturbance  continues  indefinitely  without  either  increasing  or 
decreasing.  Hence,  at  the  point  of  neutral  stability  there  must  be  a  finite, 
although  undetermined,  solution  for  the  rotor  amplitudes,  i.e.  the  8  end  values, 

0  .,  8  ,  ,  . ,  y  ,  cannot  all  vanish.  This  implies  that  the  determinant  on 

the  left  hand  side  of  Eq.(35)  must  be  zero  for  the  system  to  be  neutrally  stable. 

On  this  basis,  a  calculation  procedure  can  be  dev  sed.  Select  a  sufficiently 
low  value  of  the  rotor  speed  that  the  system  is  known  to  be  stable  and  scan  the 
entire  frequency  range  to  obtain  the  value  of  the  determinant  for  each  frequency. 
Repeat  the  calculations  for  an  increased  rotor  speed  and  proceed  in  this  way 
until  a  speed  is  encountered  at  which  the  determinant  becomes  zero.  At  that 
particular  speed  the  rotor-bearing  system  is  on  the  threshold  of  instability 
and  any  further  increase  in  the  rotor  speed  will  make  the  system  unstable. 

Although  the  method  is  quite  simple  in  principle  certain  difficulties  arise  in 
applying  the  method.  Considering  the  matrix  in  Eq.(35)  it  is  seen  to  be  an 
8  b  8  matrix  in  which  all  elements  are  real.  Actually,  it  can  be  written  as  a 
4  by  4  matrix  with  complex  elements.  Therefore,  its  determinant  will  always 
be  positive  and  in  applying  the  outlined  calculation  procedure  the  aero-point 
of  the  determinant  will  appear  js  a  minimum.  There  will  be  no  cross-over 
from  a  positive  value  to  a  negative  value,  or,  in  other  words,  at  the  onset  of 
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instability  the  chosen  form  of  the  determinant  has  a  repeated  root.  Thus,  it 
is  necessary  to  plot  the  determinant  as  a  function  of  the  rotor  speed  to  be  able 
to  establish  when  it  becomes  zero.  However,  the  determinant  also  depends  on 
the  disturbing  frequency  and  it  only  vanishes  for  a  particular  value  of  that 
frequency.  It  is,  therefore,  also  necessary  to  know  the  frequency-value  at 
which  the  determinant  should  be  evaluated  at  each  rotor  speed. 


For  this  reason  it  is  chosen  to  calculate  the  real  and  imaginary  part  of  the 
4  by  4  complex  determinant  which  is  equivalent  to  the  8  by  8  real  determinant. 
Denote  the  real  matrix: 


11  ••• 

**  bV2n 

2n, 1* " ’ 

7  b2n,2n 

(n  ■  4  in  the  present  case) 


and  let  the  corresponding  complex  matrix  by: 


(n  ■  4  in  the  present  case) 


With  the  convention  followed  in  the  present  analysis  the  two  matrices  are  re¬ 
lated  by: 

bll  "  Re{alJ  bu"Im'{anJ 
b21  “  *Inlall}  b22  "Re  Kli 


Let  these  matrices  be  the  coefficient  matrices  in  a  set  of  linear  equations: 

(50)  B2n  •  X  -  U 

(51)  A  •  Z  -  W 

n 

where: 
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such  that: 


-  x  -  iy  j  -  1  ,n 


J  ') 


"j  '  “)  '  lvj 


J  -  1,0 


Since  the  two  equations  are  equivalent  they  must  yield  the  same  solution.  There¬ 
fore,  first  set:  w^  ■  m  ~ ~ 


w  , ■  0  and  w  ■  i  which  means:  u, ,  u ,--u  ■  0, 

n-1  n  1  2  n 


,  “  v  -  -  -v  ,  *  0  and  v  "  *1.  Solve  Eq.(51)  for  z  : 

1  2  n-1  _  n  _  ns/  n 

x  i ^***<  -  i  Ah.,  ( 


(52) 


V 


.1  , 

*  A,/ 


where  A  is  the  determinant  of  A  ,  A  .is  the  determinant  of  A  with  the  n'th 
n  n  n-1  n 

column  and  row  removed,  and: 

A^+tA,,.*  Xh 

i  .  e .  ,  A  and  *ni  are  the  real  and  imaginary  parts  .respectively  of  the  complex 
matrix  and  it  is  desired  to  calculate  them. 


Next,  solve  Eq.(50)  for  x  and  y  with  v  ■  -1: 

n  n  n 


(53) 


(54) 


X< 


.-S* 

B,„ 


*  Bjn 


where  B.  is  the  determinant  of  B_  ,  B„  ,  is  the  determinant  of  B.  with  the  n‘ 
2n  _ ,  2n  2n-l  2n 

column  and  row  removed,  and  B.  .  is  the  determinant  of  B,  with  the  (n-1)  'th 
column  and  the  n'th  row  removed.  Furthermore,  it  is  known  that: 

iA.r=  c + it,,- 

Equate  Eqs.  (52),  (53)  and  (54)  to  get: 

i  i  Aw  i )  ^  B|»^ 

&  ”  aW*. 


th 


or: 

(55) 


B,.., 
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where  B,  -  is  the  determinant  of  B_  with  the  two  last  columns  and  rows  deleted 
Zn-Z  2n 


Thus,  Eq.(55)  gives  a  recurrence  formula  where  the  order  of  the  original  deter¬ 
minant  is  reduced  by  1.  If  it  is  applied  repeatedly  the  final  result  becomes: 

(56)  Ahr  *  I]  = - 

k*l  &2K-1 


with  the  definition: 


B.-I 


(57)  is  the  determinant  corresponding  to  the  first  (2k-l)  columns  and  rows 

of  the  matrix  B_ 

2n 


(58)  B^.^  is  the  determinant  corresponding  to  the  first  (2k-l)  columns  and  rows 
of  the  matrix  B2n  but  where  the  (2k-l) ' th  column  has  been  interchanged  with 
the  2k' th  column. 


(59)  B2k2  is  the  determinant  corresponding  to  the  first  (2k-2)  columns  and  rows 
of  the  matrix  B^  .  ■ 

Thus  a  method  has  been  established  to  calculate  the  real  and  imaginary  part  of  an 
n  by  n  complex  determinant.  The  method  is  used  to  evaluate  the  determinant  of 
the  8  by  8  matrix  in  Eq.  (35). 

It  may  be  noted  that  Eq.  (35)  and  Eq.  (50) are  identical  except  for  the  change  in 

nomenclature.  Hence,  in  Eq.  (50)  X  represents  the  8  rotor  end  coordinates: 

0cl,  8g ^ ,  . .  y-  ygl  and  V  represents  the  moment  and  shear  components  at  the 

otheT~end  of  the  rotor:  M  ,M  . V  ,V _ .  Eq.  (50)  is  solved  in 

xcr'  xsr  ycr*  ysr 

Eqs.  (53)  and  (54)  for  the  case  of  v^*  -1,  i.e.  for  Vy8r"  ’1  which  means  that  a 
force: 

fore*  *  -  sinvt 

has  been  applied  to  the  rotor  end.  The  corresponding  amplitude  at  the  same 

rotor  end  can  be  computed  as  described  in  the  rotor  response  analysis.  Let  the 

y-component  be:  amplitude  ■  y  cos  Vt  +  y  sin  V  t 

cr  s  r 
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The  energy  input  imparted  to  the  rotor  motion  is: 

(60)  E cycle  “I  (~Vl£r$inyt  tVijir(oiA)  4t  : 

'# 

A  positive  energy  input  implies  that  energy  is  required  to  sustain  the  particular 
vibratory  motion,  i.e.  the  motion  is  stable.  On  the  other  hand,  a  negative  energy 
input  signifies  an  unstable  motion.  When  the  energy  input  is  zero  the  motion  is 
neutrally  stable.  It  must  be  noted  that  the  motion  itself  may  not  be  possible 
unless  the  previously  discussed  determinant  is  also  zero,  i.e.  the  outlined  energy 
criterion  can  not  be  used  alone  to  test  the  stability  of  the  rotor-bearing  system. 
However,  the  criterion  can  be  used  to  determine  the  value  of  the  Instability 
frequency. 


A  simple  example  may  serve  as  an  illustration.  Let  a  single  mass  M  be  supported 
on  a  spring  with  a  coefficient  K  and  a  dashpot  coefficient  C.  The  amplitude  is: 


ljc(cSvt  +ljtSihvt 


Apply  a  force  V  sin  Vt  such  that  the  equations  of  motion  become: 
ys 


(K-Mv')  vC 
-vC  (K-Mvx) 


r* 

m 

foj 

tvj 

in  analogy  to  Eq.  (35).  The  determinant  becomes: 

(K'Mv')  +  (vO  (note:  the  determinant  is  always  positive) 


which  vanishes  when  VC  ■  0  and  V  -  K/M,  i.e.  the  frequency  must  be  such  that 
it  simultaneously  makes  the  damping  VC  zero  and  also  equals  the  natural  frequency 
of  the  system.  Next,  solving  for  and  computing  the  energy  input  yields: 

1  _ vC 

Energy  input  per  cycle  ■  -■nv^r  fiv,,  (k-ivMvO* 


Applying  the  energy  criterion  from  above  the  motion  is  unstable  when  VC  is  negative 

and  vice  versa,  which  is  of  course  evident  in  this  simple  case.  However,  the 

system  is  only  neutrally  stable  if  in  addition  the  frequency  also  equals  the 

ir 

natural  frequency:  V  ■  rr. 
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In  calculating  the  energy  input  from  Eq.(60)  it  is  seen  that  a  difficulty  arises 
when  the  system  determinant  is  equal  to  zero  in  which  case  the  amplitude 
cannot  be  determined.  To  circumvent  this  problem,  y is  not  computed  as  its 

actual  value.  In  evaluating  0  .  ,  0  , . y  ,  from  Eq.  (35)  with  V  *  -l,  Cramer's 

rule  is  used  but  such  that  the  system  determinant  is  always  set  equal  to  l  re¬ 
gardless  of  its  true  value.  Thus,  the  resulting  calculated  energy  should  actually 
be  divided  by  the  system  deteiminant  to  obtain  the  real  value  of  the  energy  input. 

To  solve  for  the  onset  of  instability  the  procedure  is: 

a.  Select  a  rotor  speed  sufficiently  small  that  the  system  is  known  to 
be  stable. 

b.  Scan  the  frequency  range  and  determine  those  frequency  values  at  which 
the  real  part  and  the  imaginary  part  of  the  system  determinant  equal  zero. 

c.  Repeat  the  calculation  for  several  values  of  the  rotor  speed  covering 
a  sufficiently  large  speed  range. 

d.  Plot  curves  of  frequency  versus  rotor  speed,  obtaining  one  (or  more) 
curve  corresponding  to  the  real  part  of  the  determinant  being  zero  and  one 
(or  more)  curve  corresponding  to  the  imaginary  part  being  zero.  The 
intersection  of  the  curves  determines  the  threshold  speed. 

To  assist  in  searching  for  the  instability  frequency  the  single  bearing  frequency 
value  is  determined  for  each  bearing.  This' value  is  derived  by  considering  a 
stiff,  symmetric  rotor  supported  in  similar  bearings.  Let  the  rotor  mass  per 
bearing  be  M  whereby  the  equations  of  motion  becomes: 

ta-Mv'-K  vCj  (JC+ivO  f  V 

(61)  {  -0 

,  OCjk  +  ivCjJ  ( K^-Mvl+  i vCfl)  { •j 

The  real  and  the  imaginary  parts  of  the  determinant  have  to  equal  zero  separately: 

(62)  lV„  -Mv*)(lrw-Mv*)  -  If,,  If,,  -  v‘fa,C„  -C.jC^-O 

(63)  (K.„  -  My1)  Cy,  +  -K^C,.  %  C,,  -0 
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Solve  the  equations  to  get: 


(64) 

(65) 


Mv* 


CoCyy  ~  CO  C+iCyt 

Co  C»»  +  CjC ^ 


Computing  Eq.  (65)  first  and  substituting  into  Eq.  (64)  yields  the  instability 

V 

frequency  ratio  —  for  the  bearing.  Since  for  the  actual  rotor  the  bearing 
reactions  are  in  general  unequal  and  the  8  bearing  coefficients,  therefore, 
differ  in  value  among  the  bearings  the  instability  frequency  will  not  have  the 
same  value  for  all  the  bearings.  However,  the  instability  frequency  for  the 
rotor-bearing  system  will  lie  between  the  minimum  and  maximum  value  of  the 
bearing  frequencies. 


COMPUTER  INPUT 


The  input  data  is  prepared  according  to  the  instructions  given  in  the  following. 
Note,  that  unless  stated  otherwise  no  input  card  may  be  omitted. 

Card  1  FORMAT  (49  cols.  Hollerith)  Any  descriptive  text,  used  to  identify 
the  particular  calculation,  may  be  punched  in  cols.  2-49. 

Ca rd  2  FORMAT  (615)  Control  parameters: 

Word  1  (NS)  Number  of  rotor  mass  stations.  The  number  of  stations  is  selected 
according  to  the  previous  discussion.  There  must  be  a  mass  station  at  each  rotor 
end  and  at  each  bearing.  The  mass  at  a  station  may  be  zero.  The  maximum  nu...l>er 
of  stations  is  30. 

Word  2  (NB)  Number  of  bearings.  This  integer  denotes  the  total  number  of  bearings 
along  the  rotor.  A  maximum  of  10  bearings  is  allowed. 

Wora,  3  (NFR)  Number  of  frequency  ratios.  This  integer  specifies  the  number  of 
items  in  the  input  list  for  the  frequency  ratios. 

Word  4  (NCAL)  Number  of  speed  and  bearing  data  input  sets.  Each  set  of  data 
consists  of  a  speed  range  and  the  values  of  the  8  dynamic  coefficients  for  each 
bearing.  The  speed  range  is  specified  by  an  initial  speed,  a  final  speed  and 
a  speed  increment.  For  each  speed  value  in  the  speed  range  the  frequency  range 
is  scanned  and  the  corresponding  values  of  the  system  determinant  are  calculated. 
There  is  no  limitation  on  the  value  of  the  input  item. 

Wocvl  5  (NPST)  Pedestal  flexible/rigid.  If  this  integer  is  zero  the  program 
assumes  the  pedestals^to  be  rigid  and  no  pedestal  data  can  be  furnished.  If  the 
integer  is  1  the  pedestal  has  both  flexibility  and  damping  and  tfte  pedestal  data 
must  be  given* 

Word  6  (INP)  Input.  If  this  integer  is  zero  the  program  will  return  to  read  in 
a  new  set  of  input  upon  completion  of  the  computation.  For  the  last  set  of  input 
the  integer  should  be  1. 
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-VOHM-\~  (lXKli.6) 

2 

This  curd  a>nt  lins  i  single  word.  Youngs  modulus  E  in  lbs/in*-.  It  is  constant 
throughout  the  rotor.  Since  the  program  never  usesE  by  itself  but  always  in  the 
product  El  (1  -  cross-sectional  moment  of  inertia)  any  actual  variation  in  E 
can  be  absorbed  in  a  corresponding  change  of  I. 

Rotor  Data:  FORMAT  (4(1XE13.6)  ) 

The  rotor  data  consist  of  as  many  cards  as  there  are  mass  stations  (card  2, 
word  1).  Each  card  has  4  words: 

Word  1:  the  mass  at  the  station  in  lbs. 

Word  2:  the  Length  of  the  shaft  section  to  the  right  of  the  station  in  inches 

Word  3:  the  cross-sectional  moment  of  inertia  of  the  shaft  section  to  the  right 

~  ,  .  4 

of  tne  station  in  in  . 

Word  4:  the  polar  mass  moment  of  inertia  minus  the  transverse  mass  moment  of 
inertia,  lbs. in2. 

For  the  last  mass  station  the  shaft  lettgth  and  the  cross-sectional  moment  of 
inertia  has  no  meaning  and  may  be  set  equal  to  zero. 

LOCATION  OF  BEARING  SUPPORTS:  FORMAT  (10(1X14)  ) 

This  list  of  integers  provides  the  number  of  each  mass  station  at  which  there 
is  a  bearing.  The  stations  should  be  listed  in  sequence,  beginning  with  the 
lowest  number. 


PEDESTAL  DATA:  FORMAT  (6(1XE11.4)  ) 

If  the  pedestals  are  rigid  no  pedestal  data  are  required  and  item  5,  card  2, 
must  be  zero.  Otherwise,  the  data  for  each  pedestal  oust  be  provided.  There 
is  one  card  for  each  pedestal  and  they  should  be  In  the  same  sequence  as  the 
bearing  station  numbers  in  the  previous  list.  Each  card  contains  6  words: 

Word  1:  the  vibratory  mass  of  the  pedestal  for  motion  in  the  x-dlrection,  lb3. 

Word  2:  the  pedestal  stiffness  in  the  x-direction,  lbs/in 

Word  3:  the  pedestal  damping  coefficient  in  the  x-dlrection,  lba.sec/ln 

Word  4:  the  vibratory  mass  of  the  pedestal  for  motion  In  the  y-dlrectlon,  lbs. 

Word  5:  the  pedestal  stiffness  in  the  y-direction,  lbs/in 

Word  6:  -he  pedestal  damping  coefficient  In  the  y-dlrection,  lba.sec/ln 


LIST  OF  FREQUENCY  RATIO  VALUES:  FORMAT  (4(1XE13.6)  ) 

V 

This  input  list  gives  the  values  of  the  frequency  ratio  ~  at  which  the  program 

evaluates  the  system  determinant  (V  =  disturbance  frequency,  radians/sec,  tu  * 

angular  speed  of  rotor,  radians/sec).  The  values  should  be  given  in  descendir 

v 

order,  for  instance:  ~  =*  .55,  .50,  .49,  .43,  .40  .  The  program  automat 

inserts  in  the  list  the  "eigen-instability"  frequency  ratio  from  each  bearing 
determined  from  Eqs .  (64)  and  (65).  In  most  cases  these  values  are  equal  to  a 
ximately  .5  but  may  be  less  for  heavily  loaded  bearings.  Unfortunately,  the  " 
instability"  frequency  ratio  is  very  sensitive  to  even  small  deviations  in  the 
dynamic  bearing  coefficients  from  their  accurate  values.  It  is,  therefore,  re 
commended  that  the  irput  values  of  the  bearing  coefficients  be  checked  beforeh 
by  means  of  Eqs.  (64)  and  (65).  If  the  thus  calculated  frequency  ratio  value 
differs  much  from  .5  the  bearing  coefficients  should  be  checked. 

ROTOR  SPEED  AND  BEARING  DATA 

The  following  input  data  should  be  repeated  as  many  times  as  specified  by  word 
card  2.  First  comes  a  card  giving  the  speed  range  for  the  calculations.  The 
speed  range  is  defined  by  an  initial  speed  value,  a  final  speed  value  and  a  sp« 
increment,  all  in  RPM.  Thus,  if  the  initial  speed  is  given  as  3000  RPM,  the  fl 
speed  as  9000  RPM,  and  the  speed  increment  as  1000  RPM  calculations  are  perfon 
for  3000,  4000,  5000,  6000,  7000,  8000  and  9000  RPM. 

Next  follows  the  8  dynamic  coefficients  for  each  bearing.  There  are  4  values 

per  card,  hence,  there  are  2  cards  per  bearing.  The  first  card  gives  Kxx,  coCx 

K  and  o>C  ,  and  the  second  card  gives  K  ,  cuC  ,  K  and  cuC  .  All  the  co- 
xy  xy  yy  yy  yx  yx 

efficients  are  measured  in  lbs/inch.  The  cards  should  be  given  in  the  same  se¬ 
quence  as  the  bearing  station  numbers  in  the  previous  input  list- 

COMPUTER  OUTPUT 

An  example  of  the  output  is  included  in  Appendix  B. 

The  first  page  of  the  output  lists  the  input  values  for  checking  and  control  pu 
poses.  Next,  follows  the  results  of  the  calculations  with  the  results  for  each 
rotor  speed  given  separately.  The  first  line  specifies  the  particular  speed 
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mass 


and  then  foil  ws  the  "eigen"  -  instability  frequency  In  RPM  and  the  "eigen”  - 
in  lbs  for  each  bearing  ds  determined  from  F;s.  (64)  and  (65).  They  are  labeled: 
"INST. FREQ,  RPM"  and  "INST. WEIGHT" ,  respectively.  Thereafter  the  results  of  the 
calculations  for  each  frequency  are  listed  in  5  columns.  The  first  column,  labeled 
"FREQ. RAT.",  gives  the  frequency  ratio  ~.  The  second  column,  labeled  "DETERMINANT", 
gives  the  square  root  of  the  system  determinant  (i.e.  of  the  matrix  in  Eq .  (35)  ). 
The  third  and  the  fourth  columns,  labeled  "RE(DET)"  AND  "IM(DET) "  gives  the  real 
and  the  imaginary  pact  of  the  system  determinant,  respectively  (i.e.  and  A  .. 
fiom  Eq.  (56)).  The  last  column,  labeled  "ENERGY",  is  proportional  to  the  energy 
input  given  by  Eq .  (60). 

It  should  be  noted  that  in  order  to  determine  if  the  rotor  is  stable  or  unstable 
it  is  necessary  to  find  at  which  speed  it  becomes  unstable.  Hence,  results  must 
be  available  over  a  range  of  speeds.  To  determine  that  speed,  at  which  instability 
sets  ir,  the  results  must  be  plotted.  One  method  is  as  follows:  for  each  speed, 
plot  the  real  and  imaginary  part  of  the  system  determinant  against  the  frequency 
ratio.  Find  those  frequency  ratio  values  at  which  the  two  functions  become  zero, 
(there  are  usually  several  values).  Next,  plot  the  "zero-point"  frequency  ratios 
against  the  rotor  speed,  obtaining  a  curve  corresponding  to  the  Imaginary  part  of 
the  determinant.  Where  the  two  curves  intersect  is  the  threshold  speed.  Thus, 
for  the  output  example  given  in  Appendix  B  the  threshold  speed  is  found  to  be 
8,350  RPM  at  a  frequency  ratio  of  .4937. 
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FIG.  I  SHAFT  SECTION  BETWEEN  TWO  MASS  STATIONS 
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FIG.  2  CONVENTION  AND  NOMENCLATURE  FOR  ROTOR  CALCULATION 
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FIG.  3  BEARING  AND  PEDESTAL  SYSTEM  FOR  TRANSLATORY  MOTION 


55 


APPENDIX  a 

SAMPLE  CALCULATIONS  AND  INPUT  PON®  FOR  THE  COMPUTER  PROCRAM 
"UNBALANCE  RESPONSE  OF  A  ROTOR  IN  FLUID  FILM  BEARINGS" 
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i!  3MC  lot.)*  M94.xA7,NOOEC<*LlST 

C  MECHANICAL  TECHNOLOGY. INC.  JORGFN  W.  LUND  j^l* 

C  IUNiAlANCu  RESPONSE  UE  ROTOR  IN  NON-JNIEORK  OtAKlNGS  PNOGll  wti) 

OIMCNSION  d'-xruo.  UO I  •ciMXi.l  10.  60I.BMYCI10.  UGl.uMYSIlO,  dUl.VXC  ju3j 

llt0.41),»xJ(10.4ii.vfC(10.41).Y»SU0.41)  *DXU  10 .40)  .Das <  10. 4U)  .L'XC  jj<*  w 

2  (  10.4  0) .or J I  10.401 .AC ( IO.hC) .as  I 10.401  .VC  1 10.40  )  .rSl 10 .40 ) »dmx  A  1 4  C  jjbC 

3  )  .UMAtJI  40  I  .uMXC  I  40  I  •  IMXOI  40)  .DMYAI4U  )  *DMY3<  401  .DMYCI40)  ,DMYL( 4j  |  . 

4DvxA(40I.Oyxd(»0).L.XC(40).i5VXDl40)  .OVYAI  40  )  .  OV  Yb  ( 40 )  .0  VYC  ( 40  >  .  JO  70 

bDVYDI  401  .RM<  40 » .RL <  40) .RSI  40 I ,K IP  140  I  .HI T I  40  I .0 1 AI40 ) »D  1 3 < 40 ) ,l>1C (  uOoG 
6401 .0101401 .ANI40 ) -INI  401 .ON  140 ) .UX (40  I  .U Y 1 40 > • LU( 40 1  .LB (  25  1  .dXXX  I  uO >0 
73.2b) .dCxAl 3.?b I .PKx* I  3.2b) .bCXYl 3,2b)  .BAY Y I  3.2 5 ) .BCYY I  3.2b »  »BK yx  (  JlOO 

8  3  ..?bl  »BCYX(  3. i  .BjMXX  <  3,2b  )  .BDMXX  (3.2b)  .8SMXYI  3.2bi  .BCMXY  1  3.2b  )  •  0.110 

98SMYV I  3.2b  I .BOMYVI 1,2b : .RSMYXI 3.2b ) .BDMVX I  3 ,2b  I  C120 

DIMENSION  PMXl2bl.MMYi2bl.PAXl2>».PCX(2bl.PAYI2bl.KYI25l.P|Xl2bl.  0130 

1PJ Yl 2b  I .PSMX ( 2b  I  .PUMA  I  2b  I .PSMYI 2b  I ,°DMY I Ab  I  ,LC 120) .Ald.d)  ,b(  0.4  |  ,  j 1 40 

2C  I  8.4  I  .01  8. 1  I  ,F  (4.4 ,20)  .OM(  O.o  I  .ENT  110) ,  KHS (8.1  )  .DVUXI80)  ,  ulbO 

3  DYUYI80I .Dummy |  300  I .SmAJ 3  I .SHHI 3).SHC  I  3  >  »S*1D( 3  I  .AUXI2S .4 , 3  I «  u 1 60 

4B'JXI  2  b.  4. 3  )  .  OUM)  I  100)  .CL  NR  I  «)  0170 

CUMMON  A  •  a  ,  C  •  D  «  CFM  •  EN1  0200 

Common  RHS  ,  OuM 1  .  CF9  .  MAT  t  KFN  *  CLNH  0210 

COMMON  PRN3  0220 

COMMON  NS.NB.mU.NC.NPST.NMOM.NGYR.NCAL.NDIAG.INPUT.YM.sCF.RM.RL .  u230 

1RS.M  T.DGYR.RJP.ftl  I  .lB.L'J.UX.U “.LC. PMX.PKX.PCX. PMY.PKY, PC Y.KST,  0240 

2PIX,PSMX.P0i«X.P1Y,PSMY .POMY . SPS  T  .  SPF N  ,  SP 1 NC »6KX X  *RCXX .BKX Y ,BCX Y  >  02 bO 

30KYY.8CYY.8AYX  *BC  YX  .OSMXX  .BOM)' X  .  BSMX  Y  »  BDMX  Y  .  dSM  Y  Y  .80MYY  .  BSMYX  » 
4BDMYX.CF1K.CF1C.CF1U.CF lE.uF2<.CF2C.CF2D.CF2t.CF2A.CF2tf.CF2M,CF2N, 
5CFlA.CFia,CFlM,CFlN.CFl,C-2»CF3»CF4.CF9,CF6.CF/.CFB.STF.ANSP2. 
6ANSP.SPCAL.OVXA,OVXB.I)VXC,DVXO«DVYA.DVYB.DVYC»DVVO.DMXA.DMXB.DMXC.  028 


7DMxD,0MYA»D^YB«0MYC  «0M YD, AN .8N.DN. I | | 

COMMON  BMXC.aMXS.BMYCxBiMYStrfXC.VXb.VYCxVYS.DXC.DXS.DYC.DYS.XC.Xi, 
1  YC  .YS.Di  A.DIB.DIC.  •>  I  >.' •  OVUX  .  D  VUY  .DUMMY  .bHA.6Hd.8HC.8H0.AUX.BUX 


COMMON  JCAL.Amj.nITC 

201  111*1 

2 00  CALL  SUdR  0310 

C  ROTOR  CALCULATION  3000 

400  CFl.o.O  1010 

CF2"0.0  3020 

CF  3*0.0  303u 

CF4*0.0  3040 

CFb*0.0  30b0 

CF6*0.0  3060 

Cc7*0.0  3070 

Cc8*0.3  30o0 

N  t  TC •  1  30V0 

I  ST^  1  3100 

I  rN*9  3U0 

401  J  >  T  *  1  3120 

NCC* 1  3130 

IrlNC:  403.403.402  3140 

402  JFN*lC 111  JlbO 

GO  TO  404  3160 

403  JFN*NS  ii/u 

404  DO  40b  U*  t .NS  31o0 

DVUX(J)»0.0  31 Vo 

40b  DVUY I j ) *0*0  J2UU 

DO  4 22  |*|  T.JFN  3210 

<VT*|  32tO 

IFlJST-1)  411,411*406  3230 

406  <  F  N* JS  T* JK  T  32,0 

bMXCll .KFN)*0.0  329U 

bM*sil.KFN)*0,0  3260 


ri Mv C  (  I  »vl  I  =  0.0  it  .-J 

dWSI  I  \  I  -0.0  3200 

,)*f  {  T  •  J S  T  i  *0.0  U  tO 

/XSI  r  •  J4  I  1*0.0  3300 

„VC  (  I  .-1ST  )  *0.0  3310 

.GYM  !  .jST  I  =0.0  3320 

GO  10  ( <.0  <  .40“  .<.00*410  .41  8.4  la.x.  IB. *18  .416.418  1  .<ST  3330 

40  7  OXC ( 1  .  JS  r  1*1  .0  3340 

00  To  41d  33a0 

4"8  DXSI  2.JM  I  *1  •"  3  3e>0 

GO  TO  4 1  d  33?0 

409  >xct  1.JST  )  *  1 .0  33oO 

oO  TO  418  33x0 

410  DYSI4  .JST l*l.o  3400 

GO  TO  418  3410 

4  1  1  oMxCI I .1 ) *0.0  3420 

rt'TX.G (  I  .1  '  *0.0  3430 

O^VCI 1.11=0.0  •  •  3440 

d  ivsi  I  ,  1  |  *0.0  3450 

VXCU. 11*0.0  3460 

V  <91  I  .1  1*0.0  ,j  3470 

VYC' 1  .1 1*0.0  f  *  3400 

VYSII. 11*0.0  f  .  3490 

Xvfl. 11*0.0  ‘  3600 

xsu.  11*0.0  ‘  3510 

YCII. 11*0.0  ’  •  3520 

r S (  I  •  11*0.0  •  3530 

0X01.11=0.0  *  3540 

OX S I  I  .1  I *0.0  ■  ’  3550 

OVCIl.l 1*0.0  •  •  ■  3560 

t>YM  !  .1  1*0.0  3570 

00  TO  1407. 408.409.410. 412. 413.414, 413.416.4181 .KST  *  3580 

412  XCI4.1  1*1. C  3590 

GO  TO  4 1 d  •  3600 

4 1 X  XS(5.1)*1.0  3610 

GO  TO  4ia  3620 

414  YC<  7,  1  1*1.0  3630 

GO  TO  41a  3640 

4 1 5  YS(d, 11*1.0  3650 

GO  TO  418  36o0 

416  JO  417  J*1,NJ  3670 

<ST*UUIJI  3680 

DVUXIKSTl *UXC Jl»ANSP2  3690 

4 1  y  DVllvcxSTl  *  JY(  j:  .ANGP2  3700 

418  00  422  J.JM.JFN  3710 

<ST*J.J  3720 

XFN*«.ST-1  3730 

KMD*<  ST ♦ 1  3740 

1F( JST-ll  420. 420,419  ‘  3750 

414  IF(J-JST)  421  ,421.420  3760 

420  6MXCI  l  ,<5T  I  *t»V<C  I  I  .XFNI.OMXAI  J)«I;XCC  I  «  J  I  .OHXol  J I *LX8(  I  .  J). OMXCI  J>*  3770 

1 OYC  (  I  .  J  l  *OMXU{  J  I  *L‘V  S I  I.JHUI.IJI  3780 

a*xS<  1  .1ST  1  .O^XSI  ;  .4FM-0MXMI  J)»()XC«  I  »JJ.OMx«(  J1  »DX!>I  1  .JI-OMXOI  J)»  3790 

10  YCII  »J)4JMXC|  J>»t)VSI  1  »  J )  ♦•  l  r)  1  J I  3800 

1MVC(  1  «<5  T  I  *-3MYC  '  1  .xFM.lMY^I  J|  »  JXCIt  ,  J 1  *OMYu I J I • JXS I  l  ,J1  .OHYmI  j|*  3610 

1CYCI  I  .J  j  ,  *  >f',|  |  »  J )  *0 IC  ( J I  3820 

84»S(  I  ,<S  r  1  *-;VY  ,  I  1  .<FN1-0»YI>|  JI#UXC  (  1  «  J  I  ♦OV  VC  I  J  I  «i'XS  I  I  .JI-OMVtM  Jl4  3o30 

II.  'Cl  1  ,  J  I  *  )«Yi  I  J|  * ,5Y  M  !  ,J  I  .I'l  J1  3460 

WXCI I ,J*1  l*vxri 1 .J». )VXA( J1*XCI I .JI-DVXrll J»*X6( l . J 1 -UVXC I J ) • YC ( 1  * J  3890 

1  I  -OVXDI  J  1  •  Y  > I  I  ,  J  1  IX  C  J  )  *  3860 

vxs«  I  ,j*i  »«vx*  1 1  ,ji  1  jvxmjinui  .ji.uvxai  ji*xsi  1  .ji.jvxoi  j:  *rtU.j  jm/o 

1 »-OWXC{JI*Ydl I..M-JVUY IJI  *  3880 

V  *c I  I . J*1  I  * VYf  C I .  JI-OVYCI Jl«xci  1 tJI-0WYnl.i|»X5l  1 .JI.UVYAI J)»TCI I.J  3«90 

J 


t  l-OV¥J(  JI<YSI  !  .JUOVUYIJI  390  "i 

VYSII. J<1I<VYS<I  » J»<0VY0l Jl <XCI  I  .Jl -DVYCI Jl<XS<  I  *J  l<OVYB< J>*VL  UiJ  391  C 

iuovvai ji »ysi i .j'+nvum ji  3920 

IMJFN-JI  *22, *22, *21  3930 

*21  *JW*Clt  .<MOI»9yi*C«  1  .XST  |<RlI  Jl  »V*U  !  .J<ll  39*0 

Byxs(  I  .xmoi  •ewxsi  i  ,xst i*rl<  ji<vxsi  i  .j<i  i  393  o 

&MY  C  I  I  .X'ADUdMVCI  t  .XST  l<RL<  JI*WYC«  1  ,J<1  I  3960 

BMySI  I  t**OI«8MYS(  I  •  <  St  l*RL<  J»»VYS  I  1  »J*1 1  3970 

OXC  I  I  » J*1  I  aOXC  I  t  .JI<ANIJI<BMXC(  1  .X S Tt <bN < J I • VXC < I  .  J<1 1  39o0 

OXS(  |  » J*1  I  <UXS(  I  •  JUAN<JI<BMXS<  I  *XST  l*aN<  JI<VXSI  I  .J<ll  3990 

OYCII .J<1  I  ■  JVC  I  I  •  J  >  <AN  <  J  I  »!1VYCI  I  »XSTI<dNlJI<VYC<  I.J<U  *000 

OYSI  |  »J<1  l-OVSI  I.JUANI  JI*8MYS<  I  .XSTI<BN<JI*VYSI  I  *J<ll  *010 

XC<  I  » J*1  ) »XC (  1  »J !♦«<.< JI<0XC I  1 • J I ♦BN  <  J  I  <6MXC< 1 .AST I  <UN<JI<YXC  I  1 .J<1  *020 

1)  *030 

XS(  I  ,  J<1 1  »XS< | i J I <Hl I J) *0XSI 1 . J I <BN I J I  <BMKS< l  .AST  I <DN(  J I <VXS  < I  * J*1  *0*0 

1,  *030 

YC  <  I  •  J<1  I  <YC  ( |  •  J1<RU<J1<DYC  I  I  .JUBNt  J I  <BMYC  1 1  .XST  1  <DNU>*VTC  I  I  .J<1  *060 

1,  *070 

YSII»J<1)»Y3<I <  J ) <HL ( J ) <0 Y  S II • J I <BN I J I  <BMVS< I  .AST  1 <ON<JI<VYS<  I  »  J<1  *0«0 

1 |  *090 

*22  CONTINUE  *1°° 

I F I  NO  I  *G I  *23. *2*. *23  *130 

C  DI AGNOST I C  2  *120 

*23  WHITE  <*.1361  *130 

WHITE  <6.1161 |St.IFN.JST.jFN»XST.SFN*XMO.NCC  *1*0 

WHITE  <6. l**l< <0*0  I .JI.OXSI 1 ,J! *DVCI  I  »JI »OYS( t .Jl  .  XCtI.J),XS<  *130 

11  .Jl.YCI I . J I *YS  < I .Jl. 1*1. 101 »J*1.NS|  *1*0 

WHITE  <6.1*3)<OfA(J|.01H<JI.OICIJI.OIDIJI .OVUX I J  I  .  OVUTIJI.J-1  *170 

l.NSI  *1W 

KST<NS<1  *190 

W  1 1 TE  <6.10*l<<VXCI|.JI.VXS<I.JI.VYCII.Jl.VYSII.JI.  I-1.10I.J-I  *200 

l.XSTl  *210 

X  j  T <NS<NS  *220 

WHITE  1 6. 10* I < ( dMXL ( I . Jl .BMXS  < I . J I »8MVC I  I  »JI  .BUYS  I  I.JI»I<1»10|.J<1  *230 

l.XSTI  *2*0 

*2*  lr<N.3-JFNI  ***.***.*23  *230 

C  ROTOR  HAS  COUPLING  STATIONS  *2*0 

*23  XS T ■  JFN*,lFN  *270 

IFMST-IOI  *2*. *3*. *3#  *2*0 

*26  00  *27  J<1,*  *290 

F< l.J,NCCI»dMxC< J.XST1  4300 

F 1 2 • J .NCC I <dMxS  <  J.XST I  *310 

F<3.J.NCCI<aMrCIJ.XSTl  *320 

F  I* ,  J »NCC I <dMYS <  J.XST I  *330 

XFN<J<*  *3*0 

B<1.JI<-8MXCIKFN.X5T|  *330 

a<2.JI«-aMXS<xFN.XSTl  *3*0 

a ( 3 . j i •-bmyc  <  xfn.xst  r  **io 

B  <  * , J I «-8MYSIkFN.X ST  I  *3»0 

DO  *27  |<1>*  *390 

A<|,J»-F<I,J,nCC»  **oo 

*27  C< I  ,JI*F< I.J.NCCI  **10 

CF9<5CF  **20 

MAT<1  **30 

IF < NO I AG I  *29, *30. *29  ***0 

C  DIAGNOSTIC  3  ***0 

*29  WRITE  <6.1371  **6J 

•RITE  <*.10*MICII.JI.I<l»*».J<l»*l.l<8ll.Jl.l<l.*I.J<l.*l.HA<  I.J  **70 
II  .I<1.*I.J<1.*I.< <F| I.J.NCCI .|<1.*I.J<1.*I  ***0 

**0  XFN<XFN  **90 

CALL  COS  *300 

GO  TO  <*31.Slfl.311),«FN  *310 

*31  Dl’.l!<-aMXC<9.<ST»  *»20 


'(2.1  I  =  -:)■«*  I  4510 

i.  M  * !  I  =-  )«rc  ( -».<  5f  i  4540 

V'.il  4530 

c  ~  C  ^  4560 

■  <♦  5  / 0 

*  *‘i  =  «.cN  45o0 

cm  F  25  45*0 

C?  to  H>WOl5oUlt<fN  **600 

4  12  nj  4  3  1  t  »  1  .4  4610 

4' 1  A  (  J  ,6  I  «C(  :  .  1  I  4620 

oo  434  j.jst.jfn  46J0 

00  4  14  <  =  1»5  4640 

<sr*>c»4  46S0 

00  434  1=1.4  4660 

VXC<<St.J*l)*VXCUi>T,J4l  )4VX'C  I  I  .J*n*At  1  .<)  46  10 

VXSIKST  » J*1  )  *VXSI K6l  .J+1  l*VXS( I .J*l )«A(  I  .K  I  4660 

VYCIKSf  »J*1  »  «VYC  IK.. it  ,  J4 1  I  ♦  VY  C  (  1  »  J*1  1*4  I  1  »K  |  4650 

VYS KSt  .J»l  I  *VYS UiT  .J*l  (♦VYill ,J*1 l*AI I «K)  4700 

XC(<ST . J) *XC UST . Ji*XC  1  I .JI*A( I »< I  4710 

XS1XST  «JI*A5(<St » J I 4XS  t 1  »  J I  •  A  ( 1  »<  I  4  720 

YCtKS  t  ,*)  *YC  (<:>r  .JI4YC  t  I .JI»A(  I  .M  4710 

YStKST.JI *YM«ST . J! ♦YSC I .J»*A( I ,K »  4740 

OXCf <ST ,J>  *0xr ( < jT .Jl 4DXCI 1 «J)*A{ I «K I  4  750 

DXSUST  .  J|*DXSt<Sf  .  JIaDXSI I  .JI*A(  J  ,<  I  4760 

OYcncsr. ji«DYctxsT.j) 40yci i »j>*Afi .< i  4770 

434  DYStKST.J  I*0yM<ST.J)40YS(  I  .JI»A(  I  ,IC)  4760 

XFN*JFN4JFN  4750 

xi*D=jsr-»jsr  4600 

DO  436  J*<M0.<FN  4610 

00  436  <*1.6  4820 

(CST»<*4  4830 

00  436  1*1.4  4840 

HMXC ( X  S  T . J I *BmxC ( Kb  7 . J I *8MXC  t I ♦ J I *A I 1.0  4560 

8MXS(<iT. JI*BMXSt<5T  .  J  i  +-9MXS  (  I  . J I *A  (  |  ,M  4860 

HMYCl <ST. Jl *e«YC(X6r .JI*BMVC( | ,< I  4870 

436  RKYS  (  <ST  *  J I  •BMY'S  1 XST  »  J  I  ♦RMYSt  l.JIMtl.K)  4880 

00  to  4 h2  4850 

438  D( 1.1 I *-dMXC  < 10  »<ST  I  4900 

DI2.U*-0MXS(  10.1ST  1  4910 

01 3.1 I *-9MYC ( 10  »<ST )  .  49P0 

0(4,1 )*“9MYS! t O.XSt |  4930 

00  428  J* 1  * 4  4940 

0)  428  1*1,4  4930 

428  CtI.JI*F( I.J.NCCI  4960 

C  r9*SCF  49/0 

M4  T*  3  4960 

KrN*<FN  4990 

CALL  EOS  6000 

GO  TO  1439*510*511) *KFN  5010 

439  DO  440  J-JST.JFN  5020 

00  440  1*1.4  5030 

VXC110.J4  1  l*VXC(  10.J41  l>VXC(  l,J*l  l*CU.ll  6040 

VXS(10.J4ll«VXS(10.J4l 14VXS1 l.JM l*C( I* 1 1  5050 

VYC  1 1 0  *  J*  1  I  *VyC  1 10»  J*1  I  ♦  VYC  1  1.J41  !•<.«  I  .11  5060 

vvsno.j4i»*vYSiio.j4i»4wsii.j+ii6Cci.n  5070 

XCI10.J)»XC(10.JI*XC< I »J)»CI 1,1 J  5060 

XS(10.JI*Xi(l0.Jl4XS(l ,J)*CI I.ll  5050 

YCl 10.J)*YC( 10.JI4YC1 I *JI»Ct 1.1 1  6100 

VS(10»J1*YS( ln.JItYM I ,JI*Ct 1.1 1  8110 

DXC(10.JI*DXC(10.JI*0XCI I.J1*CI 1*11  5120 

DXSC10.J 1*0X5(10. Jl40XStt.Jt6Cll.il  6130 

OYCtlO.JI*OYCtlO.Jl4nvC(l.JI*CII»ll  *  5140 

440  OYSt 10. Jl *0Y5( 1 0 . J I *OYS I  I » J I *C I  1 .1 1  \  5180 


<f S*JFN*JFS 
*Mn»JST.JST 
00  <.41  J.KMD.KFN 
00  441  1 •  1 .4 

flM*C ( 10.  / 1 «BMXC ( I0«J1 ♦BM*C II  ,J|»CI I  .1  I 
RMxK  I  I  0  ♦  J I  *HMX S  (  10*  Jl  .RMyKI  I  *j|  »C  (  !  .  1  ) 
0«vr  (  10.  JMrfMvc  f  10.JS.Rmyc  1 1 .Jl *C<  I  .1 ) 

441  3*751  1  0  *  J  »  *  (  1m.JI.RMy SI  1  .J|«C(  t  .1  1 

442  JST*J*N 
NCCOCCM 

IF(SCC-NC)  444.444.443 

443  JFN*NS 

GO  TO  444 

444  JCN*LC I  NCC  I 
443  GO  TO  404 

C  END  CONDITIONS 

446  <  ST  *NS*NS 
x fN»NS*l 


D  )  447  J 

*1.8 

Ccw(l 

.Jl 

•dMXC( J.KST  | 

Cr'M  1 2 

.Jl 

•BMXSl J.KST j 

CPMI  3 

•  Jl 

•OMVCl J.KST  | 

CPMI  4 

•  Jl 

■dMVSI J,<sr» 

CPMI  5 

.Jl 

•VXCl J.KFNI 

CFM(  6 

.Jl 

■VXSI J.<PNI 

CPMI  7 

•  Jl 

•VYCI J.<FNI 

447  CFMI8 

.Jl 

•WYSI J.KFNI 

RHS(  1 

.11 

•-BMXCI9.KSTI 

RHSI2 

.11 

■-BMXSI9.KSTI 

RHSU 

.11 

■-3MYCI9.KSTI 

RHSI4 

.1  1 

•-8MVSI9.KSTI 

RHSI5 

.11 

*-VXC(9,KFN| 

RMS  ( 6 

.11 

*“VXS| 9.KFN | 

RMSI  7 

•  1  ' 

•-VYCI9.KPNI 

RMS  1 8 

.11 

•-VYSI9.KFNI 

IP  (  IS 

T-101  449.448.448 

448  RMS ( 1 

•  11 

*RmS 1 1  .1 l-BMXCUO.KSTI 

RMSI2 

•  11 

•RMSI2.1 1 -dMXS ( 10.KST  1 

RMSI  3 

•  11 

•RMSI3.1 1-dMYCIlO.KSTI 

RMSI  4 

•  11 

’RMSI 4.1 l-JMYSI 10.XSTI 

RHSI5 

.11 

■RMSI3.1 l-VXCI 10. KPN) 

RMSI  6 

.11 

■RMSI 6«ll -VXSI 10. KPN  1 

RMS  1  7 

.11 

•RHSI7.1I-VYCI10.KFN) 

RMSI  8 

.1  1 

•RHSI8.1 1-VYSI 10. KPN) 

440  CF9-SCF 
MAT  *4 

IP  I  NO ! AG  1  4S0.4S1.4S0 


c  DIAGNOSTIC  4 

430  WRITE  (6.1361 

WRITE  (6. 10411 1  Cl  I »J)» I ■ 1.4 1 ,J*1 .4) »( ( A | I .Jl .!■! .4) • J«|.» 1 f ICFMI I • 
1J 1.1*1 .6 1 .J* 1,8). (RMSI 1, 11.1 -1.6  1. CP  9 
WRITE  I6.144M  IOXCI  I  »JI .0X51  |  .J|  .OYCI I  »J!  ,0Y3I  I  .J»  .  XCI  I .JI.XSI 
1 1  .Jl .YC( I *JI .V$( l .J| .1 *1,101 .J*l .NS I 

431  <*N*<FN 
CXLL  EOS 

GO  TO  (452, 310, 3111. KPN 

432  E5T 191*1.0 

lc  (NO  I  AG  I  496,497*426 
<  DIAGNOSTIC  3 

496  WRITE  (6.1391 

WRITE  16. 144I| (CPMI I .Jl, 1*1 ,61 .J*l.i) 

497  IP( IST-10 I  453.500.300 

433  IPINCAC-ll  474.473.474 


5160 
31/0 
5160 
5  19C 
5200 
52  :o 
3220 
5230 
5240 
5<  50 
5260 
52/0 
3260 
5290 
5300 
3310 
5320 
5330 
5340 
5350 
5360 
5370 
5360 
5390 
5400 
5410 
5420 
5430 
5440 
3450 
5460 
5470 
5460 
5490 
5500 
5510 
5520 
5530 
3540 
5550 
5560 
5570 
5560 
5390 
5600 
5610 
5620 
5630 
5640 
5650 
3660 
5670 
5660 
5690 
3700 
5710 
3720 
5730 
3740 
5730 
3760 
5770 
5760 


••  »  <  a  a  !  I-.  I  6.  I  ,t  I  '  PC  At 

.  * fc  «•  \ n  :  01-0.0 

!  r  (  \ '  v  }  I  tSiiii  ■"?  ,484 
<.■<.  aw i r “  if.un 

H  v*  I  ’ 

;-C  T'  4  S  6 
..*•  i  <v-»=l 

c  *-?!U  :jf->L’r 

446  » £?  :  T  ?  (6,1-11 
aR|TE  (6.180) 

^0  447  I  =  1  •  -S 
4-1  r  ENT  (  I  I  *CFM<  I  .1  1 

iriNCi  4->y,46o«-4''!< 

488  <FN*LC  1  NC  I  ♦  1 

no  *  S 1  1  *  1  *  •« 

4  9  ENT ( I  I =0.0 
GO  TO  461 

460  <FN*NC*1 

461  00  716  j=i»m; 

VND=LUI J) 

OVI.IX  (  VN3  I  ='.;x  (  j)  *  4  N  8  P  2 
716  0VUY('*N0)  »J>  .  Jl  *ANGP? 

.01*5*0.0 

ENGY*0.0 

9NB»1 

V3R.1 

0190*1.3(11 
L  J3L*LU( 1  I 
O')'  470  J*l.NS 
I  "  J  J-<  FN  I  464,462.64', 

462  PD  463  1*1.4 

463  ENT  (  I  )*CFM(  I  .1  ) 

464  <jT  *J.J 
<C=<ST-1 

TO  466  <*1.3 
00  466  1*1.4 

466  H ( I  ,< 1 *0.0 

00  466  I  *  l  »  J  F  *4 

HI  It! 1 *«l  1 •! I.XCt ! .JI«INTI 1  ) 

rt  ( 2 .1 I  *d ( 2  » 1 )  *X*  ( 1  ♦ Jl *rVT f  t  ) 

31 3.11 *B( 3.1  I .YC< I . Jl »ENT( 1  ) 

(4.1  l*9(4.ll*Y>(  I  ,J)»£NT(  1  ) 

•’(1.21*3(1  .2  l+JMXCI  l.<C)*rNT(ll 
0(2»2)*3(2«2).<3MXS( I  .<C)*CNT(I; 

:1  (3.  21*3(3. 21  .avYCI  I  »KCl»E\f  (  II 
H (4.21*3(4. 2) »dMYS( 1 ,<Cl*ENT: 1 1 
0(1.31*3(1.31 .oMXCf I . <ST ) *r  NT ( I  I 
3(2.31*3(2.3)  ,JX*M  I  ,K5T  )»i  NT  *  I  ) 

3( 3.1 )*d( 1.3).b>YC( l »<3T l*ENT I  1  I  ' 

466  0(4.31*3(4. 3l.fi'1Y5(  I  ,<6T)*‘NM1  I 
1 F ( J-IBRG  I  70*. 703. 700 
C  TRANS'*  1  T  T  EO  FORCE.  PEJECTAL  Vl.’TICN 

7P3  CF1A.RV(J|»AN'-P? 

SMA(1  )«CF1A»H(  1.1  ’♦r'VUX(  Jl 
SH<3  ( 1  3 *CF  1  A*4 (2.1  ) “OV'JY (  j • 

SHC  I  l  I  *CP  1 A •  fl 4  3,1  I.OVDYI  J) 

SHO  ( !  1-CF  14*ii (4,1  )>C7UX(  J) 

00  704  1*1, IFN 

SHA  (  1  )  *ShA  (  1  I.lvxcn  .Jl-vxci  i  .J»l  1 1*ENT  1 1  I 

SH3(  1  l*SH->(  1UI/XM1  »J1-VX.(  !  .J»l  1  I*IKTI1  » 

SKcm*5HCi  1 i«i  wei  i.ji-vvgi  :..*n»»i  nth  1 

’04  SHO  ( 1  I  •<<')(  1  )  ♦(  VYS(  1  .JI-VVSI  1  ,J.l  #  I  •»  NT  1 1  I 
IF  INPSTl  706.706.706 


3  l~t‘. 

sboc 

^810 

3820 

6830 

8040 

8830 

8860 

8870 

88OU 

889Q 

8900 
3910 
3920 
39  30 
8940 
6960 
3960 
8970 
8980 
3990 
6000 
6010 
6020 
6u3G 
6040 
6030 
6060 
6070 

6OB0 
6090 
6100 
61  10 
6120 

61  30 
6140 
0180 
Albo 

0l  70 
6160 
6190 
8200 
6210 
6220 
6230 
6240 
6230 
6260 

62  7  o 
62o0 
6290 
6300 
6310 
6320 
6  330 
6340 
6330 
6360 
6370 
6380 
6390 
6400 
64  *0 


’**!  (  "  !«■  ji'vk  (  |  /  64cO 

C‘*l  V.D-/V(  Wjni/ias.OM^AN'P.’/Ay.S  64  30 

C  •'  l  A  *.-><x  (  Vrt'M /AV$-rc  1<  644C. 

r  y;jM  i /Aws-tt  h  6430 

;e  1  '  'PC*  (  *•!.«  ) /AV$«AN3P  6466 

Cr  l:>ii’C  Y  (  vh  <| /A‘4&*ANSP  64 /l/ 

CHtrtM<*C1U»CH‘-»0  1C  bboJ 

*fF  1 A-SHJI  1  )  »(.P  1C  l/CF  IS  6490 

SMBI  3  |  >  I  ->>4A  ( 1 1  »CPK+SM3<  1  )  *»>■  1A1/I.PU  6600 

CMr:*Cf  1  l»CFlR«CPll>»Cf  ID  6310 

SHCMI*li-4Cm*CFJtl-SHCC  1  1*0  ID l/CF It  6620 

sni»  )  «  (SHCI  l|»CFl')*SHni  1  !*CF!d>/me  6530 

SHA(.?|».<HA(  1  UCFU»SHAOl  6640 

AHH  { 2  I  *  Vl'J  ( 1  )  *CF  1<  »SH9  I  J  |  6630 

SHC(2l*SHcm*CFiy»6HC(3l  6^60 

SHDI2  l«SMDt  :  l*CFlM»SHf)r‘ 1  6670 

CFlA.mi#n-5HA»3»  6500 

CF1b«B(2i1 I-SH813I  6660 

CFl<«niJ*ll-SHC(?»  6600 

CFlV»d(4,l)-SH0«3|  6610 

GO  TO  707  6620 

706  CFlA*B(l»ll  6630 

CF13*:3(2,  1 1  6640 

CflK«fl(}ill  6660 

CFlv,fl(4,n  6660 

SHAI2  l-SHAd  |  6670 

5*0(21 -SM31 1 1  6660 

SMCI2I*3MC(1I  6690 

SHO(2l»5Hn(ll  6700 

SMA ( 3 ( *0,0  6710 

SHB«3)*0.0  '  6720 

SMCI3M0.0  6730 

SHD(3I*3.0  6740 

CF1C-0.0  6730 

CF13«0.0  6760 

7C7  CF2A.-1ICXX  (  1  «M"R)*BC5XI  ?,MHR  I  *ANSP*><CXX  I  3  »  v-!R  I ‘ANJP?  6770 

CF?.1»«CXV(  1  ,MuRI-»OCXY|  2»F'dai«A,V‘>P4»  CXYC  »  »9;  ■<  1  *AUW>2  67»C 

CC2C »HCYX  (  1  »»fim*jCV*<2.(YSk»»A!<SP4eCVXl'»»v;»^|*AN<P2  6790 

Cc2l)»BCVV  I  1  ,MPR  I  ♦dCVV<  2,4-JH  !  •  AN3P*t*CYY  (  7  , AB*  l*ANSP2  6800 

CT2M«9KXY  ( 1  ».VBR  I  *B*X  V (  2  ,»Ml*  I  •A.'jFP^mx.XV  I  '  t^HK  I  »ANSP2  681  3 

CC2N»BXYX  (  1  .MBR  I  *BXY x  (  2  ,vqtf  )  »A35,P4B*yx  ( 3 * VBM I • AHSP2  6820 

OISSOIS»*3«l4l6927»(CF2A«(U  IA*CF  1-4CF  1  l0)4CF2O»tCFl*.*CFl**  6830 

lC'lW»CFlM)*(CF2iJ4Crx»l*(CF|M»CFl4*C7  U»CF  IP  »♦  C02H-LF2K » *  I CF  lo*  68*0 

2CF)K-CF1A«C1 ly|*GFlt*( SHA (31 *SMm (  3  I  ♦  SHb  131  *6HB 1311 ♦CFIO*  (  6HC 131*  68J0 

3S.<»3l*SHD(3l«3HOOm  6860 

IF  IN909I  717,716,717  6870 

717  CFIA.0,0  6660 

CFld»0#0  6830 

CF1K«0.0  6900 

CFlv»0»0  6910 

00  740  I » 1 , l F V  6920 

CFlA.CFlA*l5XCll,J»*tNTHl  6930 

CF18»CCH*0X.61  (  ,J1*E7|T  I  1  I  6940 

CU<»CF!K*DYCll.JI*fcNTU>  6930 

740  CFiy-CFiy^pySlI.JlFCNTm  6960 

CF2A»flCl,31-B(l,21-CIAI.’|  6970 

CF2a»8l2«3)-JI2,2»-018U»  69»0 

CF2C«813#?l-:J(3»2l*’OIC  l<l  6990 

CF20«B14,1l-d(4»2l-0ir)lJ»  7000 

IF  INPSTl  741,742,741  V  7010 

741  CF2\UlPS'UC9U<n-P!XI'«BR»/186.CM«Ay,32l/A*S  7020 

CF2N.|PSvriKM»l-P(y (9rtRI/3»6,069*A\SP2)/AyS  7030 

CFK«P™x(W!j9i/«KS*A9$P  7060 


CFl  =P')Vf'(MdR|/4vs»4NSO  7U50 

CF  lt»Cf  2M«Cf- 1C*CF  1C  7C6u 

CF3A*(CF2A*CF2M-CF2B*CF1CI/CF1£  7070 

CF  3M*  (CF2  A»CF  K  +  <-F2  i »CF2M | /CF  IE  70*0 

CF1E*CF2N»CF2m*CF1D»CF10  7090 

CF3C*(CF2C»CF2',-CF20*CFlf;i/CFl£  7100 

CF3r)«(Cr2C»CF10  +  CF2D»CF2N)/CFl£  7110 

CF1A.CF1A-CF34  7120 

CF10«Cc18-CF3R  7130 

CFlH*CFlK-CF3f  7140 

C"1M«CF1M-CF30  7150 

Ce2M«CF iC*(CF3A*CF3A*CF3P*CF3B)+CF10*  1  CF3C*CF3C4CF3D*CF3D )  7160 

r,  3  TO  743  71  ?y 

7<»2  CC2M«0.0  7160 

743  0ISS»DlS$+3.l4l5927*(CF2A»CFia-CF2b*CFlA+CF2C*CFlM-tF2D*CFHt+CF2M»  7190 

716  03  708  1*1*3  7200 

CF1A*SHA(I)  7210 

CrlB*-SHB (  I  )  7220 

CF 1C*SHC 1  I  I  7230 

CF 1D*SHD( ( )  7240 

AUX(MBR,l,II»<;ORT(CFlA#CFlA+CFlo»CFlB)  7250 

AUX1MBR. 2. I )*ANG<CF1A. CF 13)  7260 

AUX1M8R.3.1  )*SQRMCF1C*CF1C4CF1D«CF1D>  7270 

AUX1M8R.4.I )*ANG«CFU),Cf  1C)  7260 

CF2A«CF1A*CF1A  7290 

CF2B-CF18*CF10  7300 

CF2C«CF1C*CF1C  7310 

CF2D«CF1D*CF10  7320 

CFlE-(CF2A*CF28-CF2C>CF2C)/2.0  7330 

CFllC«(CF2A*CF28-CF2C-CF20l/2.0  7340 

CF1M«CF1A*CF1C-CF18*CF1D  7350 

CF1N«CF1A*CF1B-CF1C*CF1D  7360 

CF2A»(CF2A-CF20*tF2C-CF2OI/2.O  7370 

CF2B«S0RT(CF1ic*CF1<4CF1M*CFIHI  7360 

CF3B-CFlA*CF10^CFld*CFlC  7390 

CF38*CF3B/A8S(CF3fl I  "  '  7400 

BUXIHBR.1.I)«S0RT<CF1E*CF28)  7410 

BUX1MBR.2.I >»CF3&*SQRT<CFl£-CF2dl  7420 

BUX1MBR.3.I )«ANG1CF1<,CF 1MI/2.0  7430 

708  BUX(*<BR.4.I  )«ANG(CF2A,CFlN)/2.0  7440 

IF1MBR-NBI  702. 701*701  7450 

701  LBRG*NS*2  7460 

GO  TO  709  7470 

702  MBR*MBR4l  7480 

LBRG-LB1MBR)  7490 

709  IF  ( J-LN01 )  714,710.714  7500 

710  IF  1MNB-NUI  712.711,711  7510 

711  CN8L«NS*2  7520 

CO  TO  713  •  7530 

712  *»NB«MNB*1  7540 

LN8L-LUIMNBI  7550 

713  £'4GY*ENGY  +  3. 141 592  7* 1 DVLX  (J)*(B(2. 11-813.1)  1*0VUV1 JJ  *IIM  1.1)4  .  7560 

1814,1  III  7570 

C  CJ8VERT  RESULTS  TO  FLLM»SIS  7580 

714  0)  469  1*1*3  ,  7590 

CF1A-BU. 11*811. 1)  7800 

CrlB*812»l )*BI2»I  )  7810 

CF1C*B(3. 1 1*813.1 )  .  7820 

CFl 0*814.1 l*BI*»l  )  .  7830 

cfie*icfia4cfib4Cfic4Cfio)/2.o  78*0 

CFl<.KFlA4CFia-lFlC-CFlP|/2.0  *  *  7*»0 

err'-1!  !  .1  I.BI2.I  >*rt(4.ll  '  7880 


CM  ill  -  -  I  •«  •!  )•  lUdl  *»*(  1,1  (•3  14.,;  > 

cr  \r  *s.:rt  i  cf  i<»cfu*cf 

CF JM»F| 1  .  I  Mrtl4.  I  I -rt*  2,1  I  Ml  7.1  I 
Cr  3 J 1 CF ’  I/A-3SICP1M 

■tit.!  ms)rt(cmf*-:mci 

■»(?.!  i-cc  3d»'-.-:'»T(cnr-cnr. i 

46  7  Hl',1  J  1 1 

468  414.1 MA\G(CM4»CF|dl/2.0 

469  CONTI  VJF 

WR  I  TE  16.  1 48  )  J ,d I  1.1  1.412.1  )  »b  (  3.1  I  <3(4,1  I  .0(1,21,  B(2«2  >  «<>(3. 

121,814,21 

<>70  w  R  I  I  E  16*14913(1,31, 612, 3). d(3, 3), U(4, 3) 

CF 1 A« AMS»AV5 

dissoiss/cfia 

FNGY.FNGY/AMS 
WR I  TF  (6. not 
00  726  Ml,? 

IE  11-21  723,720,  722 

720  IE  (NPSTI  721,726.721 

721  w  R  I  T  c  (6,731) 

00  TO  723 

722  WHITE  16. 732) 

721  WRITE  (6.7331 

eO  724  J.l.NB 
IR*L 3  (  J I 

724  WRITE  (6.  148  |m3R«HUX(  J.i.D  ,.«OXI  J.2. 1  I  .PJXIJ.3.I)  »  BOX!  J.4,  l  )  ,AU 
lXU.l.I).  VJXIJ.2.H.AUXI  J.I.II.AUXI  J.4.I) 

726  CJNTINUE 

726  WRITE  16. 734JCNGY.JISS 
GO  TO  C611. 480,613), <V0 

C  MAKE  RCAOY  FOR  GYROSCOPIC  MOVENT  calculation 

480  1ST. 10 
IFN.iO 

WRITE  16.1631 

C  CALCULATE  GYROSCOPIC  moment 

481  no  482  1*1.8 
482  ENTC I )»CEV( 1.1) 

I E  |  NC )  48  7,436 .463 
483  <EN«lC!NC) 

DO  484  1*1.4 
484  ENTIIMO.O 
GO  TO  496 
48S  <FN*NS*1 
486  DO  4R3  J* 1 .NS 

IEI J-KFNI  409,487,489 
467  00  488  1*1,4 
488  ENTMMCFMII.il 
489  CE1A.0.0 
CFIB.0.0 
CE1C.0.0 
CE1C0.0 
00  490  MI.IEN 
CFIA.CF1A.0XCI  1  r  JI*ENTM  ) 

1  CE19«CMJ.0XSM  .JMENTII  1 
CF1C*CE1C.0YCM.J)*ENT  1 1 1 
490  CElD.CFn.’'Y«M.JMENT  II  ) 

CE 1 E  *CP 1 A*CF lo 
CE  1K«CE ld-CF  lc 
CF IM.CE1 A. CF10-CFl0.CE  1C 
CE1N«Cc1E»C-'1E.CFI<8CF  U 
I F | CF INI  492*491.492 
491  guijmq.o 
o*BT  JMO.O 


tb  o  '. 
74 

770- 
7  710 
77,  - 

n  )<j 
7  /•.(. 
/  7  6  x 
77&U 
7  7  7'j 

7780 
77^0 
/OOu 
7810 
7620 
7o30 
7840 
7860 
7860 
78  70 
7boO 
7890 
79oo 
7910 
7920 
7930 
7940 
7960 
796i'. 
7970 
79o0 
7990 
oouu 
8010 
oo2u 
80  30 
0040 
oOsO 
8060 
8070 
8080 
o090 
8100 
6110 
8120 
0 130 
0140 
6l?0 
8160 
0170 
6100 
6190 
6200 
o2 1C 
022C 
o2  JO 
6240 
0200 
0260 
\«2/C 
6200 
6290 
6300 


•'  I r  (  j  i  o.?  a  310 

"I'M  .1 -0.0  0320 

GO  TO  49?  83 JC 

wii  cf  i •■  =  ;c  /cr  i *■  •  0340 

CF1*.*-*!3 1.M«A\S?2  3  3oG 

cf i •  =or  •  \’»cf  i -*rr  i v  836o 

cr  is-.>  i N « c -  1<»C- lv  a3/o 

Cf  l\*pIT(  JI*ANST-'2  o 3oO 

o  i  a  i  j )  =cF  nr-CPi\»cFi 4  83vo 

0I*3(  Jl  *CF  1<-CMN*CF  in  8400 

DIC(J)  =*CF  U~rclN»CFlC  !  0410 

/oir'(j)=csrl_-rrlN*<"erlr>  8420 

493  CONTI  V.Jr  8430 

IFI.NOIAGI  494.495.494  0440 

C  MAGNETIC  it  3430 

494  WRIT"  16.140)  8460 

•  RITE.  (6.104l(0IA(J).0In(J).JICIjl*ClD(J).J*l.NS)*  CFlA.CFlb.C  o470 

1F1C.CF ID.CFl" ,CC 1<*CF1V.CF1N»CF 1 .CF2.TF3.CF4.  CF5.CF6.CF7  o4o0 

2«CF8»(CF>!(I»1)*::NT(1), 1  =  1.81  8490 

495  GO  TO  401  8500 

C  GYROSCOPIC  ->OVfNT  MOTION  a51C 

500  CF2AsAa$(CFV(l .1 J-CFl )*ABSI'>  <12.1 I-CF2 )»A')S(CFV(3.1 I-CF3 ) ♦  AbS  8520 

llCFMI4.ll  -CF4)*'.  3S(CFM  (5.1  )-CF5)  ♦AdStCFVt  4,1  )-CF6)+  Ad.GICF  8  530 

2^( 7,1 I-CF7) .ATS(CFM( S, 1 |-CF3 t  8340 

CF2H*0.0  8530 

DO  501  1*1.8  i  8560 

501  CF2i3  =  CF23*Art.6(CF'1(  I  ,1  )  )  8570 

IFICF23)  502.603.502  8380 

502  CF2A=CF2A/CF23  d590 

303  WRITE  (6. 152 )M TC.CF2A  o6b0 

IFICF2A-DGYRI  506.506.604  8610 

604  NITC*N1TC*1  8620 

t  F ( Nl TC-N  I  T  )  505.506.506  8630 

506  CFl*crM( 1.1)  8540 

Cc7=CFM(2,i  )  8630 

Cr3 *CFV ( 3 , 1  )  8660 

C*4sCFV(4,l)  88/0 

CC5*CFV(5,’|  o6e0 

C  r6 * CFM ( 6  » l  I  8690 

C-7«CFV|7,1)  8700 

CF8=CFV(d.l)  8710 

GO  TO  481  8720 

606  WRITF  (6.1841  8730 

<MD  =  3  8740 

<NT(10)*1.0  8730 

GO  TO  456  8760 

C  ADVANCE  GPCFn  8  770 

613  IFINCAL-II  612.603,612  8 7a0 

512  JCAU*JCA(  *1  8  7Y0 

I  F  (  NC  A  t.  -  JC  AL  I  608.901.901 

901  111=2 

GO  TO  200 

50T  S3CAl *SPC AL ♦SP INC  8610 

AN6P*C.19471976*SnCAC  6820 

IFISoFN-jPC/.L)  60H.902  .902 

902  111*3 

GO  TO  200 

608  IF(INPOT)  609,201*309  6630 

C  9R0.'.7AY  EM*  9  6640 

609  WRITE  (6.9991 

999  FORMAT  MHl.;7MLA6T  SE  T  OF  INPl'TI 
STOP 


C 


XSI^r-jK  ;  1  AGNOGT  1 C 


6660 


/ 


s  1  0  Wfil  Tf  (6*  166  »mat 

go  n  m) 

ill  WRITE  i6.li4)».\r 
00  n  mi 

100  FORi'ATl  7?mo 

1  I 

101  forkati 72« 

HEARING  SUPPOK,'  PNOoilANC£  *LSP°HU  UF  *0T°k  *«T 
107  F)RV4T(120H  11  » 

1*.  T FFhNOLOGV ,  |  NC.  MECHANIC 

to*  f  JRva T I 1P*E15.7|  1 

105  f ORVAT  <10151 

106  FJRVATt l5.lPE25.6t 

107  FQRVA  r  (  '•oHOYCtittGS  ^COIJLUF  .SCALe  FACTOR  I 


loo  for, -41(120.  ;  NO.BRuS.  :.o,. 

1LEX«  3RG.M0v£nT  CvRO.KOM.  no.cases 

109  F0RM4 T ( iP5f 1 *, 7 1  *c 

110  F OR-AT (17.0112) 

111  FORVATOftHO  RQja 

1)2  F ORMAT ( 1P1E14,6) 

111  F  OR./A  T  1 1P5F  1  *fc6  ) 

11*  FORVATII7.1P1E20.7I 

115  FORMAT! I 7.1P5F20. 71 

116  Format  I  72rt  STATrot,  NO.  vASS 

1  SECT. INERTIA  > 

117  FORma T 1 1 2 OH  STATION  NO.  MASS 

15  *>lCT.  INERTIA  POLAR  MOV,  INERTIA  TRANSV 

118  FODMUMHIM  ’  ikansv 

119  FORMATn6M01TrRAT.  I  Tf  PAT.CONvERG.UMl  T 

120  FORMA’I  ISHOOEARIR'-  WATI-t'S  t 

121  fORMA  T ( 71  HQ 

1AT  STATION  NO,  |2| 

122  FORMAT  I  12  0*4  XXX  CXX 

1  <VV  cvr 

12J  FORMAT  I  1 36E 1 5,6  I 

12*  format i j 20R  Mxx  nxx 

1  MW  ->ry 

125  FORMA TI1P6E12.AI 

126  F0RmaT(I7,1?6ci6.a, 

127  Forma  Til  02H  bRi,  STATION  VASS.X-i;  IR, 

l  mass.v-oir.  xr 

128  FORMAT  I  78^0  ‘ 

1 TRANSl ATORy  motion  » 

129  FORMA  T  |  7HHJ) 

l ROTATIONAL  VOTIO.J  I 

150  Format  I 10XH  8RG, STATION  INERTIa. X 

1  INERTIA. Y  MV 

151  FORMAT  I l5.1P2F15*7t 
112  F  IRMA T 1 1  7 • 1P2E21 . 7t 

III  F  JRMATI5AM0UNML*NLv  ST,  *-UNBALAnCE 
11*  FORMATUSmOCOuPlINO  STATICN51 
115  FORMAT!  !8i<0OI  AGNOST  lc  1  | 

•I*  FORMAT  1 18*30 1 AGNQST 1C  ?  I 

117  FORV>|iA-lOO|AGNO.Srir  v  i 

118  FORMATIl«MOOI**NOS»:c  *  I 

119  F0RMA1 I 18H00I AGNOST Jc  •  » 

1*0  FORMATI18HOOIAGNOSTK  6  I 

1*1  FORmaTI*2moin|T|Al  SPEED  FINAL  S»EED  ' 
1**  F0RMATUP8E15.7) 

1*5  FORMAT! l^EEll.Ti 

1*«  FORMAT  1 11HQR0T0R  S»E F p. f )Pfc | *,7.}hRPM ) 


ROTOR  DATA  I 


length 

LENGTH 

TRAt.SV.KOM,  INtPT  I A 


CROSS 


HEARING 


<xx 

CYY 

CXX 

X  YX 

xxv 

cvx 

CXY 

) 

MXX 

ovv 

OXX 

MVX 

mxv 

0  YX 

OXY 

1 

INERT  I  a, x 

MV 


•“UNBALANCE 


REOESTAL  OATA. 

peoestal  data. 


^-unbalance 


FINAL  S»EEO  rnUF>  I  NCR.  I 


8tJ9U 

89uO 

0910 

8920 

8930 

89*0 

8950 

8960 

8970 

o9oO 

89  vO 

9000 

9010 

9020 

9030 

90*0 

9050 

9060 

9070 

9060 

9090 

9100 

9110 

9120 

9130 

91*0 

9160 

9160 

9170 

9180 

9190 

9200 

9210 

9220 

9230 

92*0 

9260 

9260 

92  70 

9280 

9290 

9300 

9310 

9320 

9330 

93*0 

9330 

9360 

9370 

9380 

9390 

9*00 

9*10 

9*20 

9*30 

9**0 

9*30 

9*60 

9*  70 

9*«c 

9*9u 

9500 

9510 


1W  OTITMOm  rflTHO'JT  GYPOSCC/PIC  •'J'-e \T  I  98x0 

laa  cr>o/.\r  |  I4.1?E18.6.1P3C14.  .*1PE1  ^■3i0 

u>  c  -jsvA  t  ( : ->c  n,* •  i pieu.e  i  954c 

1<<i  ro«?'-'iT  I  ]  10*  ST4TIC\  .va  JCP  AXIS  MINOR  AXIS  ANGLE  X-MAJOH  PHASE  9850 

’  A\:,Lr  "AJO?  AXIS  1  NOR  AXIS  A \GL t  X-lVAJOR  PHASr  Ai>CLt  *  VSoG 

it-  Fnpvjtnrj  amplitude  9570 

:  RE NOlNG  VC«tNTI  95  =  0 

FORMAT ( I 7,’PE 1 7. 7)  5890 

laj  f okma r ( 24mo 1 7 fra r •mg*  i  *6pu 

154  E  ORwA  T  (  24H  .VlTrt  G7K0SG0PIC  VOMLNT  )  5610 

ltt  format  <  30H0GVFP/UYH  RFLO*  1%  XSU'ECF  AT  III  9620 

155  F JRMA T (  -X.HOVATP I  X  IS  SINGULAR  IN  XSIMfcGF  AT  III  *6 iJ 

1 S  7  FORMAT! 1PAE14. 6)  9640 

710  F0RVATI16H0  F~9CE  T R A NS* I T T t J  TO  ri*  AR  I  N$  HOUSING)  96:>0 

711  F0pvAtn3H0  FC<CS  TRANSVITTF''  70  FOUNDATION)  •  7660 

7  ■»?  FQRVA  T  (  j  8H0  PFOEFTAL  MOTION  I  9670 

FCRVATI120H  HPG.NO.  V AJOR  AXIS  MINOR  rfx I S  ANGLE  X-MAJOR  PhASE  96o0 

I  INGLE  X-AV°LITJDE  X-PHASt  ANG  Y-AMPLITUDE  Y-PHASt  ANGI  96*0 

7 14  FOPVATIISHO  ENERGY  INPUT., IP L  14.7.21 H  ENERGY  DI SSIPATEC*  » 1  Pc.  1 4.  7  9700 

II  9710 

END  9720 

tIHrTC  ‘.UdSS  Y94.xR7.Ll  iT.NOOEOC 
SUBROUTINE  SURR 

0  I  ME  NS  ION  JMXfllO.  ROI *BMXS( 10.  80>.RMVC<10»  BO  I  .OVYSI 10.  BOI.VXC  o030 

1 <  lO*4l > .VX5( 10*41 ).. 7C ( 10.41 1 .VY5( 10.41 I .DXCI 10.40) .0X5(10.401 «OYC  0040 
21 10.401  .On  (  10.401  .XC(  IC.40  I  »XS(  10.40 1.YC1 10.40)  .YS  (*10.40)  .0MXM40  0050 

1 )  .PVXRI40)  .PYXC  (40)  OMXD(40>  .0'* Y  A  (  4 0 )  .D»Y1  (  40 )  . DMYC I  40 )  .DM  YD  ( 40  )  .  OJ60 

4PVXA(40)  .PVXRI4C)  .OVXC  ( 4C  )  •  '  VX^  (  40 )  .PVY  A  (  40  )  »DV  YB  ( 40  I  .OWC  ( 40 )  •  0070 

5DVYD(40> .RM(40) .RL!4QI »R S( 40  I »R IP ( 40 1  •« ! T ( 40 1 .01 A 1 40 ) »Dl B ( 40  I .0  I C (  0080 

640) .010(40) .AN (40) »oN< 40 ) .ON <40  .UX(4CI .07(40) .LUI40) .Li*  (26) ,dKXX(  u090 
7  3.25) .8CX xi 3.25 l.b<xY( 1,251 .oCXY (?. 25). RKYY I  3.25  1. &LYY  I  3,25). dXYXI  uluu 
81.2S),UCYX(3.25).RsVXX(3.25I  .30MXX ( 3 .25  I  .BaMXY ( 3 .25 ) .BDMXY ( 3 .25 ) ,  ollU 

9aSMYY  (3,25)  .d-'MYYI  3,25  )  , dSMYX  (3.251  .HPMYX  (  3.25  )  0120 

01 Mt NS l ON  pyx  I  2 5  I ,PVY ( 26 ) .P<x(25 ) .PCX ( 25 ) »PK Y ( It  ) .PC Y ( 2 5 ) »P I  X ( 2 5 ) .  «13u 

IP  I Y (25  I .PSMXI25) *POMX( 25) .PbMYl 25  I . PCYY (25) «LC( 20) ./*( 8.8) »d( 8.4) ,  0140 

2C(8.4 I .p( 8.1 ) »F (4 .4.20) »CFM( 8,8  I *EN7 ( 10) . RhS ( 8  *  1 ) ,IK'UX(80) .  ,  ul5U 

3  OVUYiaOl  .DUMMY  (  300)  ,3HA(  3  I  ,  SHd  (  3  I  »SHC  (  3  )  .SHL(  }  )  .AUXI29.4.3  )  .  J16U 

48UX  (25,4.3  )  .P’lVl  (  ICO)  *CLNR(  8  )  01  /C 

COMVON  A  ,  a  ,  C  .0  ,  CF*  •  cNT  0200 

COMVON  PHS  ,  Ou**l  ,  Cc 9  ,  ME T  ,  XFN  •  CLNR  0210 

COMMON  PRN3  0220 

COWON  NS.Nd.NJ.NC.NPST.NMOM.NGYR.NCAL.NPIAG.INPUT.YM.SCF.RM.RL.  0230 

IRS, NIT, OG  YP.PjP.r; 7 .LR.LU.UX >UY  »LC  »PMX ,P<X .PCX  *PMY .PKY.PCY.XST.  0240 

2PIX.PSYX.P0I'X,PI r.PSYY .PpVY .SPST.SPFN.SPINC.dXXX.HCXX.ttKXY.BCXY.  0280 

38K  Y  Y  ,  RC  Y  Y  .HSYX.RCYX  »)3S  VXX  .KOM  XX  .  d.SMX  Y  .  3nMXr  .dS.M  Y  Y  .dOMYY  .BSMYX  » 
4dOMYX.CFU.CF  1C.  OF  10, CF  1E,CF2<«CF2C  ,CF  2P.CF2t  .CF2A.Lf2d.CF2M.LF2N, 
5CFlA,CFld.CF  l.V,  CF  1  N. CF  1  ,CF  2, CF  3. CF  4. CFS.Cf  6  »CF  7. CF8.STF.AN5P2. 
6ANSP.SPCAL,OVXA.OVXd.OVXC.Dyxt)..OVYA.nvYR.OVYC»PVrO,OMXA»OMXB*DMXC»  028 


7DMXD.0MYA  .Pi'  YR  .PM YC  .pMYD,  AN  »‘JN  .PNi  I  I  I 


COMMON  dMXC.l3yXS.9MYC.BMYS.VXC.VXS,  YYC.yrSfOXC.CXS.OYC.OVS.XCJ.XS. 
lYC.YS.OIA.OId.OIC.  7  I  n,  PVDX  .OVIJY  .OL’MMV  .SMA.SHd.SHC  .SHO.AO’X  ,BU< 
COMMON  JCAL.AyS.NITC  j 


100  FORMAT! 72HO 

/ 

8930 

1  i 

i 

8940 

101  FORMAT  I 72h 

1 

8980 

1  1 

1 

8960 

102  FORMAT (120H1 

UNBALANCE  RtSPONSt  OF 

RqTOK  »it 

ov /O 

1M  NON-UNIFOKV  HEARING  SUPPORTS 

PN001 1 

/  * 

8980 

103  F0RMATI120H 

mechanic 

8990 

1A*.  TECHNOLOGY, INC. 

/ '  ». 

9000 

1 04  F  )RMAT< 1P4E15.7I 

/ 

VvlU 

105  F'IRMA T  1 101  5  1 

9020 

106  F JRMATII5.1PE23.6) 

9030 

1q7  FORMAT OOHOVOuNGS  MODULUS 
108  FORMAT! J20H0  STATIONS 


scale  factor: 

NO.URGS,  N0,,UN3Al. 


NO.COUPL.  PED.F 


HEX*  8R0. MOMENT 

109  FORMAT! 1P5E15, 71 

110  FORMAT! t 7,91 121 


GYRO, MOM,  NO, CASES 


DIAGNOSTIC 


INPUT 


) 


111  FORMAT ( 35^3 

112  FORMAT ( 1°3E1*,6I 

113  FOR**\T (  1P5£1*,6! 

11*  FORMAT( I7.1P3E20.7I 

ROTOR 

DATA  I 

115  FORMAT IIT,1P5F20«7I 

116  FORMATI72H  STATION  NO, 

MASS 

length 

cross 

1  SECT. INERTIA  1 

117  FORMAT! 120H  STATION  NO, 

MASS 

LENGTH 

CROS 

IS  SECT. INERTIA  POLAR  MO* 

,INER(|A  TRANSV 

.MOM, INERTIA 

1 

118  FORMAT (1*151 

119  FORMA H36H0ITERAT,  I TERAT ,CONv£RG,L 1  MI T 

120  FORMAT (18HQ9EARING  STATIONS  ( 

1 

121  FORMAT ( 71H0 

BEAR  I  ' 

1 AT  STATION  NO.  I 2 1 

122  FORMAT! 120M  XXX 

CXX 

XX  Y 

Ci.  f 

1  XYV 

CYY 

KYX 

CYX 

) 

123  FORMAT (1P8E15.6) 

12*  FORMAT! 120M  M*x 

DXX 

MXV 

DXY 

1  MTV 

OYY 

MVX 

DYX 

1 

125  FORMAT <1P6E 12.61 

126  FORMAT  1 I T . lP6r16,*l 

127  FORMAT I102M  faRG, STATION 

MASS.X-OIR, 

xx 

cx 

1  MASS.Y-OIR, 

XV 

CY 

1 

128  FORMAT! 78H0 

PEDESTAL 

DATA. 

1  TRANSLATOR Y  MOTION  1 

129  FORMAT  I 78*0 

PEDESTAL 

DATA, 

1R0TATIGNAL  MOTION  1 

130  FORMAT ( 102H  URG. STATION 

INERTIA, X 

MX 

DX 

1  INERT  I  A, V 

MY 

or 

» 

131  FORMAT! 15. 1P2F 15,71 

132  F JRMATII7.1P2F21.7I 

133  F  JRMAT 1 56M0UNRAL ANLc  ST, 

x-unbalance 

V -unbalance 

1 

13*  FORMAT  I18H0C0«JPLING  STATIONS! 

135  FORMAT! 18HODIAGN05TIC  1 

1 

1*1  FORMAT(*2HO!NIT(AL  SPEED 

final  speed 

SPEED 

I NCR.  1 

1**  FORMAT (1P8E15.7I 
1*5  FORMAT I1P6E15.7I 

1*6  FORMAT! 13HOROTOR  SPEED* » 1P£ l*»  7.AHRPM I 
156  FORMAT! ]*HOMATRIX  IS  SINGULAR  IN  XSiMEQF  AT  III 
157  FORMAT  ( 1 P*E 1*,6 I 

GO  TO  1900,227*3001*111 
900  READ  ( 5, ICO  I 
REAO  15*1011 

READ  I  5* 105 INS*NR*NU*NC , NP£T ,NMOM.NGYK,NCAL *NOI AG, INPUT 

READ  15* 10* I YM,SCF 

WRITE  (6*1021 

WRITE  16*1031 

WRITE  (6*1001 

WRITE  (6*1011 

WRITE  (6*1081 

WRITE  l6*110>NS*Ne,NU*NC*NJ>S?«NMOM,NGVM»NCAL*M3!AG, INPUT 
WRITE  (6*1071 
WRITE  !6,I0*IVM,SCF 
IF(NGYR)  202*201*202 

201  REAO  (5*112l(RM(jl,RL|Jl,RSIJ},J*l*NSI 
WRITE  (6*1111 
WRITE  16*1161 

WRITE  (6«ll*l(J,«N(Jt,RL( Jl*RS( JI*J*l»NSI 


90*0 
9050 
9060 
9070 
9U«0 
9090 
91 UO 
9110 
9120 
9130 
91*0 
9150 
9160 
91  70 
9130 
9140 
9200 
??10 
9220 
9230 
02*0 
9230 
9260 
9270 
9280 
9290 
9300 
9310 
9320 
9330 
93*0 
9350 
9360 
9370 

mD- 

9390 

9*00 

9*10 

9*20 

9**0 

9*90 

9500 

9510 

9630 


0*00 
0*10 
0*20 
0*30 
0**0 
0*50 
0*60 
0*70 
0*80 
0*90 
0500 
" "0310 
0520 
0530 
o»*0 
0550 


>0  '0  -01  ./3ti() 

2V  -ha:)  t  b  •  I  c^i  » V  I  T  .  u5/U 

A  <!Tf  16.1 I  V  I  •  u5oj 

WRITE  0599 

WRITE  IF.  Ill)  U6u0 

READ  I 6  , 11 3'  f  R  v  (  J  I  »RL(J)»W5(J) »« I P ( J  •  »  B  I TIJI *  J  *  1 , NS  I  J6H 

A { T  r C  IF,  1 1  » i  062C 

%  <1  rr  I  i.  1  17  I  (  J.  i'M  jl  .RLf  J)  .R5I  J)  ,4  IP  I  jI.HITI  J)  ,J*1  .NS)  06  if. 

2M  RIAD  <5.1191  IUIIJI.  1*1  *  K’  =3  »  0t>60 

w<ITE  I  6 ♦  12  01  0650 

WRITE  (6.U9)tL3IJ).J*l.N3)  0660 

212  1cln  15.  m  |  (LtJ(  Jl  .0X1  J)  .'JY(  Jl  »J*l  ,NUI  0670 

a  W | T  E  (6.1^11  >  o690 

WRITE  16. 1 32 1  (UK  J»  ,'JX(  J>  .JVC  Jl  »J*1  .NU)  •  0650 

216  Iff  NC  I  215.207.215  0700 

216  READ  (5.110) (LCtJ».J«l. NCI  *  0710 

WRITS  (6.116)  o  720 

ViRItr  (6.119|(LC(J).J*1.NCI  ,  .  ‘  u7  JO 

?r  7  ici'.ncri  20a.229.209  0740 

v  R  t  T  E  I  a. 1 ?9 I  0750 

WRITE  (6,1271  0760 

DO  2C9  J* 1 .NU  0770 

Rt  AD  (6.125)  3M<  ( J I  *0<X  ( J)  .PCX ( J  )  »PT‘Y  ( J  I  «  PK  Y  ( J )  .PC  Y  ( J I  07bO 

K5T=!_S(J)  u  790 

209  WR 1  T  E  (6.  126  KST.P'fxi  JI.PXXl  Jl  ,PCX(  Jl  ,P»Y(  Jl  »P<Y(  J)  *  Pf.Yl.ll  vpOO 

I f (  lvOv 1  210.229.210  oalO 

210  WRITE  (6.125)  '  *  0020 

,  .RITE  (6.1301  ‘  oOJO 

DO  211  J* 1 ».NH  0040 

READ  (  5. 1 25 IPlX ( Jl .PSw* ( j) ,P2FX( J) .  =  1 Yl J) »P5MY| J) »  PDMVCJ)  0050 

<ST*l_R(JI  0060 

211  WRITE  (6.!26lxe.T.PIX(J).P^,*Y(J).P0MX(J)tPlY(J),oSMY(J),PDWY(JI  0870 

228  IFC.CAl-ll  27l.220.27t  OOOO 

220  READ  ( 5 , 1 12  I rP6 T . SPFN.FP1NC  0090 

WRITE  (5.1 41)  '  0900 

W R  I  T c  (  6. 1  12  I  r-P5T  .6PF\»5PJ  \C  .  0910 

READ  16, 1 12  » (  M<XXI  ’  .j)  ,  !*1 .3)  .(.)CXX(  I  ,JI  .1*1,31  ,  (oFXYll.J  o920 

11.1  =  1.1).  (•)(«’'(  I  .Jl  ,  1*1.31  .(2KYY  I J  ,j)  ,  I  *1 ,3)  ,  ItiCYYl  |  ,«,)  ,  1*1.31.13  -9  JO 

2KYXI  I  .  J),  I*  I  ,7)  ,  l:.*CvX<,l  .J)  .  1*1.1)  .J«l  ,N«|  0940 

•;0  204  J»  1  ,\R  .  3950 

<ST*Ll(JI  0960 

WRITE  (6.121  I Y  ST  0970 

WRITE  (6.1221  09o0 

2r 4  WRITT  ! F , 1 77  I ( R< AX ( ) . J ) Tx ( I . J | ,fKXY ( | » J | ,hCX Y ( I » J ) .  PKYYM.Jl.f)  0990 

1CWY( I « J I »  ?<  f X ( I , J ) »  JCYX  t I ,JI .1*1 ,3)  -  .  10C0 

I  ~  ( NMC'I )  205.216.205  1010 

205  READ  (5,1 12) ( (3SYXXI I, J). 1*1* 7) , <  Hl'YX  X ( l.J). 1*1.3 ),( G5MX  Y ( I .  1020 

1  J)  .  1*1  .3)  .  (dD«XY(  I.  II. 1*1, 11.1  WIYI  I  *J)  .  1*1 .31  «I4DMVY(  I  ,J  1.1*1  1030 

2,JI.(8S"7VX(I,J),|*1,1),( BDMYX (J,J)«1*1»?),J*1 ,Nh )  1040 

DO  206  JM.NB  1060 

<3T«U0(JI  1060 

WRITE  ( 6.  12 1  I x  6T  1070 

WRITE  (6,124)  lwo0 

206  WRITE  (6.173  I  (HSwxxi  1,J),RDMXX(  |  .Jl  ,HSMXY(  I  .Jl  .iiiJWXYU.JI  .BbRYVII  ,  1090 

tJI.BIRYVf  I.JI.HS'AYXf  !.J),»D»YX(  |.J|, |*1,1)  •  *  HOC 

f.O  TO  216  1110 

221  JCAL-1  '  1120 

22  T  READ  15.167)  OST  1130 

RE AO  ( 6. 157  I ( 0<XX( 1  ,J) .RCXX ( 1 , J) ,H<XY ( 1 ,J| .dCX'  l ,J I .  B(YYIl.J)  1140 

1  ,aCYY(l*J)  .J*YX(  l.J  I  .ncvx(  l.J)  »J»l.Ntl)  1150 

WRITE  C6.14AIXPST  1160 

00  221  J*l.NH  ....  1J70 

<  ST*l 8( J I  1180 


I6,122> 

IMNMOVI  22*. 256. 22* 

. - 

*ST»L8( Jl 
*»ITC  (6.121 I  xST 
«IT5  (6.12*1 

00  22 6  1*2.) 
noon  i.ji. o.o 
acxxti.ji.o.o 
Btxvf |,J|«0.0 
8'xrri.ji.o.o 

B'.YYI  I.JMO.S 

acrvi i.ji.Q.o 

a;yx(i,j|.o.o 

8CVX(  i.jlo.o 
asMxx(i,j|.0.o 
aoMxxn.jj.o.o 
85«xy(i,j|.0.0 

aoMxvn.jj.o.fl 

asMYvn.ji.o.o 

aoMrvn.ji.o.n 

asMYxd.jt.o.o 

BDMYX ( 1 . J 1 *0,0 
226  CONTINUE 
216  GO  TO  25 0 

250  ANS., 000.0  *0mi€*'  ,M,Ut  U*m 

CFl*366.06V«Aw$ 

CF2-CF1/2.0 
IIS  ( NS  I *NSI 1 1 
00  251  J*1,NS 
RM( Ji.ftMIJi/cri 
B1P( Jf.NlPf  H/CE2 

B1TU|.R1T|J,/Cri 
STr.rN/AMS*RS(JI 
*N(J|.Ru ji/stf 
*"‘j»*«uj>/2.o*aniji 

251  0N( J| *Ri( J ) /3(o*UNI j I 
00  252  J>l,NU 

...  U*»J7"U*(J1/*I77,1 

252  UY| J|»UYI JI/61 /7#l 
SPCAU-SPST 

ANSP*0.1067I«76«SPCAl 

NITC-1 

300  ANSP2.ANSR.ANSPS,>eE0 
00  101  J. 1 ,NS 

ONXAl J)«o,0 
OMxSUl.o.o 
ONXC(J»«0.0 

0«XO|J)>0.0  ,  . 

ONyaiji.o.o 

04Y8I JI-O.O 

ONYCUco.O 

0«YO(J»«o.O 

SCF.RWt U»*6NSP2 


li  vo 
1200 
1210 
1220 
1230 
12*0 
1250 
1260 
1270 
12B0 
1290 
1300 
1310 
1320 
1330 
13*0 
1350 
1360 
1370 
1360 
1390 
1*C0 
1*10 
1*20 
1*30 
1**0 
1*50 
1*60 
1*70 
1*60 
1*90 
1500 
1510 
1520 
1530 
1 5*0 
1550 
1560 
1570 
1560 
1590 
1600 
1610 
1620 
1630 
16*0 
1650 
1660 
1670 
1660 
1690 
1700 
1710 
1720 
1730 
17*0 
1730 
1760 
1770 
1760 
1790 
1600 
1610 


cvxai  j)*srf 
DVYA ( J  |  *S  TF 

f)VX3<JI*0.0  •  f  ... 

0VXC(JI»0#0 

DVXCHJI-0.0 

ovYaui»o.o 

0VYC<JI«0.0 

dvyoi Jl»0.0 

DIA«J)*0.0 
DIBIJI'O.O 
DIC(JI«0.0 
3M  D 1 0 ( J I *0*0 
3n?  00  111  J«l,N8 
«ST»L8< Jt 
XPN»0 

CF1K*B*XX< l»J)+dXXX(2*J)  «ANiP*UXXX( 3  t J ) *ANJ»P2 
CF1C«  dCXXI  l»J)+aCXX(2»J)*AN.SP+oCXX(3tJ)»ANiP? 

CF  1D»BXXY ( 1 •JI*8XXY(2t.'l*ANSP+6XXY( 3 t J)«ANSP2 
CF  IE*  dCXYl  1»j)-*-bCXY(2  • Jl aABSPybCXY  I  3  •  J  I  *  ANSP2 
CF2X«8XYY(  1  .J|A8XYYI2»  J»*ANi.P>dfcYYI  3,J)»A.NSP2 
CF2C*  BC Y  1  ( 1  * Jl ♦BCY Y( 2  * Jl *ANJP*ot YY I  3. JI»AN5P2 
CF20*8XYX  11  »J)+0XYX(2.J)  •Ai'«SPa8xYX  (  3  *  J  I  *ANSP2 
CF 2t*  8CYX1  l.j)*dCYX<2,j|»A»tSP*dCYXI3.J|«ANiP2 
IF(^P$TI  303*336*303 

303  CF1 *»Pxx( Jl -Pyx ( J 1/386 *069*ANSP2 
CF  1N-PCXI J)»ANSP 
CF1A«CFU*CF1m 
CF1B«CF1C»CF1n 

Cc2M«PXY 1 J I -Pmy ( j )/386. 06 9*ANSP2 

CC2N*PCY<  j ) *ansp 

C“?A«CF2<*CF2m 

CV2b*CF2C+CF2N 

03  TO  307 

304  X-'N«1 

C  F 1 X  aBSMX  X ( 1 , JI+BSMXXI 2* J)*AN*>PfBSMXX(3.J)«ANSP2 
Cf  1C*  BDMXX ( 1 • J ) tSOMXX ( 2  • J I • ANSP  +  BOMXX ( 3. Jl «AHSP2 
CFlD*«SMXYU.JI»B3MXY(2*JI*ANSP4dSMXy(3.JI*ANSP2 
CF1E-  BOMXYt 1  ,J»+80MXY (2*J)*ANSP*H0MXYI 3. J)»ANSP2 
CF2X-BSMV Y ( 1 , J | ♦daMY Y  < 2  .  J I  »ANjP*dSMY Y I  3  »J I *ANSP i 
CF2C«  80MY Y ( 1  * J I *B0MYY ( 2  * J I *AMSP>dOMYY (  3*  J) *ANSP2 
CF20*B 5MY  X  I l*J)+d$MYX< 2. J )  •A,NSP-*ttSMYX  I 3*J)*AN£P2 
CF2E-  BOM  YX  (  1  •  J  I  ♦BQMYX  (  2  • Jl •ANSP+8DMYX  I  3*  J I  6ANSP2 
IF(NPST)  306*306  •  306  ,  . 

305  CF1M«PSMX ( JI-PlXI Jl  /  386.069»A4SP2 
CF1N«PDMXI JI»ANSP 
CF1A*CF1K«F1» 

CFld-CFlC^CFlN 

CF2M-PSMY ( JI-PIY ( J 1/386. 069* ANSP2 
CF2M«P0MY( J|**NS° 

CF2A-CF2X*Cr2M 
CF2B-CF2CACF2N 
GO  TO  307 

306  CF1-CF1X 
CF2.CFIC 
CF3«CF10 
CF4*CF If 
CF3-CF20 
CF6-CF2E 
CF7-CF2X 
CF8-CF2C 
GO  TC  308 

30*  CF6«CF?A»CF2A*CF2B»CF2B 

CFl«lCF?A*CFir*CF2d»CF2t)/CF* 


18/0 
1630 
1840 
1650 
i860 
lo  70 
1600 
1690 
1900 
1910 
1920 
193* 
1940 
1930 
1960 
1970 
1960 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2060 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
21oO 
2190 
2200 
2210 
2220 
2230 
2240 
2230 
2260 
2270 
22*0 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2360 
239V 
2400 
2410 
2420 
2430 
2440 


a  /•‘(Cr2A*CF2t-CF/o*CF2L>l/CF4 
Cc  )« (C*-  2A»(.F  2v*CF  2.>*CF  2M/CF  A 
2A*Cr2N-CF/o»CF2MI/CF4 
C**'»CF  l*-tU»Cf  1  J*U  2«CF  U 
C-6-CF  1i1-CF2»CF  lU-«.r  1  »C*-  It 
C  :7«-CF  **CF  1  L»*CF  4*C  ^  It 
C*8—Cr4»CF  10-CF  l*t»  It 
CF2N.CFS.CF5.CF64CF6 
Cr2A«CCFS*CF  Im*CF6*CF 1 N I /CF 2* 

CF  2b*  I  CF  j*CF  iN-CF**CF  1MI/CF2N 
CF2v*ICF5*CF  /*CF6»4Fbl/CF/N 
CF2N*«CF5*CF0-CFS»4F  71/CF  2N 
CF  1A.-CF 1*CF2A*CF  2*CF2d 
CFlb«CFl*CF2.i*CF2*CF2A 
CF 1M-CF J-CF l«CF2M*wr2»CF2N 

CF  1N*CF4-CF1*CF2N-CF2«CF2M 

CF1.CFU*Cf2A-CF1C*CF2BfCFIu»U  1  A»CF  lt'CFlO 
CF  2*CFlX»CF2b*CF  1C*CF2A-CFI  J»v.f  1j*CF1E«CF1A 
CF1«CF1A»CF2M-CF1C»CF2MCFU*CF  IM-CF11*CF1N 
CF4*CF1A»v.F2N*CF1C«CF2M*CF  1  i«Cr  \N*Cf  1E*CF  1M 
CF5*CF20»CF2A-CF2E»tF2tl*Cr2A*CFlA*CF2t*CFib 
CF6*CF2U*CF2o»CF2£*CF2A-CF2A«CF lb*CF2C«CFlA 
CF  F«CF2D»CF2M-Cr2E«tF2N*CF2A»CFlF-Cr2C*CFlF» 

CF0*CF2D»Cr2N*CF2t*CF  2M*CF2<*CF 1N*CF2C#CT 1* 

100  IF  (AFN  I  110*309* i 10 
109  OVXAlASt  l»0V*AUiT  )-CF  1/AMS 
OVXBi AST  I »CF2/AM$ 

OVXC ( AST ) «CFi/AMS 

DVX0(ASTI*CF4/ANS 

DVYAIASTI *0VY4I<ST l-CF 7/ AMS 

DVYB|XST»*CF8/AMS 

l)VVCtASTI«CFS/AMS 

:JVY0|XST)*CF6/AMS 

1F(NM0M|  104  *111*  304 

no  umxa(ksti*cfi/ams 

DMX8C  AST  I *CF  2/AMS 
OMXC (AST  I «CF 1/AMS 
OMXOI AST  I *CF  4/AMS 
OMVAI AST  t  *CF //AMS 
OMVBf  ASTI *CF8/AMS 
DMYC«ASTI»CFS/AKS 
OMVQIASTI-CFS/AMS 

111  C1NT1NUE 

IF'INDIAQ)  112.311*112 
C  31 AGNOSTIC  1 

112  WHITE  1 6*  1  IS  I 

WHITE  I6*10»ICF1X.CF 1C.CF lu.CF  It »CF 2K .CF2C.CF 20 *CF2t •  CF2A«CF2B*C 
1F2M»CF2N*CF1A(CF1o«CF1M»CF 1N.CF 1 *CF2 *CF )*CF**Cf  a*CFft*  CF7*CFo*STF 
2*AN$P2*ANSP*SPC4L 
WHITE  1 6*  1 10 ) AST *AF N 

WHITE  <*•  144II0VAA1JI*  uVXtl  (JI.OyXCIjI  .OVXOI  J I  ft)  VTA  I J  I  •  DVViH  Jl  *OVY 
1CIJ1 (OVYDI Jl *J*1 «.>■»* *I3MXAI Jl *OMXb< Ji *OMXC( JI.IMAUIJI*  OMVAI J»*OMY 
2BI JI.OMYCI  Jl  .OMYDI  Jl.J-l.VSI 

WRITE  I6.145IIRWI  J»*R|P1  J|b«|T(  J)  *AM  u  >  *B»<  J)  *ONI  J I  •  J*1*NSI 

111  RETURN 
CNO 

SIBFTC  ARCT  M94.XR7 
C  ARCTAN  ROUTINE 

FUNCTION  ANO) AFCS.AFSNI 
ACSoAFCS 
AJN-AFSN 

•00  IFIASN1  00A.0ftl.OO4 
.  A01  IF (ACS)  802*«0l*«0t 


A4S0 

2460 

2470 

24o0 

24*0 

2S00 

2S10 

2S20 

2S30 

2S40 

2SS0 

2S60 

2S70 

2SoO 

2590 

2600 

261u 

2620 

2630 

2640 

2650 

2660 

2670 

260o 

2690 

2700 

2710 

2/20 

2730 

2740 

2750 

2760 

2770 

2700 

2790 

2000 

2010 

202U 

20Jv 

2040 

2050 

2060 

2070 

2000 

2o90 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2900 


M020 

0010 

0030 

0040 

u0>0 

0060 


802  A\G*180.0 
GO  TO  812 
ANG*0.0 
GO  TO  812 

804  (FfACSl  80V. 805,808 

80^  IF(ASN)  808 ,803.80 7  >' 

806  AN G*-90.0 
GO  TO  812 
802  ANG*90.0 
GO  TO  812 

808  ANG*0.0 
GO  TO  810 

809  ANG*-180.0 

810  A  SN*A.SN/ACS 
A'tS^ABSlASN) 

A:s«ATANJ ACS) 

MG»ANG4ACS»57.2V5  780 

f-^SN)  811,812.812 
oil  AHG«-ANG 
912  RETURN 
t-O 

41BFTC  EOSS  M94,*R7 

SUBROUTINE  ECJS 

1CUM1U00)  ,CLNR*8)  ,','8,4,*C'a»*>»Ol8.n.CFM(8.8>,fcNTI10),RHSt8,U, 

COMMON  A  f  u  r 

COMMON  RMS  ,  0UM,  ;  •  0  •  CFM  ,  £*T 

COMMON  PR.N3  *  •  *1AT  ,  KFN  ,  {.LNR 

<FN«*CFN 
KFN«1 

3*0  oS  ^4'J:?;ffl*«0.280».MAT 
PRNR*l.o 
PHN2-0.0 
PRN3*1.0 
PRN4«0.0 

A  4*0 

00  248  1*1, A 
PRN1.AU,  JJ 
1PIPRN1)  242,2 41,243 

241  PRN2*PRN2«1 .0 
GO  TO  244 

242  PRN1.-PRN1 

243  PRNR*PRNR»JPRfc,i,#0  ,5J 

2A4  PRNl.ail, J»  ' 

.  IPfPRN ll  246,243,247 
243  PRN4*PRN4U*0  • 

•C4.IC441 
GO  TO  248. - 

246  PRNI.-ornI 

247  PRN3.PRN3* (PRnI* *0.231 

248  CONTINUE  '  • 

249  ORN^H^.f50’250*^9 

2*»  ®W*A.0/|*,o.P)li,?,  •* 

P  INR*PrnR*#prn2 

230  C5.NRI  J)*PRSR 

23.  «2.252,231 

231  l*(<4>4|  237,252,2!,^ 

257  p'*NA«4,0/<4,0-PRN4) 

PRN3*PRN3*«PRN4 

232  0UM1(ji*prn3 
00  231  1*1,4 

All  »J)*A(I (JI/PrnR 


U070 

0080 

0090 

0100 

0110 

0120 

ol30 

0140 

0120 

0160 

0170 

0180 

0190 

0200 

0210 

0220 

0230 

0240 

0230 

0260 

0270 

uoio 

uOZo 
0030 
0060 
00  7u 
0080 
0120 
0130 
0140 
0150 
0160 
ol/u 
0180 
0190 
0200 
0210 
0220 
0230 
0240 
0250 
0260 
m2  70 
o2oO 
0290 
0300 
0310 
0320 
0330 
0340 
0330 
03*0 
0370 
03«0 
0390 
0*00 
o*10 
o*20 

0*30 

0**0 

0*50 


y  ,»t-  i'U 


t 


254  CONTINUE 

CALL  «ATINVI*,4.a.4,CR9» 

00  256  J«l,4 

PRN  t -DUM1 1  J  > 

00  255  1-1,4 

ill  t i 1  *Jl/CLN«m#PRNl 
256  CONTINUE 

00  TO  29/ 

260  PRN3-1.0 
PRN4-0.0 
K4-0 

00  271  J-1,4 
PRNR- 1 ,0 
PRN2-0.0 
00  26A  1-1,4 
PRNl-Cl I *J) 

*«,PWn  2*2  *261  *263 

261  PRN2-PRN2+ 1 «0 

00  TO  264 

262  PRN1--0RN1 

263  PRNR-PRNR- ( PRN1 -*0-2S t 

264  CONTINUE 

n,  i^RN2>  266,266.265 

265  PRN2-4,0/(4cg.pD«|2| 
P‘TNR-PRNR#-P«N2 

2 66  C'.NRl  JI-PRNR 
0)  267  1-1,4 

267  C ( I ,J I -C( 1  * J) /PRNH 
PTN1-01J,1| 

1-1PRN1I  269,266,2/0 

266  PRN4-PRN4,} ,q 
K—K4,l 
GO  TO  271 

269  PRN1--PRN1 

270  PRN3-PRN3 > I PRN1 **0-291 

271  CONTINUE 

273,273,272 
222  1FTK4-4I  277*271,273 
227  ‘*£n*-4.0/<4,0'*PRN4I 

PRN3-PRN3—PRN4 

273  00  274  J-1,4 

274  01  J,1 |-0| j#1 )/pRN3 

274  rA,jA>!5l,*l,/CUW,*»#w‘»»3 

«u  TO  297 
280  PRNl-1,0 
PRN4-0#w 

K4-0 

00  291  J-1,6 

PRNR-1,0 
PRN2-0.0 
00  264  1-1, « 

PRN1-CFN1 |,j| 

IF  IPRN1I  262,261,263 
261  PRN2-PRN2-1.0 

GO  TO  264 
282  PRNl—PRNl 
261  PRNR-PRNR-(PR,v|,,Qa  123  | 

264  CONTINUE 

J'WI  266,264,265 

265  PRN2-6-0/ 1 4«0»PRN2 ) 


04/0 
0460 
04  90 
0500 
U510 
052C 
0530 
0540 
05*0 
0560 
0570 
0560 
0590 
0600 
0610 
0620 
0630 
0640 
06*0 
0660 
0670 
0680 
0690 
0700 
0710 
0720 
0  730 
074Q 
0750 
0760 
0770 
0760 
0790 
0600 
081C 
0820 
0830 
0840 
0850 
0660 
0870 
I860 
0890 
0900 
0910 
0920 
0930 
0940 
0950 
0960 
0970 
0960 
0990 
1000 
1010 
1020 
1030 
1040 
1050 
1040 
1070 
1060 
1090 


p <\K-PrtN*»»DHN2 

286  CLNRI JI*PRNR 
0)  28?  1*1*8 

28  7  C-Mll  .Jl-CFMI  l.j)/).*  W 
PRN1«Rh$(  J*1  ) 

Ir  (PRN1I  289,288.290 

288  PRN4«PR*'<.*1  ,o 
<  9  *K4>  1 

GO  TO  291 

289  PRNl »-PRN 1 

290  PRN3-PRN3»(PRn1«*0.125) 

291  CONTI NU£ 

IF  (PRN4I  293,293,292 

292  IFHC4-8)  298,293.2*3 
298  PRN4-8.0/(8.0-PRN4l 

PRN3-PRN3—PRN4 

293  00  294  J* 1 , 8 

294  RMSI J,1 l-KHSI J,1 I/PRN3 

C*U  MATlNV(CFM»8»8HS»l»CF9) 

295  DO  296  1*1.8 

*96  CFM< I .1 ) *RHSII»1 J/CLNRl I »*PRN3 
297  RETURN 
END 

IBFTC  MATRIX  M94.XR7 

MATRIX  INVERSION  WITH  ACCOMPANYING  SOLUTION  OF  LINEAR  EQUATIONS 
SUBROUTINE  MAT  I NV I A.n.B.M.OE TERM) 


DIMENSION  IPIvOTId  I,  AI8.8  ).  BIS. 4  ), 

EQUIVALENCE  (IROW.JROW),  I  I COLUM « JCOLUM | . 


1NDEXI8.2  I ,  PIVOT (8  I 
(AMAX*  T,  SWAP) 


INITIALIZATION 


10  DETERM-1. 0 
15  DO  20  J-l.N 
20  I  PIVOT  I J ) -0 
30  DO  550  1*1. N 


SEARCH  FOR  PIVOT  ELEMENT 


C 

c 

c 


40  A MAX-0. 0 
95  01  105  J-l.N 

50  I*'  (IPlVOTIJi.u  60.  105,  60 
60  DO  100  <»  1  »N 

TO  Ir  I  I P I VOT ( X  I -1 )  80.  100,  740 

85  !Ioi2!S‘AMAX‘*M5,A,J*,t,»»  85»  10°*  100 
90  IC0LUM-X 
95  AMAX -A | J ,K ) 

100  CONTINUE 
105  CONTINUE 

110  1 P I  VO  T ( I COLUM I ■ I P 1 VOT ( IC0LUM|+1 


INTERCHANGE  ROWS  TO  PUT  pivot  element  on  diagonal 


130  IF  1 1  ROW- 1 COLUM | 

1 90  DETERM--DETERM 
150  00  200  L- 1 »N 
160  SwAP-AURO-.d 
170  A( IROW.L ) *AI IfOLUM.L) 
200  *  I ICOLUM.l ) -SmAP 
205  IF|M|  260.  260,  210 
210  DO  250  L-I*  m 


190.  260,  140 


lluU 
1110 
1120 
1130 
1190 
1150 
1160 
11  to 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
12pO 
1260 
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210  a( IRO^il  I  «dl ICOLUM.l) 

2*0  8 ( ICOLUM.LI»SwAP 
260  INDEXII.ll-IROW 
2 TO  INDEX  (  !  .2  >■  IC^L'IM 
110  PIV3TI I >  > A ( ICCLJM.ICOLUMJ 
320  DETERN«D£T£RM*P1V0I ( I  I 

DIVIDE  PIVOT  ROW  B»  PIVOT  ELEMENT 

130  All  COLUM.  I  COLllM  I  ■  1  •  0 
3A0  DO  1*0  L-l.N 

3*0  A ( I COLUM*  L I "A | I COLUM ,L  I /P I VOT III 
3**  I?(M|  380 .  380.  J60 
360  0)  370  L-l.M 

3  70  HI  ICOLUM.L  I -HI  ICOLUK.L  I /P l VOT I 1 1 

REDUCE  NON-PIVOT  ROWS 
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400  T*A I Ll » ICOLUMJ 

420  A  I L 1 . I COLUM I *0*0 
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550  CONTINUE 

INTERCHANGE  COLUMNS 
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660  SWAP*AIX» JROWI 

670  AIK.JROWI-AIK.JCOLUMI 

700  au.jcolumi-swap 

70*  CONTINUE 

710  CONTINUE 

740  RETURN 

7*0  END 
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SAMPLE  CALCULATION  no. 2 

ROTOR  IN  FLEXIBLE  “EDESTALS  AND  GYROSCOPIC  MOMENT 
27  4)0111101 

1.1 300000E  *0  7  1.0000000E*00 

7  1.000000E-03 

7. 6000006*01  3. 750000E+00  5.250000E+01  1.282000E+04  7.110000E*03 

2 . 1 AOOOOE+O I  A.C90n00E*00  4.880000E+0 1  O.OOOOOOE+OO  0.000000E*00 

3 .OOOOOOE  *01  5 . 310C00E.00  A.ttSOU^OE  +  Ol  0.000000E»00  O.OOOOOOE^OO 

1 .250900E»02  6.620000E+00  1.440000E.02  0.000000E.00  O.OOOOOOE*00 

1 .240 100E+02  2.2SOOOOEOO  3.02000CE+02  O.OOOOOOfcYOO  O.OOOOOOE+OO 

3.820')00E*02  5.910000E.00  7. 788000E.02  O.OOOOOOE«UO  0.000000£*00 

8 «  760  300E+02  9. 780000E.00  ?.  580000E.03  0.000000E*00  0.000000E*00 

7  #  7  70')00f  *02  7. 190000L03  8.380000E*02  0.000000E*00  O.OOOOOOE+OO 

S.70C3n0E»O2  3.57COOOEtOO  3.950000E.02  O.OOOOOOEOO  O.OOOOOOEyOO 

A . 5 10000E»02  7 . 190000E .00  3.859000E402  O.OOOOOOE^OO  O.OOOOOOEaOO 

2 • 350000E  *02  5.000000E.00  3.210000E+02  O.OOOOOOEyOO  O.OO0OOOE*OO 

3. OOOOOOE  *02  6.620000E.00  6.130000E+02  O.OOOOOOEtOO  0.000000£>00 

3.1 30000E+02  4.870000Et00  1.930000E+02  O.OOOOOOE+OO  0.000000E*00 

3.200000EAU1  2.1 3O00OE.00  6.800000E+01  0.000000E*00  O.OOOOOOE^OO 

3.200000E-01  3. 190000E.00  4.180000E>01  O.OOOOOOE+OO  O.OOOOOOE+OO 

3.6  72  000E+01  7 .420000E+00  2.000000E+01  O.OOOOOOEaOO  0.000000£*00 

2.821000E+01  8 .650000E+00  5.476000E*00  l.«20000E»02  9.700000E+01 

9.043000EO1  1.237500E*01  3.120000E-01  1.0»0000E*03  S. 450000002 

4.2970006*00  1.237500E*00  3.120000E-01  O.OOOOOOE+OO  0.000000£*00 

4.297000E*00  1 .237500E*00  3.120000E-01  0.000000E*00  0.0000006*00 

7 . 852000E*01  4.750000E*00  2.0000006*01  1.650000E+03  8.800000E*02 

3.4230006*01  5.3200006*00  1.092000E*02  0.0000006*00  0.0000006*00 

9.8200006*01  6. 190000E*00  8.540000E*02  0.0000006*00  0.0000006*00 

2.4310006*02  6.5100006*00  1.6020006*03  0.0000006*00  0.0000006*00 

3.3000006*02  6.210000E*00  8.360000E.02  O.OOOOOOE+OO  0.0000006*00 

2.B43000E+02  4.3300006*00  l.564000E*02  O.OOOOOOE*00  0.0000006*00 

6 . 3  70000E  *01  0.0000006*00  0.000000E*00  -0.0000006*00  0.0000006*00 

3  16  22  27 

1  1.0000000E*01  0.00000006*00 

10  1.00000006*01  1.00000006*01 
18  -1.00000006*01  0.00000006*00 

3.2000E+01  8.20006*05  2.20006*01  5.2000E*01  5.30006*05  2.20006*01 

5.20006*01  S.2000£*05  2.20006*01  5.20006*01  5.30006*05  2.20006*01 

4.40006*0 1  1.3000E*06  3.40006*01  4.4000E*01  9.00006*05  3.40006*01 

4.4000E+01  1.30006*06  3.4000E*0l  4.40006*01  9.00006*05  3.40006*01 

9.50006*01  1.70006*06  1.80006*01  9.5000E*01  1.3000E*06  1.80006*01 

9. 5000E*01  1 .70006*06  1.80006*01  V.5000E*01  1.30006*06  1.80006*01 

8.5000E-01  1 .9000E *06  2.70006*01  8.5000E*01  1.90006*06  2.7000E*01 

8.50006*01  1.90006*06  2.7000E*01  8.5000E*01  1.50006*04  2.70006*01 

1 .0000006*03  5. 100000E.03  2.000000E*03 

1.5420006*06  0.000000000  0.0000006*00 

3.3620006*06  0.0000006*00  0.0000006*00 

-8 .0109005*05  0.0o0000E*O9  0.0000006*00 

1.3027006*06  0.0O0000E*0O  0.000000E*00 

1.542000E>05  0.000000£*00  0.000000£*00 

1 .040  300E+06  O.OOOOOOEaOO  O.OOOOOOEaOO 

9. 720900E+0S  O.OOOOOOt.OO  O.OOOOOOEyOO 

1 . 302000E+06  O.OOOOOOCaOO  O.OOOOOOE^OO 

1.342000E«06  0.000000£*00  0.000000000 

3*  362000E+06  O.OOOOOOt.OO  O.OOOOOOEyOO 

•8.010000CY03  0.000000000  0.000000000 

1.302000E406  O.OOOOOOEaOO  O.OOOOOOEaOO 

1 .542000005  O.OOOOOOE.OO  0.000000E>00 

1.040000E«06  0. OOOOOOE *00  O.OOOOOOEaOO 

9.720000E409  O.OOOOOOE^OO  O.OOOOOOEyOO 

1.302000008  0.0O000QE«00  O.OOOOOOEaOO 

1.463000006  O.OOOOOOE4QO  O.OOOOOOE^OO 

3.296000E.06  O.OOOOOOEyOO  O.OOOOOOE^OO 


-1  ,0<.80C0t*C6 
1. 7o  1 0C96*Ce 
-1.1970006*39 
1.7660005*06 
1.1  760005*06 
1.7810005*06 
1.66  300^5  *06 
3.2760005*06 
-l.o<.50',05*06 
1.  7810005*06 
-1.1970005*06 
1.7660005*06 
1.1760005*06 
1.7810005*06 
1 . 070000E+06 
2.2600005*06 
-6.7000005*06 
9.1600005*08 
1.0700005*06 
6.8000005*06 
6.3200005*09 
9.1600005*06 
1.0700005*06 
2.2600005*06 
-6.7000005*05 
9.1600005*09 
1.0700005*09 
6.8000005*05 
6.3200005*05 
9.1600005*09 
9.8000005*09 
2. 1600005*06 
-7.0000005*06 
1.3  lOCOOL *06 
-8. *000005*06 
1.2000005*06 
7. 6000005^09 
1.3300005*06 
9.8000005*99 
2.1600005*06 
-7.000000t*09 
1.3300005*06 
-8. 7000006*06 
1.2000005*06 
7.f.0n00O**06 
1.3300005*06 
6 NO  05  5  i 


O.OOOOCOt  *00 
0.0000005*00 

o.oooooot*oo 
o.oooooot  *00 

0.000000t*30 
0.0000006*00 
0.0000005*00 
0.0000005*00 
0.0900005*00 
3.0000005*00 
0.0000005*00 
0.0000005*00 
0.0000005*00 
0.OO9900E*00 
0. 0033305*00 
0.0000005*00 
0.0000005*00 
0.0000005*00 
0.0000005*00 
0.0000005*0.3 
0.0000005*00 
0.0000005*00 
0.0000005*00 
0. 0000005.00 
0.0O0r0CL*30 
0.000000t*00 
0.0000005*00 
0.0000005*30 
0.0000005*00 
0.0000005*00 
0.0000005*00 
0.0000005*00 
0.0000005*00 
0.0000005*00 
0.0000005*00 
0.0O00005+00 
O.OOOOOOE*OC 
0.0000035*00 
.>•0000005*90 
0.0000005*00 
O.OOOOOOt *00 
O.OOOOOOt *00 
0.0000005*90 
0.0000005*30 
0.0600005*00 
0.0000005*0'' 
L5 


0.0000005*00 

0.0000005*00 

0.0000005*00 

0.0000005*00 

0.0000005*00 

0.0090005*00 

0.3000005*00 

0.0300005*00 

0.0000005*00 

0.0000005*00 

0.0000005*00 

0.000000t*00 

0.0000005*00 

0.0000005*00 

0.300000E+00 

0.0000005*00 

o.ooooooe*oo 

0.0030005*00 

0.0000005*00 

0.0000005*00 

0.0000005*00 

0.000000t*00 

o.oooooce*oo 

0.0000005*00 

0.0000005*00 

0.000000t*00 

0.000000t*00 

0.0000005*00 

0.0000005*00 

0.0000005*00 

O.OOOOOOE*OC 

0.0000005*00 

''.0090005*00 

0.0300006*00 

0.0000005*00 

0.0000005*00 

0.000000t*00 

0.3000035*30 
0.003030t*00 
0.030000t*00 
0.3000005*00 
0.000000t*00 
''.0000006*00 
0.0900005*90 
8. 0009006*00 
0.0000005*00 
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INPUT  FOR  PM0Q11: 


UNBALANCE  RESPONSE  OF  A  FLEXIBLE  ROTOR  IN  FLEXIBLE.  DAMPED  BEARINGS 

Card  1  Text  Col.  2-72 

Card  2  Text  Col.  2-72 

Card  3  (1015) 

_____  I.  NS.  Number  o<T  rotor  mast  stations  (<  80) 

____  2.  NB.  Number  of  bearings  (<  25) 

_____  3.  NU.  Number  of  unbalance  stations  (<  80) 

_____  4.  NC.  Number  of  coupling  stations  (<  20) 

_  5.  NPST.  0:  Rigid  Pedestal  1:  Flexible  Pedestal 

_  6.  NMOM.  0:  No  bearing  resistance  to  moment  1:  Moment  resistance  l 

_  7.  NGYR.  0:  No  gyroscopic  moment  1:  Gyroscopic  moment  calculation 

_  8.  NCAL.  1:  1st  type  of  bearing  data  input  >  2:  2nd  type  of  bearii 

data  input. 

______  9.  0:  no  diagnostic  l:  diagnostic  given 

_____  10.  0:  More  input  follows  1:  last  set  of  input 

Card  4  (1P4E15.7) 

_ 1. 

_ 2. 

IF  NGYR  -  1 
Card  (15,  IPE23.6) 

_ 1. 

_ 2. 


2 

E,  Youngs  modulus,  lbs/ln 

Scale  factor  in  simultaneous  equation  solutii 


NIT.  Nus&er  of  iterations  in  gyroscopic  mom.' 
Convergence  liartt  for  gyroscopic  moment  calo. 


97 


If  NGYR  -  0,  use  only  first  3  columns,  FORMAT  (1P3E14.6) 
If  NGYR  -  1,  use  all  5  columns,  FORMAT  (1P5E14.6) 

Give  one  card  for  each  rotor  station,  in  total  NS  cards 


Rotor 

Station 

(don't 

punch) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


Station  Mass 


Length  of 
shaft  sec¬ 
tion 


lbs . 


iH£ll 


Cross  sec¬ 
tional  Moment 
of  Inertia 
_ In* _ 


Polar  Masai 
Moment  of 
Inertia, 
lbs .in* 


Transverse  Mass 

Moment  of  Inertia 

lbs. in2 _ 


Rotor  stating.  vUh  Pouring  Support 
(1415) 

Give  NB  items 


Vnbglpnce 

(I5.1P2E15.7) 

Give  one  card  for  each  rotor  station  with  unbalance.  In  total  NO  cards 
.1.  Rotor  station  number 
.2.  X-component  of  unbalance,  os. Inch 
.3.  Y-component  of  unbalance,  os. inch 
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Rotor  Stations  with  Courlln 
(UI5) 


Applies  only  if  NC  i  0.  Give  NC  items. 


Pedestal  Data  for  Translatorv  Motion 
(IP6E12.4) 

Applies  only  when  NPST  •*  1.  Give  one  card  for  each  bearing,  in  total  NB  card 


Pedes .Mass 
x'dlrectlon 
lbs. 

Pedes .Stiffn 
x>direction 
lbs/in 

Pedes . Damping 
x'dlrectlon 
lbs. sec/in 

Pedes .Mass  1 
y-dlrectlt 
lbs. 

1  Pedes. Stiffn. 
y’dlrection 
lbs/in 

Pedes. D 
y-direc 
lbs.se 

I 

Pedestal  Data  for  Tilting 
(1P6E12.41 

Applies  only  when  NPST  »  1  and  NMGM  “  1.  Give  on  card  for  each  bearing. 
In  total  NB  cards. 


Mass  Mom. 
of  Inert. 
X’dlrectlon 
lbs . ln^ 

Angular 
Stiffn. 
x-direction 
lbs. in/rad 

Angular 

Damping 
X’dlrectlon 
lbs. in. sec/rad 

Mass  Horn. 

of  Inert, 
y-direction 
lbs. in* 

Angular 
Stiffn. 
y-dlrection 
lba . in/rad 

Angular 
Damping 
y-dlrec  tli 
lbs.ln.se 

Speed  Data  (1P3EU.6) 


1.  Initial  speed,  RPM 

2.  Final  speed,  RPM 

3.  Speed  increment,  RPM 


Bearing  Coefficients  for  Tranalatory  Motion 
(1P3E14.6) 


Give  8  cards  per  bearing, 
in  the  form:  K  -  K  +  K 

XX  XX, o 


in  total  8*NB  cards 
2 

,oi  +  R  .ui  ,  coC 
xx,  1  xx,  2 


xx 


Each  card  gives  one  coefficient 

2 

1  C  +C  ,oi  +  C  „ui  ,  etc. 
xx,o  xx, 1  xx, 2 


K 

XX 

cuC 

XX 

K 

xy 

cuC 

xy 

K 

yy 

cuC 

yy 

K 

yx 

coC 

yx 


K 

xx 

03C 

XX 

K 

xy 

cuC 

*y 

K 

yy 

<uC 

yy 

K 

yx 

uiC 

yx 


100 


Bearing  Coefficients  for  Tilting 
(1P3E14.6) 

Applies  only  when  NMCH  »  1.  Give  8  cards  per  bearing,  in  total  8*NB  cards 

Each  card  gives  one  coefficient  in  the  fora: 

2  2 

M  ■  M  +  M  .a:  +  M  „a>  ,  tuD  ■  D  +  D  ,co  +  D  ,a>  ,  etc. 
xx  xx  ,o  xx,  1  xx  ,2  xr  xx,o  xx,l  xx,2 
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Type  2  Bearing  Data,  NCAL  >  2. 

Repeat  the  following  input  aa  aany  time*  aa  given  by  PCA'.,. 
Si-eed  Data  (1P4E14.6) 

_ _ .  Speed,  RPM 


Bearing  Coefficients  for  Tranalatorv  Motion 
(1P4E14.6) 


K 

JLD  C 

,K  ,CD C 

■  XX  XX 

xy  xy 

K 

.coC 

,K  ,<uC 

yy  yy 

yx  yx 

K 

,uC  , 

K  ,cuC 

‘  XX 

XX 

xy  xy 

K 

,u£  , 

K  ,cuC 

•  y> 

yy 

yx  yx 

Bearing  Coefficients  for  Tilting 
(1P4E14.6) 

Appliea  only  when  NMGH-1.  Give  2  carda  per  bearing,  In  total  2*NB  carda 


M 

cuO 

,M 

,<uD 

XX 

XX 

xy 

xy 

M 

,cuD 

,M 

,wD 

yy 

yy 

yx 

yx 

M 

0)0 

,M 

,(o0 

XX 

XX 

xy 

xy 

M 

,u)0 

,M 

,o>D. 

yy 

yy 

yx 

yx 
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APPENDIX  b 
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103 


7^ 


$ J08  0.5.5000  65-3*1  .FREEMAN. SESCD 

sexecute  i a jos 

$  I  8  JOB  rSfVAN  W4P  .. 

SIHFTC  SOTOR2  «9*  .xR3-...SoECX 

COMMON  XC  <  9  •  70  )  «  X  S  (  9  •  30  I  •  VC  I  9 . 30  >  .YS  (  9  •  30 1  tDMXA  (  30  I  .DMVAOO)  .DV 
1 X A  I  30  I • DV X3 ( 30 • .OVXLX  301 .OVYAI 30 ) .DVYb I  30 1 .OVYC ( 30  I .DVYD( 301 .OVXC( 

2  30) .RM( 30  I .Rl( 30) *HS( 30) .RIPt 30) *AN( 30) *8NI 301 »DN( 30) * 

3  L8I 10) .dKXX ( 10)  .BCXXI 10)  .8KXY 110) .BCXYI 10) .BXYYI 101 »8C YY ( 

430) .BXYXI 10> .BCYXt 10> .PRMSI 10) »8RFRI 10) .FRQ1 1 100 1 .FREUI liO) *  PMX( 
510) .PMYtlOl .PKXt 10) »®CX(  10) .PKYI 10) .PCY( 101  *  CFM(8.ai .ENH9) 

6  «CLNR ( 8  I *  82N 1 (8.8) *d2NP ( 8»8I »82N2 (3*8) .CF9.FRW.0T .OR.OE .EN0Y1 *821 • 
7B2P.822 

2 00  REAO  (5.1001 

RE AO  I  5. 105  INS.N8.NFR ,NCAL .NPST.INP 
RE AO  15.103)  YM 
WRITE  (6.100) 

WRITE  ( 6.  10d  ) 

WRITE  (6.1 10 INS.NB.NFR.NCAU  .NPST.INP 
WRITE  (6.1 0* ) YM 
WRITE  (6.1111 
WRITE  (6.116) 

00  201  J*  1  .NS 

RE AO  (5.103) RM(J).Rl(J).RS(J1.R1PIJ) 

201  WRITE  (6.114IJ.RMI Jl.RLI J).RSIJI.RIP(J) 

READ  (5. 118)  (L3(J)  .  J*  1  >N8 I 

WRITE  (6.120) 

WRITE  (6*116)(l8(JI«J*1»NB) 

IF(NPST)  208.228.208 

208  WRITE  (6.128) 

WRITE  (6.127) 

DO  209  J» 1 »N8 

READ  (5.125)  PMX  (  J )  .PAX  l  J I  .PCX!  J I  .PMY  I J )  .PAY I  J)  .PCY  (  J ) 

KST  *LB  (  J ) 

209  WRITE  (6.126IKST .PMX (J).PKX(JI.PCX(J) »PMY (  J  I  <PA  Y  ( J  ) .PCY ( J) 

228  READ  (5. 103) (FRCl(J) .JM.NFR) 

WRITE  (6.101) 

WRITE  (6.103)  (FRQKJl  ,J*1,NFRJ 
C  CONVERT  INPUT  UNITS 

250  A  AS» 1000.0 
C‘*l«386.069*AwS 
RS(NS)»RS(  1  I 

00  251  J-l.NS 
RM( J)»RM( JI/CF1 
RIP(J)«RIP( JI/CFl 
STF-YM/AMS*RS( J) 

AN(J).RU JI/STF 
8N(J)*RL( J)/2.0*AN(J) 

251  0N( J)»RL(J)/3.0*BN(J) 

229  READ  ( 5 . 103  I SPST .SPFN .SP INC 
WRITE  (6.1*1) 

WRITE  (6. 103ISPST.SPFN.SPINC 

00  20*  J-1.N3 

AST*LB(J) 

READ  (5.103 )RAXX( J I .BCXXI J)  .RAX Y ( J I ,BCXY ( J ) 

READ  (5.103) 0AYYIJ) .RCYY(J)  .BAYXIJI.BCYXIJ) 

WRITE  (6.121) KST 
WRITE  (6.122) 

WRITE  (6.103) RAXX ( J )  .8CXX ( J )  .  3<XY ( J  )  .BCXY(J) 

WRITE  (6.1231 

WRITE  I6.103IRAYY(J1.BCYY(J)  .SKYXI J) .BCYXI J) 

C 1 *BCXX ( J I ♦OCVY ( J I 

Cl«(BAXX( JI*BCYY( J)*BAYV( JI»BCXXI JI-BAXYI J )  *8CYX I J l-BAYX( JI»BCXYIJ 
II )/Cl 


O010 
oQ2U 
uO  JO 
uQ*0 
0050 
0060 
00  70 
0080 
0120 
0130 
01*0 
0160 
0160 
0170 
0180 
0190 
0200 
0210 
0220 
0230 
0240 
0250 
0260 
0270 
0260 
0290 
0300 
0310 
0320 
0330 
0340 
0350 
0360 
0370 
0380 
0390 
0400 
0410 
0420 
0430 
0440 
0450 
0460 
0470 
0480 
0490 
0500 
0510 
0520 
0530 
0540 
0550 
0560 
0570 
0580 
0590 
0600 
0610 
0620 


1  __ _ vJi 


V*. 


i^,C2>  206. 205.205 

205  C2*S3?T(C2| 

C1«C2 

GO  TO  ?o? 

206  C1.-1.0 

207  WHITE  (6.1241 

"1,lE  I6.10JIC2.C1.C4 

pH*S(  j|.C4 
0/XA( j).c J 
2 04  8-»PR(J|«C3 

*dR«0 

iu  <sr«o 
nbR«mbr»1 

ERW»FRQ1(\«fb, 

00  213  Js|.N3 
C 1 »Dyx  A ( J | 

,F,CH  213.213.210 

211  Hift'**1  »'•«*•»» 

*ER*WFR-l 

212  K$T«  J 

213  CONTINUE 

*F.< <ST  »  216.216.213 
213  0VXA(<sri«-J,0 

216  FREQ(»BR)  sFRm 

217  HI™" 

21S  VFR«NFR«| 

KST-0 

00  222  J»UNd - - - - 

C1.DVKAIJI 

„A  ”2*«2.220 
II?  »»•«»•»»! 

KST«J 

222  CONTINUE 

IFMC3TI  ??4.??4.J»1 
221  NBR-M9RM  '”*'**' 

P«CQ|N3R|.C2 
OVXAUSTl  —  l.o 
00  TO  219 

224  N*ftl»/BR 
S’CAl-SPST 

225  WHITE  (6.146KPCAL 

A  <5PR«0«  1  0*  7 1 974#j(>r »l 

C  2*ANSPR*ANSPR 

**ITE  16.1481 
DO  226  J« | (nb 
TST«tN( J| 

Cl.PRMSl J|/C2 
CJ«8RFR|j|«SPfAL 

226  uoIIS  “*»*7KSr.C3,Cl 

WRITE  16* 147 | 

<0C»0 

VFR.| 

300  pRW»FRE3l.viFR*C0UE,KV  °CMCNDtNr  pa«AMETE*!» 


U6J0 

0640 

0650 

46b0 

06/0 

06o0 

0690 

0700 

0710 

0720 

0730 

0740 

0750 

0760 

0770 

0780 

0790 

0800 

0810 

0020 

0830 

0840 

0850 

0860 

0870 

0860 

0890 

0900 

0910 

0920 

0930 

0940 

0950 

0960 

0970 

0980 

0990 

1000 

1010 

1020 

10J0 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

1140 

1130 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

124U 

1230 


A\<;p»ANSPR»F.7'.v 

A\*P2*  A\S°»  A,\<-P  1260 

ro  vi  j* i  .no  12 '0 

*^TF  I?  I  3 1  J  |  »  VNCP?  1200 

n^x.-.i  ji  »ftf  1270 

r'vrA  (  j  i  .<;  rr  1300 

STe»1V|  JJ.A.NSt'?  1310 

!'VX  A  (  J  I  »o  TF  1320 

;vya(jixstf  liiu 

-VX3(JI*0.0  1340 

f'VXC(J)*0.0  133  0 

'vx:i JI.J.O  1360 

DVYB ( J ) *0.0  13/0 

DVYCIJMO.O  1 3oO 

vi  dvyoijvo.o  1 3vo 

302  00  111  JM.NH  1400 

<  ST  »UQ ( J I  1410 

CF1<=?<XX(J!  1420 

CF  K  =  3CXX  (  J  l»FR*  1430 

Ccn  =  .3XXY(JI  1440 

Cclt’*',CXY  (  J  !  »  r  R  2  1490 

C'2<*«<YY(J)  1460 

CC2C*BCYY<  J)*FRA'  1470 

C:20*dKYX(JI  1400 

C-'2E*0CYX»  14*0 

IF(NPST»  303. 306. 303  1500 

303  CrtM=.p<X(  J)-°"X{  J|/3,J6.06?»ANSP2  lsl° 

CF1N*PCX(J)*ANS°  1520 

CFiA*CclK>CFlv  1630 

CF1H*CF10CF1N  1540 

CF2M«P<y(  JI-Pmyj  J  1/386. 069»ANSP2  1550 

CF2M*PCY»  J|*AX!SO  1560 

CF2AiCF2<*CF2v  I3/O 

CF2t3*CF2C+CF2N  l96u 

GO  TO  307  1 6 VO 

306  CF’  •  CF 1 <  1600 

CF  2*CF 1C  1610 

CF3«Cc10  1620 

CFa.CFIF  1630 

CF6«Cc20  1640 

CF6*CF2F  1690 

CF7*CF2<  1660 

CFa«CF2C  1670 

GO  TC  308  1660 

307  CF4*CP2A»CF2A*CF2B*CF2B  16V0 

CF;«(CF2A#C<  jo+C'2.i*CF2t  I/CF4  l700 

CF2«(CF2A»CF2r-CF2o*CF2D)/CF4  1710 

CF3«(CF2A*CF2m*C^2o«CF2N)/CF4  1720 

CF4*(CF2A«CF2»!-cF2c»CF2V|/(;fa  *  1730 

CF8«CF1 A-CF  1  »f F  1  )*CF ?»CF 1 r  1740 

CF6«CF13-rF2*cFtO-Crj»CMf"  1790 

CF7.-CF3»CFio*CF4»CFlE  1760 

CF8*-CFA»CF10-CF3»CF1C  1770 

CF2N»CF6*CCI>  +  Cr0*CF6  •  1760 

CF2A*iCF3«CF  1m*cF6*CHM /C*- 2  s  1 /VO 

CF2o»(CF6*CF1  J  F6*C.  1'3)/CF  2iN  1600 

CF2m.(Cf^#Cc /♦CF6»(.F6»/CF2S  1*10 

CF2N»(CF6*tFb-CF6»'.F  73/Cf  ?\  1**0 

C«-1A.-CF1  »CF2a*ci  2»CF2i  1*»30 

Ceia«CFl*CF2H*CF2»Cr2A  1640 

C_lY»CFl-CFi»^F2w.>.F2»CF?N  1*50 

C‘1N«C14-CF 1»CF2N-;f2»CF2v  i860 

C-!*CFl<.CF2S.CF|r.-r2n,rr,,„cri^cril#{.F1H  jd7o 


CF2»CF1<«CF2o»CF1CACF2A-CF1D*CF1BaCF1E*CFIA  i«*u 

CF  3«CF J«»CF2m-CF 1C»CF2N  +  CF1D»CF1M-CF1£«CF1N  190  0 

CF^.CF 1K»CF2NaCF 1C»CF2**CF JO»CF INfCF 1 E*CF IM  1910 

CF5«CF2D»CF2A-CF2£»CF2d*CF2K»CFlA*CF2C*CFlt>  192^ 

CF6*CF2D»CF2d*CF2£»CF2A-CF2<»CFlD^CF2C»CF  1A  Oju 

CF  7'tF20*CF2H-CF2£*CF2M»CF2A#CFliH-CF2C*CFlN  lq<*0 

CF8»CF2D*CF2N*CF2£,CF2M+CF2<»CF1N+CF2C*CF1M  1950 

308  0VXA(<5T|«DVXAIKSTI-CF1/AMS  1960 

DVXBHCSTI-CF2/ANS  19?3 

0VXC(<ST»«CFJ/AMS  „  _ 1,00 

0VXDI<ST|«CF4/A.mS  1^90 

DVYAI<ST|«0VYAKSTI-CF7/AMS  2000 

OVYB(XSTl«CF8/AMS  2010 

DVVC(<ST)«CF5/AM$  2020 

0VY0(KSTl«CF6/AMS  2030 

311  CONTINUE  2u'»C 

ROTOR  CALCULATION  '  2050 

00  428  I  - 1  •«  2060 

KS  T* I  2070 

BMXC»0.0  2080 

BMXS-0.0  2090 

BMYC-0.0  2100 

8MY5-0.0  2110 

VXC»0.0  2120 

VXS-0.0  2130 

VYC»0.0  .  »  2140 

VVS«0.0  I  2150 

xcn»n-o*o  2160 

XS 1 1 • 1 1 ■0«0  2T7J 

vcii.n-o.o  2100 

YS 1 1 • 1 ) »0»0  2190 

DXC-0»0  2200 

OXS-O.O  2210 

D'C-0.0  2220 

DYS-0.0  7270 

0)  TO  1 407  #408  »4Q9  »4 10  *4 1 2 1 4 1 3  «41 4»41 5  )  .XST  2240 

407  OXC-1.0  2250 

GO  TO  418  2260 

408  DXS«1*0  2270 

GO  TO  418  2280 

409  DYC-l.Q  “2270 

GO  TO  418  2300 

410  DYSM.O  2310 

GO  TO  418  2320 

412  XCIltlXl.O  2330 

GO  TO  418  2340 

411  XSI6.1I-1.0  2550 

GO  TO  418  2360 

414  YC(7tll*|.0  2370 

GO  TO  418  2360 

415  YSI8.11-1.0  2393 

418  00  424  J-l.NS  2400 

'  BMXC-BNXCaONXAI JI*OXC  *  2¥T 0 

BHXS»BNXSaDMXaUI«OXS  2420 

BMVC»BMVCA0MYAtJl*0YC  2*30 

8MYS-3NYS*0MY*U1#0YS  i  2440 

ClaDVXA (JI»XC(l»JI-OVXBIJI*XSII*Jt-OVXC(JI*YC(I • J I -OVXDl J ) *YSI I  •  J I  2450 
C2*OVXB(JI*XC(l » J I *OVX A 1 J| *XS ( I » JlaOVXDI J l*VCI I » J I -OVXC ( J I *YS ( I  « J)  2460 
CJ—OVYCI  JI*XCn»J*-3V.^I  JI«XS<  ItJlAOVYAC  JI*YCI  I.Jl-OVYBl  '  7*70 

II  2480 

C4«0VY0( JI4XCI 1 .JI-OVYCI JI»X.M  | «J|60VYBf J)«VC| 1 1 J •♦DVYAl Jl »YS ( 1 t Jl  2490 
VXC-VXCK1  2500 

VXS«VXS*C 2  7510 


/YC=VYC*C3 

VvS-VYS*C4 

42?  I  c  (  N'-J  I  424.4  24.423 

421  XC( 1 ,J*1 I *XC( t .J)*RU J)*DXC  ♦BN(J)*8MXC  ♦DN!J)*VXC 

XS( I ,  J  +  1 ) *XS( ! »J)*RL(JI*DXS  ♦BN ( J  )  *B*XS  ♦ON(J)*VXS 

YC  (  I  ,  J»1  )  *  YC  (  I  »J)*7L(J)*DYC  ♦BNIJ)*BMYC  ♦DN(J)*VYC 

Y i( I »  J*1 ) *  v  S  ( I » J ) *RL ( J )*DYS  ♦BN(j)*BMYS  ♦DN!J)*VYS 

2XC*DXC»AN!J) *BMXC+9N(  Jl *VXC 
R<S»DXS*AN< J>*BMXS*"N( J)*VXS 
DYC*OYC+AN( J)*tiMYC*wN( J»*VYC 
:y.S«DYS  +  AN(  j  )*BMYS*UN(  J)  *V7S 
B'XC*QWXC»RL(  J)*VXC 
8MXS*8MXS*RU JI*VXS 
dMYC*9MYC+3L(JI*VYC 
BMYS->BMY$*RL(  J)*VYS 
424  CONTI  NUE 

C Fm  ( 1 . 1  I *8MXC 
CFM( 2 • I ) *8MxS 
CFm( 1,1)  *9MYC 
CFm ( 4  ,  I  )*9vyS 
CF V ( 5  ,  I  )  *VXC 
CFMI6.I  )’VXS 
CFM( 7, I  I *VYC 
CFM(B,II*VYS 
428  CONTINUE 
CALL  EOS 

5 06  VFR*WF7*1 

IF(MFR-NFRl)  300  »  30G  .507 
AOVAnCE  speed 
SOT  SPCAL*St>CAL  +  SOINC 

IF(SOFN-SPCAL)  511*225*228 
PROGRAM  £N0 

511  <DC*<0C*1 

IF(NCAL-KDC)  S10.S10.229 
S10  IF(INP)  S09.200.S09 
S09  STOP 

ICO  FORMATI49HO  •  ) 

101  FORMATI 1 THOFRFOuENCY  RATIOS) 

1 0  3  FORMA  T  (  4  (  1XE13.6)  I 

104  format ( 16HOYODNGS  MODULUS*  *  £14.71 
lrs  F0RMATI7IS) 

108  FORMA T ( 58M0  STATIONS  NO. BROS.  NO.FHtO  NO. SPEED  PtD.FLEX  I 

1 NP  I 

110  FORMA  T(4Xl4.6ISXl4)  ) 

111  FORMAT! 14H0  ROTOR  DATA) 

1  14  F  1RMAT  (  4X  I  4 . 7x£  l  3.6 . 1 XE 1  3.6 . 1 XE  1  3. ft  .  1  Xt  1  3.6  I 

116  F  >RM A T 1 6  7h  STATION  NO.  MASS  LENGTH  CKO. SECT • I NER 

IT  (IP-ITI) 

117  F)RMAT(4X14.7XE13.E.1XE13.6) 

118  FORMAT! 10(1X141  ) 

120  FuRMAT! 18H0BEARING  STATIONS  ) 

121  FORMAT (24So  H'ARING  AT  STATION  NO.. 13 1 

122  FORMAT! S2H0  XXX  «CXX  XXY  XCXY) 

!?3  FORMAT  (  S2H0  <YY  V.CYY  XYX  WCYX) 

124  F0RMATI43M0  INST.FRE3.RT  MASS*! FREO I **2  i*F  l  GHT*  ( W)  *»2  I 

125  F  ORMA  T I  6 ( lxEll.41  ) 

l’S  FORMAT! IX I4.4XE1 1 .4. 1XEI 1 .4, 1 XE 1 1 .4 • 1 XE 1 1 .4 , 1XE 1 1 .4, 1XE 1 1 .4 ) 

1 2 T  F ORMA T  (  ?6M0  BRG.  '.T.  MASS.X-OIR.  XX  CX  MASS.Y-DIR 

1  .  <  Y  C  Y I 

128  S0RMATI1SH0  P"OFr,TAc  OATA) 

141  F0RMATI42M0INJTIAL  SPEED  FINAL  SPEED  SPEED  INCR.) 

14F  FORMAT!  13HlR(.ro<  SPEED*,  E13.6.3HHPM) 

147  r  ORMA  f  t  6“  ‘■*0  CR£  3.R AT  »  DETERMINANT  REIDLT)  IM(DET) 


252u 
25.JJ 
2540 
25SO 
,  2560 
2570 
2580 
2590 
2600 
2610 
2620 
2630 
1  2640 
2650 
2660 
2670 
26o0 
2690 
2700 
2710 
2  720 
2730 
2  740 
2750 
2  760 
2770 
2  780 
2  790 
2800 
2810 
2820 
2a30 
2840 
2850 
2860 
2870 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3010 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
1110 
3120 
3130 
3140 


I  ENERGY  I  31^1 

148  ECRmaT  14?ho  890.  STATION  I NST ,E REG «RP“  INST  ,»E IGhT  )  3lt>( 

ENO  317C 

KRETC-SEOG  M94  »XR7 

SURRC  JT I NE  EOS  J«-U 

'OWWi-,N  XC  (9.30!  *XSI  9.301  «YC(  9,30!  *S  Srj.30I.OMX4(  301  .OMYAI  30)  «0V  uOpC 

1X4(301. DVXjl 30 ) «DVXO( 30) .OVT 4  I  30 ) .DVYB ( 30 )  .DVYC  f  30  I  .OVYDI 30!  .OVXCt  -Cfc 

2301  »3M(  30!  «RL  (  301  .MSI  30!  ,  R  1  P  (  301  .ANl  30!  »8NC  301  »0N( 30  1  .  -0  n, 

3  U8I  10' .8KX  X ( 101 .H<"xx  I  10  I .3<XY ( 101  .BCXTI  lul .B4YYC lu)  .8CYYI  --B4 

4  30)  .BKYXI  10*  .ACYXI  lOl.PRMSI  101.8HFRI  lOl.ERClUOOI.MEOI  UC*.  PMX(  00^ 

510) . PM Y I  101 • P y x  f 101 .PCX ( 101 .PKY I  10! »PCY( 101  *  LFMI8.8 1 .ENT I  9  I  -loU 

6 .CL NR  I  8  1 .32 Nl  (  j  •  8  1 . 82NP I  8 .a  I .82N2 I  8 .8 ! »CF9  .FRtg.DT  .UR.oE .ENSYl .821 .  -11. 

782P.822  012' 

PRN3-1.0  013C 

DO  850  T»1.8  0140 

ENT ( I  I *0.0  0150 

M  850  J*  1 .8  Olfec 

850  R’Nllt.Jf-CFMII.JI  017L 

CILl  mat INVI82NI .8.ENT  .1,01)  oldC 

□  r  »S3RT ( PRN340T I  01VC 

v).7  0200 

D1«CEM( t . 1 |  0210 

OE-CFMIl ,21  0220 

831  DO  8  51  1*1, MO  0230 

00  81?  J* 1 »V0  0240 

B2N1I  I,J!-CPM»|,J1  0250 

92NPC I,J!»CFM|1,J|  0260 

832  B2N2( I ,JI *CEM| I.J)  0270 

833  82NP(  I.mO|*CFmU,mo*1I  0280 

CALL  M4TINVI8?Nl.M.-),ENTf  1.821  I  02V0 

CALL  M4T|NV(B2NP.M0,ENT.l,tt2P|  0300 

I F  ( MO-7 1  853.834.854  0310 

854  ENGYlO. 1413927  *B2P  0320 

853  MO-MO-1  0330 

CALL  M4T1NVI82N2.MJ, ENT, 1.8221  0340 

B2 1 *821/822  035C 

B2P.82P/822  0360 

C 1 "821 •0R-H2P»DE  0370 

0E.821*DF  ♦82P*0R  0380 

0R«C1  0390 

“D«Mf)-l  0400 

I F (MO- 3  I  856.851  .851  0410 

836  WRITE  I6.880IERM.0T .DR.DE.ENGY1  0420 

RETURN  0430 

880  FORMA  T I  3 (  1XE13.6!)  m44u 

END  0430 

SIBETC  SMATIN  M94.XR7 

SUBROUTINE  MATINV  IA.N.6.M. DETER!  UUIO 

C  MATRIX  INVERSION  WITH  ACCOMPANYING  SOLUTION  Of  LINEAR  EQUATIONS  0020 

DIMENSION  A ( 8 « A  1.819  I.IPIV0I9  ).PtV0T(9  I  0030 

DETER  *1.0  0030 

00  17  J«!.N  0060 

KC«0  0070 

KR.Q  0080 

00  14  1*1. N  0090 

IeIA{ I.J! I  11,12.11  0100 

11  KC«1  0110 

I *1 A| J, | t I  13,14.13  4120 

13  KN»1  0130 

14  CJNTINUE  0140 

I-'UCI  15*16.13  0160 

15  lEIXRI  17.16.17  0160 

16  OcTER.O.O  .  017C 
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00  TO  600 
IT  CONTINUE 
00  20  J*1  ,  N 

20  IPIVOI J)»0 

00  550  I » 1  »N 


SEARCH  FOR  PIvOT  ELEMENT 


4MAX*0»0 
DO  105  J»1,N 

nZ  1(rPrv0‘-'»-l»  60,105.60 
60  00  100  <*1,N 

If  UPlVOMCI  -n  80,  100.  600 
85  85.100.100 

ICOLU  «< 

AM4x»AI J.K) 

100  CONTINUE 
105  CONTINUE 

IP!VOI  ICOLUI-.IPJ  V0(  IC0LU>*1 

c  M“s  10  «<r  »i«r  <*  0IW0B>L 

ua  MjM'SSefS*”  ««• 

DO  2 00  L*1»N 
A  MAX  *  A  I  1  ROW  1 1.  I 
A <  t  RCW  »L I  *  A  ( ICOLU, LI 
<00  A ( ICOLU »L ) *AMAX 

aiax*8( iaowi 

o( IROW1-BI ICOLU! 

0 I ICOLU) *  AMAx 

260  PIVOT! I |*A( ICOLU. 1COLUI 
D" TER  *0ETER»PIV0T i i j 

0 1 V I OE  PIVOT  ROM  BY  P,voT  ELEMENT 

A<  ICOLU.  IC0LUM1.0 
DO  350  L*1,N 

Hl<!r21LU,L,1,A"col-u’‘-i/Pivom) 
d< ICOLU I *d ( ICOLUI/PlvOTt I  I 


REDUCE  NON-PIvOT  Runs 


38J  00  550  L  1  s  1  <N 

Ann  iMixl;!C?LU'  i,0°*  ^<50 

Ann  AMA X x A | L 1 .ICOLU) 

4  IH  ,  I  COL  J  I  «0«0 
00  450  L*1.N 

440  A,LJ  •<->«A(Ll.U-AI  ICCiLU.L)*A«AX 

550  C(Snt!nue‘'1''O<IC0LU,*A*'M* 
poo  RETURN 
END 

END  Of  FILE 
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INPUT  FOR  PNOQ17; 


STABILITY  OF  A  FLEXIBLE  ROTOR  IN  FLUID  FILM  BEARINGS 

Card  1  Text  Col.  2-49 
Card  2  (615) 

_  1.  NS  Number  of  rotor  mass  stations  (<  30) 

_  2.  NB  Number  of  bearings  (<  10) 

_  3.  NFR  Number  of  frequency  ratios  (<  100) 

_  4.  NCAL  Number  of  rotor  speeds 

_  5.  NPST  0:  Rigid  Pedestal  1:  Flexible  pedestal 

_  6.  0:  more  input  follows  1:  last  set  of  input 

Card  3  (1XE13.6) 

2 

1 .  E,  Youngs  modulus,  lbs/in 


ROTOR  DATA 

FORMAT  (4(1XE13.6)  ) 

Give  one  card  for  each  rotor  station,  in  total  NS  cards 


Rotor 
Station 
(don  't 
punch) 


Station  Mass 


lbs. 


Length  of  Shaft  Sec 


*nch. 


Cross  sectional  Moment] 
of  Inertia,  in^ 


Polar-Tra  nsvei 
Mass  Moment  oi 
Inertia  - 
lbs. in 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 
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Rotor  Stations  with  Bearing  Support 
(10(1X14)  ) 

Give  NB  Items 


Pedestal  Data 
(6(1XE11.4)  ) 

Applies  only  when  NPST  ■  1.  Give  one  card  for  each  bearing,  in  total  NB  cards. 


Pedes. Mass 
x-directlon 

lbs . 

Pedes . Stiffn. 
x-direction 

lbs /in 

Pedes. Damping 
x-direction 

lbs . sec/in 

Pedes. Mass, 
y-direction 

lbs. 

Pedes.  Stiffn. 
y-direction 

lbs . /in 

Pedes . Damping 
y-direction 

lbs . sec/ in 

FREQUENCY  RATIO  VALUES 
(4(1XE13.6)  ) 

Applies  only  when  NFR  >  1.  List  as  many  values  as  given  by  NFR,  4  values  per  card. 
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BEARING  DATA 

Repeat  the  following  input  as  many  times  as  given  by  NCAL 
Speed  Data  (4(1XE13.6)  ) 

_ 1.  Initial  Speed,  RPM 

_ 2.  Final  Speed,  RPM 

_ 3.  Speed  Increment,  RPM 


Bearing  Coefficients 
(4(1XE13 . 6)  ) 

Give  2  cards  per  bearing  with  4  coefficients  per  card,  In  total  2  NB  Cards 


K 

,u>C 

,K 

,ojC 

XX 

XX 

xy 

xi 

K 

,GUC 

,K 

,'i£ 

y> 

yy 

yy 

yx 

K 

,tuC 

,K 

,uC 

XX 

XX 

xy 

y> 

K 

,cjuC 

,K 

,tX 

yy 

yy 

yx 

y* 
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1  ABSTRACT 

This  report  is  a  manual  for  using  the  tvo  congmter  programs: 

1.  "Unbalance  Response  of  a  Rotor  in  Fluid  Film  Bearings." 

2.  The  Stability  of  a  Rotor  In  Fluid  Film  bearings." 

The  report  gives  the  analysis  on  uhich  the  programs  are  based,  and  the 
instructions  for  preparing  the  computer  input  and  for  interpreting  the  computer 
output. 
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